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PREFACE 


This  is  the  second  of  two  Workshops  on  turbulence  held  by  the  Panel  in  the  1986—87  time  period.  The  first  Workshop 
concentrated  on  the  measurement  of  turbulence  and  methods  of  data  collection.  The  theme  of  the  second  Workshop 
concerned  the  criteria,  analysis  methods  and  regulations  involved  in  the  design  and  certification  of  aircraft  for  turbulence. 
Taken  together,  these  two  Workshops  will  provide  invaluable  guidance  in  the  formulation  of  an  AGARD  Manual  on 
turbulence,  scheduled  for  publication  in  late  1988  or  early  1 989. 

The  authors  ate  to  be  congratulated  for  the  interesting  and  valuable  presentations  at  the  Workshop,  and  the  appreciation 
of  the  Panel  is  hereby  extended  to  them.  A  special  note  of  thanks  is  offered  to  the  FAA  and  the  AGARD  Flight  Mechanics 
Panel  for  the  contribution!:  sponsored  by  those  organizations. 
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Le  present  rapport  rend  compte  des  travaux  de  la  deuxiime  reunion  de  travail  organisee  par  le  Panel,  sur  le  th&me  de  la 
turbulence  pendant  la  periode  1986—  1 987.  La  premifcrr  reunion  fut  consacree  aux  mesures  des  turbulences  et  les  procedures 
de  recueil  des  donnees.  La  deuxieme  reunion  de  travail  conceme  les  critires,  les  methodes  d’analyse  et  les  rfcglements 
intervenant  dans  l’etude  et  1'homologation  des  aeronefs  du  point  de  vue  de  leur  aptitude  au  vol  dans  les  turbulences 
atniosphfriques. 

Les  travaux  de  ces  deux  reunions  de  travail  devraient  deboucher  sur  des  directives  qui  seront  d’une  grande  utilite  lors  de 
('elaboration  du  Manuel  AGARD  de  la  turbulence,  dont  la  parution  est  prevue  fin  1 988  debut  1 989. 

Nous  tenons  it  feliciter  les  auteurs  de  l'interet  et  de  la  valeur  des  exposes  qui  ont  etc  presentes  lors  de  la  reunion,  et  de  la 
part  du  Panel,  nous  leur  presentons  nos  vifs  remerciements.  Nous  ne  saurions  conclure  ce  resume  sans  exprimer  notre 
reconnaissance  pour  les  contributions  effectuees  sous  I’egide  de  la  FAA  et  du  Panel  AGARD  de  la  M6canique  du  vol. 
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by 


Terence  J.Bames 
Federal  Aviation  Administration 
ANM-105N 
P.O.  Box  C-d8966 
Seattle,  WA  98 168,  USA 


INTRODUCTION 


This  paper  presents  an  FAA  overview  of  the  gust  criteria  and  analysis  methods 
used  In  the  various  types  of  flight  vehicle  certified  under  the  FAR's 
(Reference  1).  The  current  criteria  for  small  airplanes,  transports,  and 
rotorcraft  are  presented,  and  the  status  of  proposed  criteria  for  the  tilt 
rotor  and  aerospace  plane  are  discussed.  The  amount  of  discussion  on  each 
class  of  vehicle  depends  on  the  significance  of  gust  loads  as  design  loads, 
and  the  importance  of  vehicle  flexibility.  Transport  airplane  gust  criteria 
development,  usage  and  problems  are  discussed  in  some  detail.  Analysis 
methods  used  oy  U.S.  Industry  are  covered  in  a  separate  paper. 

Units  used  throughout  this  paper  (Standard  English)  are  those  accepted  world 
wide  for  certification  to  FAR-2b  and  JAR-2S  standards. 

PHILOSOPHY  OF  REGULATIONS 


Structural  design  limit  load  criteria  are  chosen  such  that  there  is  sufficient 
structural  margin  around  normal  operating  loads  that,  in  combination  with  the 
factor  of  safety,  the  probability  of  catastrophic  failure  is  acceptably  low. 

The  structural  criteria  used  to  establish  the  required  strength  levels  have 
changed  over  the  years  to  reflect  changes  in  airplane  configurations, 
knowledge  of  the  atmosphere  and  ability  to  analyze. 

To  ease  the  analysis  burden,  criteria  are  presented  in  the  simplest  form 
consistent  with  obtaining  acceptable  safety  levels.  However,  it  is  necessary 
to  recognize  the  impact  of  the  simplifications  on  a  particular  vehicle.  A 
typical  example  is  the  definition  of  gust  velocities  without  spar.wise 
variation.  The  combinations  of  airplane  size  and  current  criteria  have  proven 
adequate,  however,  a  transport  with  a  significantly  larger  wing  span  than  the 
Boeing  747  would  likely  receive  special  attention  in  this  area. 

Similarly,  transport  airplanes  have  gradually  introduced  active  controls  and 
non-linear  systems.  Although  recognized  by  the  AC  (Reference  2)  and  special 
conditions,  the  basic  criteria  have  not  yet  been  changed. 

DISCUSSION 

For  each  class  of  flight  vehicle,  the  following  are  discussed  (if  applicable): 

°  Sunmary  of  current  regulations,  and  their  historical  development. 

°  Adequacy  of  current  regulations  for  conventional  configurations. 

°  Adequacy  of  current  regulations  and  methods  of  analysis  for  anticipated 
configurations. 

°  Identified  Problems. 

0  Actions  to  address  identified  problems. 


i 
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SMALL  AIRPLANES 


Snail  airplane  gust  criteria  (Part  23  of  Reference  i),  with  the  exception  of 
the  evaluation  of  dynamic  response,  are  the  same  as  the  transport  criteria. 

Most  configurations  of  small  airplanes  certified  to  date  have  been  relatively 
stiff,  slow  and  conventional.  Special  conditions  have  been  written  and  rule 
changes  proposed  to  cover  the  certification  of  canard  configured  airplanes. 

The  evaluation  of  dynamic  response  to  turbulence  has  not  been  considered 
necessary.  The  current  rules  are  presented  in  FAR  Fart  23  §23. 333(c)  and 
§23.341  as  Allows: 

§23. 333(c)  Gust  envelope.  (I)  The  airplane  Is  assumed  to  be  subjected  to 
symmetrical  vertical  gusts  in  level  flight.  The  resulting  limit  load  factors 
must  correspond  to  the  conditions  determined  as  follows: 

(1)  Positive  (up)  and  negative  (down)  gusts  of  50  f.p.s.  at  Vc  must  be  con¬ 
sidered  at  altitudes  between  sea  level  and  20,000  feet.  The  gust  velocity  may 
be  reduced  linearly  from  bo  f.p.s.  at  20,000  feet  to  25  f.p.s.  at  50,000  feet. 

(1i)  Positive  and  negative  gusts  of  25  f.p.s.  at  Vc  must  be  considered  at 
altitudes  between  sea  level  and  20,000  feet.  The  gust  velocity  may  be  reduced 
linearly  from  25  f.p.s.  at  2u,000  feet  to  12.5  f.p.s.  at  50,000  feet. 

(2)  The  following  assunptions  must  be  made: 

(1)  The  shape  of  the  gust  is  -  Ude  (  2«s  ) 

U  =  —  (1  -  cos  — -  ) 

2  (  25C  ) 


Where  - 

s  *  Distance  penetrated  into  gust  (ft.); 

C  =  Mean  geometric  chord  of  wing  (ft.);  and 

Ude  =  Derived  gust  velocity  referred  to  in  subparagraph  (1)  of  this  section. 
§23.341  Gust  Load  Factors. 

In  the  absence  of  a  more  rational  analysis,  the  gust  load  factors  must  be 
computed  as  follows: 


KqUdeVa 
- - 

498 (W/S) 

Where  - 


Kg  =  0. 8(^/5. 3+*^  =  gust  alleviation  factor; 

^g  =  2(W/S)/pCag  =  airplane  mass  ratio; 

U(je  *  Derived  gust  velocities  referred  to  in  §23.333(c)(f.p.s.); 

Density  of  air  (slugs/cu.ft.);  at  altitude; 

W/S  =  Wing  loading  (p.s.f.); 

C  =  Mean  geometric  chord  (ft.); 
g  =  Acceleration  due  to  gravity  (ft/sec. 2) 

V  >  Airplane  equivalent  speed  (Knots);  and 
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e  -  Slope  of  the  airplane  normal  force  coefficient  curve  Cna  P«r  radian  If  the 
gust  loads  are  applied  to  the  Mings  and  horizontal  tall  surfaces 
simultaneously  by  a  rational  method.  The  Minn  lift  curve  slope  Cl  per 
radian  may  be  used  Mhen  the  gust  load  Is  applied  to  the  Mings  only  and  the 
horizontal  tall  gust  loads  are  treated  as  a  separate  condition. 

FIXED  WING  TRANSPORTS 

°  Summary  of  current  regulations,  and  their  historical  development. 

The  starting  point  for  this  discussion  Is  the  gust  load  formula  (FAR  Part  25 
$25. 341(b)(3),  Reference  1).  A  good  overall  technical  discussion  of  gust 
criteria  development  Is  given  by  Noback  In  his  report  NLR  TR  82134U,  Reference 

4.  The  current  gust  load  formula  (frequently  referred  to  as  the  "Pratt 
Formula")  Mas  developed  based  on  an  evaluation  of  V-G  records  from  9  civil 
transports  from  1933  to  1950  and  presented  In  Reference  3. 

As  can  be  deduced  from  the  table  (Figure  1)  below,  these  were  essentially 
straight  wing,  stiff  airplanes  flying  slowly  at  low  altitudes. 


Air 

plane 

Design 
gross 
wt..  W, 
lb 

Wing 
area,  S, 
sqft 

Wing 
span,  b , 
ft 

Mean 
geometric 
chord,  c. 
ft 

Aspect 
ratio,  A 

Design 

cruising 

speed, 

Vr,  mph 

Estimated 

operating 

altitude, 

ft 

A 

13,400 

836 

74 

11.3 

6.6 

180 

5,000 

B 

18.560 

939 

85 

11.0 

7.7 

215 

5,000 

C 

41,000 

1.340 

118.2 

11.3 

10.4 

181 

5,000 

D 

50,000 

2,145 

130 

16.5 

7.9 

168 

5,000 

E 

25,200 

987 

95 

10.4 

9.1 

211 

5.000 

F 

45,000 

1,486 

107,3 

13.9 

7.8 

230 

5.0(H) 

G 

94,000 

1,650 

123 

14.7 

9.2 

271 

10,000 

H 

70,700 

1.461 

117.5 

13.6 

9.5 

224 

10,000 

J 

39.900 

864 

93.3 

10.1 

10.1 

256 

5.000 

Figure  I 


As  airplane  configurations,  speeds  and  cruise  altitudes  changed.  It  was 
recognized  that  airplane  dynamic  response  should  be  considered.  The  current 
rules  are  presented  In  $25. 305(c)  and  (d),  and  $25,341  as  follows: 

$25. 305(c)  where  structural  flexibility  Is  such  that  any  rate  of  load 
application  likely  to  occur  In  the  operating  conditions  might  produce 
transient  stresses  appreciably  higher  than  those  corresponding  to  static 
loads,  the  effects  of  this  rate  of  application  must  be  considered. 

$25. 305(d)  the  dynamic  response  of  the  airplane  to  vertical  and  lateral 
continuous  turbulence  must  be  taken  Into  account.  The  continuous  gust  design 
criteria  of  Appendix  G  of  this  part  must  be  used  to  establish  the  dynamic 
response  unless  more  rational  criteria  are  shown. 

The  airplane  Is  assumed  to  be  subjected  to  symmetrical  vertical  gusts  In  level 
flight.  The  resulting  limit  load  factors  must  correspond  to  the  conditions 
determined  as  follows: 
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(1)  Positive  (up)  and  negative  (down)  rough  air  gusts  of  66  fps  at  Vg 
must  be  considered  at  altitudes  between  sea  level  and  20,000  feet.  The  gust 
velocity  My  be  reduced  linearly  from  66  fps  at  20,000  feet  to  38  fps  at 
50,000  feet. 

(2)  Positive  and  negative  gusts  of  50  fps  at  Vc  must  be  considered  at 
altitudes  between  sea  level  and  20,000  feet.  The  gust  velocity  may  be  reduced 
linearly  from  50  fps  at  20,000  feet  to  25  fps  at  50,000  feet. 

(3)  Positive  and  negative  gusts  of  25  fps  at  Vq  must  be  considered  at 
altitudes  between  sea  level  and  20,000  feet.  The  gust  velocity  may  be  reduced 
linearly  from  25  fps  at  20,000  feet  to  12.5  fps  at  50,000  feet. 

The  following  assumptions  must  be  made: 

(1)  The  shape  of  the  gust  Is 

U<te  ZTs 

U  *  —  (1  -  cos)  — 

2  25C 


Where  - 

s  *  distance  penetrated  Into  gust  (ft); 

C  ■  mean  geometric  chord  of  wing  (ft);  and 

Ude  ■  derived  gust  velocity  referred  to  In  paragraph  (a)  (fps). 

(2)  Gust  load  factors  vary  linearly  between  the  specified  conditions  B' 
through  G',  as  shown  on  the  gust  envelope  In  525.333(c). 

(3)  In  the  absence  of  a  more  rational  analysis,  the  gust  load  factors 
must  be  computed  as  follows: 

KgUdeVa 

n  =  1  + - 

498  (W/S) 


Where  - 

0.88*g 

Kq  =  =  gust  alleviation  factor; 

5.3^g 

2 (/W/S) 

g  m - .  airplane  mass  ratio; 

Ude  3  derived  gust  velocities  referred  to  In  paragraph  (a)  (fps); 
^  »  density  of  air  (slugs  cu.  ft.);  at  altitude 
W/S  *  wing  loading  (psf); 

C  *  mean  geometric  chord  (ft); 
g  *  acceleration  due  to  gravity  (ft/sec^); 

V  ■  airplane  equivalent  speed  (knots);  and 
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slope  of  the  airplane  normal  force  coefficient  curve  Cna  per  radian.  If 
the  gust  loads  are  applied  to  the  Mings  and  horizontal  tall  surfaces 
simultaneously  by  a  rational  method.  The  Ming  lift  curve  slope  Cm.  per 
radian  may  be  used  Mhen  the  gust  load  Is  applied  to  the  Mings  only  and  the 
horizontal  tall  gust  loads  are  treated  as  a  separate  condition. 


The  first  step  In  the  development  of  the  necessary  techniques  to  Include 
airplane  dynamic  response  Mas  taken  toward  the  end  of  the  era  of  piston  engine 
transports.  At  that  time  It  became  apparent  that  the  earlier  airplanes  had 
satisfactory  service  and  safety  records,  even  though  no  provision  had  been 
made  in  their  design  loads  for  dynamic  effects  that  Mere  known  to  be  present. 
Thus  It  became  evident  that  the  d?-.ign  gust  velocities  had  been  set  high 
enough  so  that  for  these  airplane  no  Increase  In  design  loads  for  dynamic 
effects  Mas  needed.  On  the  other  hand,  it  Mas  apparent  that,  Mith  the  noted 
trends,  the  relative  dynamic  effects  might  Mell  increase.  Sooner  or  later, 
design  to  static  loads  alone  coulu  lead  to  a  structure  of  inadequate  strength. 


Consequently,  to  prevent  any  deficiency  in  strength  that  might  otherwise  have 
resulted  from  this  trend,  the  CAA  at  that  time  adopted  a  policy  which  was 
summarized  as  follows: 


"During  the  AIA-CAA  dust  Loads  Meeting  In  Washington,  it  was  agreed  that 
If  a  manufacturer  showed  that  for  his  new  model  the  percentage  increase 
in  load,  due  to  transient  effects,  was  no  greater  that  that  of  his 
previous  models,  it  would  not  be  necessary  to  design  for  the  increased 
load;  however,  if  the  increase  was  greater  than  for  the  previous  models, 
thia  increase  should  be  designed  for." 

This  policy,  reflecting  what  may  be  called  the  concept  of  "limited  dynamic 
accountability",  was  applied,  for  example,  in  the  design  of  the  Lockheed  Model 
1649  Constellation  and  the  Electra.  As  was  the  practice  at  that  time,  primary 
emphasis  was  placed  on  a  comparison  of  dynamic  magnification  factors  of  wing 
bending  moment. 

The  major  objection  to  a  continuation  of  this  type  of  approach  was  that  detail 
engineering  data  on  the  various  satisfactory  existing  airplanes  were  available 
only  to  the  manufacturers  of  those  airplanes.  Consequently,  a  manufacturer 
whose  past  airplanes  may  not  have  been  gust-critical,  or  for  other  reasons  may 
have  had  more  than  the  required  strength,  had  to  design  a  new  aircraft  to  more 
severe  criteria  than  the  manufacturer  whose  past  aircraft  happened  to  have 
less  margin.  Further,  with  few  exceptions,  no  criteria  short  of  "full  dynamic 
accountability"  were  available  to  a  manufacturer  who  had  no  previous  aircraft 
in  operation  with  a  long,  satisfactory  service  life. 

FAA  decided  to  develop  new  gust  criteria  using  the  power  spectral  technique. 

The  results  of  study  contracts  let  to  Lockheed  and  Boeing  for  the  purpose  of 
helping  FAA  define  procedures  and  criteria  are  summarized  in  References  5  and 
6  (AOS- 53  and  54). 

The  criteria  that  FAA  incorporated  into  the  regulations  as  Appendix  G  to  FAR  25 
(Reference  i)  In  September  19B0  were  the  result  of  extensive  negotiations 
between  FAA  and  U.S.  Industry.  Tho  primary  difference  between  the  criteria 
prescribed  In  FAA-ADS-53  and  current  FAA  criteria  are  In  the  specified  design 
gust  velocities  and  their  variation  with  altitude.  Tuese  FAA  criteria  pro¬ 
vided  the  basis  for  all  current  continuous  turbulence  criteria,  regardless  of 
the  certifying  agency.  For  example,  the  European  civil  regulations  as  of 
January  1967  are  specified  In  JAR-25,  Reference  9.  The  continuous  gust  design 
criteria  given  here  are  identical  to  those  given  in  Appendix  G  of  FAR  26  with 
one  major  exception.  JAR-25  makes  no  reference  to  reducing  design  gust 
velocities  for  airplanes  that  have  an  extensive  satisfactory  service 
experience  with  design  gust  velocities  that  are  less  than  85  fps.  This  will 
be  discussed  later. 


1-6 


As  specified  In  Appendix  6  of  FAR  25,  power -spectral  gust  loads  criteria  are 
presented  In  two  basic  forms;  the  design  envelope  analysis  and  the  mission 
(flight  profile)  analysis.  Provision  Is  also  made  for  a  modification  of  the 
design  envelope  analysis  that  considers  the  service  life  of  existing 
airplanes. 

The  design  envelope  criterion  Is  similar  to  past  discrete  criteria  as  well  as  to 
current  limit  design  maneuver  loads  criteria.  Operational  usage  of  the  aircraft 
is  Ignored.  Instead,  the  aircraft  response  Is  evaluated  for  a  specified 
design  envelope  of  speed,  altitude,  gross  weight,  fuel  weight  and  center  of 
gravity,  c.g.,  position. 

Appendix  6  of  FAR  2b,  Item  (b)(3)(1),  provides  for  reduced  design  values  of 
u  Hgma.  Specifically,  "Where  the  Administrator  finds  that  a  design  is 
coo*>arable  to  a  similar  design  with  extensive  satisfactory  service  experience, 

It  will  be  acceptable  to  select  U  sigma  at  Vc  less  than  85  fps,  but  not  less 
than  75  fps,  with  linear  decrease  from  that  value  at  2U.000  feet  to  30  fps  at 
8u,OuO  feet."  To  apply  the  reduced  U  sigma  values  requires  that: 

(1)  Transfer  functions  of  the  new  design  are  similar  to  the  prior 
designs. 

(2)  Typical  missions  of  the  new  airplane  are  substantially  equivalent  to 
thet  of  the  similar  design. 

This  modification  to  the  design  envelope  criterion  came  about  from  an  AIA 
proposal  to  the  FAA  after  extensive  studies  of  mid-range  to  long-range 
transports,  suen  as  the  L-1011,  Dc-9  and  DC-10,  and  the  Boeing  727,  737,  and 
747,  that  showed  U  sigma  of  75  fps  at  Vc  was  permissible  under  FAR  Appendix  G 
for  this  type  of  transport.  The  higher  U  sigma  values,  specified  by  the 
unmodified  design  envelope  criterion,  are  tore  appropriate  for  the  lower 
cruise  altitude  more-severe  types  of  operation.  The  more  severe  types  of 
operation  are  represented  by  short  range  or  commuter  operations  where  cruise 
altitudes  of  2u,Quo  to  3u,(xaj  feet  are  typical.  The  mid-.-ange  to  long-range 
airplanes  normally  have  cruise  altitude  in  the  vicinity  of  35,000  feet. 

As  originally  developed  in  FAA-ADS-53,  the  mission  analysis  approach  was  a 
"stand  alone"  method.  It  was,  however,  suggested  that  the  most  appropriate 
criterion  would  be  a  combination  of  the  design  envelope  approach  and  the 
mission  analysis  approach.  Combining  these  two  approaches  now  constitutes  the 
Mission  Analysis  Criterion  specified  in  FAR  25,  JAR-25  and  the  various  U.S. 

MIL  SPECS. 

In  addition  to  a  mission  profile  analysis  the  Mission  Analysis  Criterion 
requires  that  a  design  envelope  analysis  be  performed  similar  to  the  design 
envelope  criterion,  but  with  reduced  U  sigma  values  to  provide  a  design 
envelope  floor.  The  U  sigma  value  at  Vr  Is  specified  as  60  fps  from  U  to 
30,000  feet  with  a  linear  reduction  to  25  fps  between  30,000  feet  and  80,000 
feet.  The  Vg  value  is  still  1.32  times  the  Vc  value  and  Vq  Is  still  0.5  times 
the  Vc  value. 

There  was  a  transition  period,  during  which  dynamic  response  to  both  single 
gust  encounters  and  continuous  turbulence  was  considered.  Depending  on  the 
manufacturer,  there  was  a  different  balance  of  reliance  on  the  two  methods. 
Dynamic  analysis  of  the  single  gust  encounter  w r,  accomplished  by  requiring  an 
evaluation  of  the  1  minus  cosine  gust  shape  and  25  chord  length  of 
$25. 341(b)(1),  including  the  rigid  body  and  flexible  airplane  responses. 

The  Boeing  747  was  the  last  airplane  certified  during  the  transition  period. 

All  subsequent  new  airplane  certifications  commencing  with  the  Lockheed  L 1011 
and  Douglas  DC-10  used  gust  criteria  similar  to  those  currently  published. 
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Figure  2 

“Adequacy  of  current  regulations  for  conventional  conf iyurations . 

Sufficient  service  has  been  buiic  up  on  airplanes  certified  to  the  Appendix  G 
PSD  gust  criteria  to  be  confident  chat  the  criteria  are  adequate  for  these 
configurations.  However,  these  airplanes  have  gradually  included  more 
non-linear  systems  and  active  controls.  In  general,  the  group  of  airplanes 
identified  as  "flexible  high  subsonic  swept  wings"  in  Figure  2  included  only  a 
yaw  damper  in  addition  to  a  limited  authority  autopilot.  However,  the  last 
group,  wr.'ch  currently  includes  only  the  Lockheed  L-IOU-50U  and  the  Airbus 
A32U,  have  systems  which  interact  with  vertical  gust  response,  and  have 
non-linear  characteristics.  Even  though  FAA  prepared  special  conditions  for 
the  certification  of  these  airplanes,  they  discussed  the  technical  concerns 
without  defining  an  accepfable  means  of  compliance. 

“Adequacy  of  current  regulations  and  methods  of  analysis  for  anticipated 
configurations. 

As  more  airplanes  are  certified  with  increasing  levels  of  system  interaction 
with  structure,  the  need  for  changes  in  the  regulations  increases.  For 
example.  If  the  control  surfaces  deflect  at  high  rates  to  large  deflections  to 
maintain  a  desired  airplane  attitude  in  turbulence,  the  way  in  which  these 
effects  are  accounted  for  in  the  analyses  should  be  defined  by  the  FAA,  and 
not  left  to  the  Imagination  of  the  manufacturer.  In  addition,  concerns 
increase  about  the  need  for  consideration  of,  for  example,  spanwlse  variation 
of  gust  velocity. 

“Identified  Problems 

The  following  are  considered  to  be  problems  with  the  current  regulations: 

-  Difficulties  using  current  PSD  criteria  to  produce  design  loads. 

-  Definition  of  realistic  turbulence  to  evaluate  active  controls  and  gust 
load  alleviation. 

-  Criteria  for  reduction  of  PSD  Gust  Intensity  below  85  ft/ser. 


1-8 


The  difficulty  using  the  current  PSD  criteria  to  produce  design  loads  Is  how 
to  fit  them  Into  the  routine  by  which  design  loads  are  obtained  and  stress 
analysis  Is  conducted.  Normal  stress  analysis  practice  utilizes  design 
conditions  each  of  which  Is  defined  over  the  whole  of  some  major  structural 
component  at  a  given  instant.  Power-spectral  methods,  however,  do  not  result 
In  this  sort  of  design  condition.  They  lead.  Instead,  to  Individual 
design-level  values  of  load  of  equal  probability  at  various  points  in  the 
structure,  or  of  various  components  of  load  such  as  wing  shear,  bending 
moment,  and  torsion,  with  the  phasing  undetermined.  For  example,  it  Is  not 
determined  whether  maximum  up  shear  combines  with  maximum  nose-up  or  maximum 
nose-down  torsion  or  with  some  Intermediate  value.  This  difficulty  can  be 
circumvented  to  some  extent  by  determining  design-level  values  of  Internal 
loads  or  stresses,  such  as  front  and  rear  beam  shear  flows.  As  discussed 
earlier,  the  methods  used  by  U.S.  Industry  are  reviewed  In  detail  In  a 
separate  paper. 

The  definition  of  realistic  turbulence  to  evaluate  active  controls  and  gust 
load  alleviation  becomes  more  important  as  the  use  of  active  controls  becomes 
more  prominent  In  modern  airplane  design. 

The  criteria  for  allowing  a  reduction  of  PSD  gust  Intensity  below  85  ft/sec 
are  qualitative.  It  Is  not  possible  to  evaluate  transfer  functions  and 
typical  missions  to  arrive  at  a  value  between  85  and  75  ft/sec.  Typically  the 
jet  transports  which  cruise  at  30,000  to  40,000  feet  can  substantiate  the  use 
of  75  ft/sec,  while  the  turboprop  shorter  range  transports  which  cruise  at 
20,000  to  30,000  feet  have  found  to  require  the  use  of  85  ft/sec. 

“Actions  to  address  identified  problems. 

Two  of  the  problem  areas  Identified  above  are  being  actively  addressed. 

Most  Airworthiness  Authorities  and  Transport  Manufacturers  outside  the  U.S., 
and  some  smaller  U.S.  Transport  Manufacturers  Indicated  an  Interest  in 
developing  a  time  domain  continuous  turbulence  gust  analysis  that  would 
produce  compete  sets  of  correlated  loads  on  an  airplane.  FAA  agreed  to 
sponsor  an  evaluation  of  methods  Including  the  Statistical  Discrete  Gust  (SDG) 
analysis  method  (Reference  7)  proposed  by  Mr.  J.  Glynn  Jones  of  RAF 
Farnborough,  England.  Originally  Mr.  Jones'  Statistical  Discrete  Gust  (SDG) 
method  was  seen  only  as  a  possible  alternate  means  of  compliance  which 
produced  complete  sets  of  time  correlated  loads.  It  Is  now  recognized  in 
addition  as  a  possible  tool  for  directly  analyzing  the  response  of  nonlinear 
systems  suci.  as  gust  load  alleviation  to  continuous  turbulence.  An 
International  Ad  Hoc  Committee  has  been  formed  under  the  chairmanship  of  Mr. 
Terence  J.  Barnes,  FAA's  National  Resource  Specialist  for  Flight 
Loads/Aeroelastlclty. 

The  committee  met  for  the  first  time  on  May  22,  1986  and  agreed  that  the  first 
step  should  be  to  validate  the  "overlap"  between  PSD  and  SDG.  As  discussed  by 
Jones  In  Reference  8,  since  both  methods  reflect  a  variation  of  gust  velocity 
with  gust  length  (in  the  PSD  method  by  the  von  Karmann  spectrum,  and  In  the 
SDG  method  by  assuming  chat  gust  velocity  Is  proportional  to  gust  length  to 
the  power  1/3)  there  is  a  fixed  mathematical  relationship  between  the  methods. 
This  Is,  however,  valid  only  for  linear  systems.  NASA  has  agreed  to  assist 
FA®  *n  the  validation  of  this  overlap  between  the  two  methods.  In  the 
Interim,  several  manufacturers  Including  British  Aerospace,  Canadalr  and 
de  Havllland  Canada  are  conducting  their  own  evaluations.  A  workshop  Is 
planned  for  October  12/13,  1987  In  London  at  which  users  will  discuss  the 
status  of  their  evaluations,  and  have  an  opportunity  to  discuss  any  problems 
with  Mr.  Jones.  No  recommendations  for  use  will  be  made  until  the  method  Is 
well  understood  and  It  has  been  evaluated  on  several  representative  airplanes. 
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The  JAR  Authorities  are  considering  allowing  a  gust  Intensity  reduction 
similar  to  that  allowed  by  FAR  25  Appendix  G.  However  they  have  suggested 
that  the  qualitative  evaluation  required  by  Appendix  G  to  allow  a  reduction  In 
gust  Intensity  should  be  replaced  by  a  quantitative  evaluation  for  ease  of 
compliance,  and  to  possibly  allow  the  use  of  values  between  75  and  85  ft/sec.  _ 
A  simplified  mission  analysis,  where  the  scope  Is  reduced  to  a  small  number  of 
missions,  and  where  the  loads  are  evaluated  at  a  few  key  locations,  has  been 
discussed.  FAA  will  review  any  acceptable  proposal  by  the  European 
Authorities  with  U.S.  Industry,  and  consider  It  as  a  possible  replacement  for 
the  current  rule. 

Helicopters 

The  gust  criteria  of  FAR  Part  29  are  very  basic,  since  helicopters  are 
relatively  Insensitive  to  gust  encounters  Loads  analysis  techniques  are 
therefore  unsophisticated,  and  based  on  data  developed  in  the  1940  time 
period. 

The  concerns  regarding  helicopters  in  turbulence  relate  more  to 
controllability  and  fatigue. 

The  current  rules  for  gust  loads  criteria  are  presented  in  §29.341  as  follows: 

Each  rotorcraft  must  be  designed  to  withstand,  at  each  critical  airspeed 
Including  hovering,  the  loads  resulting  from  vertical  and  horizontal  gusts  of 
30  feet  per  second. 

Tilt  Rotor 


At  first  sight,  it  would  appear  that  a  tilt  rotor  configuration  could  use  the 
existing  helicopter  criteria  in  the  hover  mode  and  the  transport  criteria  in 
the  cruise  mode,  with  the  onl.y  area  of  concern  being  transition.  However,  due 
to  the  size  and  flexibility  of  the  blades,  and  the  configuration,  there  are 
significant  blade/wing  aerodynamic/aeroelastic  interactions  in  all  modes. 

This  subject  was  discussed  extensively  by  the  Airframe  Technical  Issues  Panel 
at  the  Powered  Lift  Conference  in  Fort  Worth,  Texas  on  23/26  June  1987.  The 
concensus  of  the  experts  was  that  these  concerns  justified  a  complete 
re-evaluation  of  the  gust  criteria.  Some  form  of  time  history  gust  criteria 
will  likely  be  proposed. 

National  Aerospace  Plane  (NASP)/0r1ent  Express 

The  Orient  Express  concept  focuses  on  sustained  supersonic  cruise  at  Mach  5  or 
6,  and  hypersonic  cruise  vehicles  derived  from  the  aerospace  plane  may  achieve 
speeds  to  Mach  10  and  beyond  in  the  altitude  region  of  100,000  ft.  Initiated 
under  OARPA,  the  U.S.  Air  Force  will  head  the  aerospace  plane  development 
program,  with  NASA  having  major  technical  responsibility.  This  is  intended  as 
a  multiple  purpose  airplane,  -  hypersonic  cruise  and  single  stage  to  orbit  -, 
propelled  by  an  airbreathing  (Scramjet)  propulsion  system.  The  ascent  dynamic 
pressure  will  be  high.  New  design  concerns  arise  for  this  airplane  due  to 
major  thermal  effects  on  static  and  dynamic  aeroelastlclty,  and  dynamic 
loadings  on  the  huge  liquid  hydrogen  tank.  1986  saw  the  release  of  funds  to 
begin  HASP  technology  development.  Aircraft  ronceptual  design  studies  are 
being  conducted  by  Boeing,  General  Dynamics,  Lockheed,  McDonnell  Douglas  and 
Rockwell  International.  These  contractors  will  recommend  gust  criteria  as 
appropriate  for  their  Individual  design  concepts.  At  the  present  time 
proprietary  aspects  prevent  revealing  the  various  proposed  criteria. 

Although  FAA  Is  not  directly  involved  In  the  Initial  development  phase.  It 
already  has  begun  thinking  about  the  job  of  certificating  a  Mach  5  or  Mach  6 
transport  for  commercial  use. 
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Cdnfinuovc  QimI  Design  Crtwria 
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Dm  eoctinuoue  gust  design  criteria  hi  tbs 
appendb  nut  be  wed  bi  estiblstdag  the 
dynamic  rwponse  cf  the  sbplnat  te  vertical  and 
lateral  eont'nuow  turbukncsa  culm  *  man  ra¬ 
tional  critorit  I*  wed.  The  fbfiowfcg  gust  lead 
rsqciremenis  ijpb  te  tiWni  annlyito  end 
Ml  n  tnvtlopt  nNjidK 
(*)  The  InH  suet  loads  atiBring  the  eeo- 
tinuoue  turbulence  concept  must  be  determined 
in  accordance  wtth  the  psvririoas  a f  either 
paragraph  (b)  or  paragraphs  (e)  aad  (d)  ef  this 
•weadk. 

04  Dm* a  metoap*  The  badt  hade 

not  be  determined  te  eeecpdknee  with  the 
Wowing: 

(1)  AD  critical  altitudes,  weights,  and 
weight  distributions,  as  specified  in 
S  26.881(b),  end  sB  critics,  speeds  within  the 
ranger  indicetsd  in  parnfnith  AXS)  ed  this  ap¬ 
pendix  imist  be  eceaidoed 

(X)  Vetoes  o?  A  (ratio  o 1  root-mean-equare 
incremental  load  rootmean-aqvars  gust 
velocity)  meat  be  determined  by  dynastic 
analysis.  Tbs  power  spoctnl  denaht  of  at- 
moaphei.  turbulence  must  be  as  given  by  the 
aquathm- 

e>x,  t+j<isselap 

r(o)«—  n^uavaiya 

where: 

0 -power-spectra)  density  ( fL/sec.)*/ 
red. /ft 

•root-mear-square  guat  Telocity, 
”  't/sec. 

0- reduced  frequency,  radians  per  foot 

L-2,500  ft. 

(3)  The  limit  loads  must  be  obtained  by 
multiplying  the  A  whies  determined  by  the 


dynamic  analysis  by  the  fallowing  taboo  of 
thcgwttilet.lt)  Ur 

(0  At  speed  V  :  Us -16  fpe  tow  gust 
releeity  to  the  btorual  0  to  10,000  ft 
aMtwb  and  to  Knaariy  daereaaed  to  go  fpe 
tow  gust  tsiodty  it  00,000  ft  altitude. 
Where  the  Administrator  finds  that  a 
design  is  comparable  to  a  shnder  design 
with  as c« nerve  satisfactory  service  ex¬ 
perience,  M  wfll  be  aoeepteble  to  sSeet  Ur 
at  V,  leas  than  SS  fpe,  Ut  not  toes  than  75 
fpe,  with  linear  decrease  from  that  value  at 
SO, 000  teat  to  00  fpe  at  10,000  lest  Th? 
Mowing  Aactoee  arts  be  token  into  aecount 
whan  eenssing  eocnparabBty  to  e  similar 

a  s  — 

OM|K 

(1)  The  transfer  function  of  the  new  design 
should  exhftft  no  unoeosl  eharacterietica  as 
compared  to  the  similar  design  which  wSl 
significantly  affect  twpouaa  to  tarimlance; 
e.g..  coalescence  of  mortal  response  in  the  fre¬ 
quency  regime  which  nan  result  in  a  signifi¬ 
cant  increase  of  loads. 

(2)  The  typical  mission  of  the  new  airplane 
it  substantially  equivalent  to  tint  of  the 
similar  design. 

(8)  The  similar  design  should  demonstrate 
the  adequacy  of  the  Ue  selected. 

00  At  speed  V,:  Uj  is  equal  to  1.82  times 
the  values  obtained  under  paragraph 
(bysyo  of  this  appendix. 

(SO  At  epeed  V„:  Ue  is  equal  to  W  the 
values  obtained  under  paragraph  (bX8X0  of 
this  appendix. 

Ov)  At  speeds  between  V,  and  Vc  and 
between  Vc  and  V„:  Ue  is  aqua)  to  a  value 
obtained  by  Sneer  Interpolation. 

(4)  When  a  stability  augmentation  system 
is  included  in  the  analysis,  the  effect  of 
System  nonlinearities  on  loads  at  the  limit 
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toad  total  naat  ba  raafatiealr  or 

MMuaftii  Jar. 


(a)  Mmim  aaa%cio.  Unit  Mi  aant  be 
<nwhi<  to  Mfto  with  Iho  Moafag; 

^s-isttrsiritts 

prMfcefawWefalheltaddMrftatieaacdlhe 


■Man  aagnaoto  or  Unto.  for  aaafeafa,  and 
•*«■••  or  affective  vafaaa  af  the  pertinent 
pwinitw  iihiifcn>tli»tt 

0)  For  web  of  tin  nfeaka  eagmaata  da- 

feed  andu  paragraph  fc)(l)  of  tNe  appendix, 
trainee  af  1  end  H.  neat  he  determined  by 
mllMt  A  b  MM  aa  the  rate  «f  rcot- 
imMtrn  IhwMI  land  to  rootnean- 
ogaare  gnat  vetocily  aad  \  fa  the  radfae  of 
gyration  M  tba  toad  power  initial  dander 
fraction  aboet  aero  frequency.  The  power 
apeotral  deMty  af  tha  ataioepberfa  tor- 
balance  onaat  ba  gfcren  by  tba  apatta  aat 
forth  la  paragraph  (N0)  af  tba  appendix. 

A  Far  aacb  at  tha  load  aad  atraaa  goan- 
tMaa  aobetod,  tha  ftognenqf  af  aaaaadaaea 
annt  ba  datoradaad  aa  a  ftmeticn  of  load  lava) 
by  ataaaa  «d  tha  aqoattja— 


t-eeieetad  time  fatarvml. 
y  -  aat  valaa  of  tha  load  or  atraaa. 
Y__.wvalne  af  tha  load  or  atraaa  la  one* 
fivalflgfat. 

N(y)>  average  number  of  exceadanoee  of  tha 
tadkmtad  ralua  of  tha  load  or  atraaa  fa 
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The  knowledge  of  tho  actual  wind  and  turbulone*  situation  along  tho  flight  path  of  an  aircraft  Is  an  important 
factor  in  tho  aroa  of  meteorological  and  aoronautlcal  research. 

In  this  paper  different  flight  tost  programs  for  the  on-hoard  Implementation  of  off-llna  and  on-lino  wind  and 
turbulence  measuring  systems  ars  prassntsd.  The  theoretical  principle  of  tho  determination  of  all  throe 
components  of  tho  wind  vector  Is  stated.  A  summary  a,*  the  installed  sensors,  the  data  acquisition  systems  a  id 
computsr  equipment  Is  represented  and  the  aassntlal  affects  of  ssnser  errors  on  the  accuracy  of  wind 
dotormlnatlen  are  discussed. 
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sensor  parameter 
Inertial  Navigatlen  Systam 
altituda,  bright 
vortical  spood  of  tho  aircraft 
vortical  acceleration  of  tho  aircraft 

distance  batwson  flight  log  and  Inertial  navigation  system 
Pulse  Coda  Modulation 
static  prssaura 
dynamic  prsssuro 

angular  vslocity  about  tho  aircraft  y-axes 
angular  vslocity  about  tho  aircraft  z-axes 
time 

total  temperature 
North  component  of  true  airspeed 
North  component  of  Inertial  velocity 
North  component  of  tho  wind  vector 
true  airspsad  voctor 
East  componant  of  true  airspssd 
East  component  of  fnartial  vslocity 
inertial  velocity 

East  componant  of  tho  wind  voctor 
Wind  voctor 

vortical  componant  of  true  airsposd 

vortical  component  of  inertial  velocity 

vortical  wind  component 

aircraft  angle  of  attack 

flight  leg  angle  of  attack 

aircraft  angle  of  sideslip 

flight  leg  angle  of  sideslip 

fllgh  path  angle 

differential  operator 

sansor  error  or  wind  componant  error 

pitch  angle 

longituda 

standard  deviation 

sum 

latitude 

bank  angle 

true  heading 

wind  direction 


Wind  and  atmospheric  turbulence  ars  oeasntial  parts  of  ths  weather  process.  Oonssquontly  ths  measurement  of 
the  wind  vector  and  its  tubulent  fluctuations  is  ths  dsmain  of  metssrologist  for  several  decades.  From 
aeronautical  point  of  view  wind  and  turbulence  ars  disturbance  variables  affecting  the  aircraft  dynamic,  ths  lead 
of  aircraft  and  pilot,  ths  passenger  comfort  end  and  wind  shear  and  downdraft  may  restrict  flight  safety 
especially  during  teko-off  and  landing.  Hones,  this  msteereleglcal  phenomena  have  to  bo  subject  of  intensive 
research  In  ths  flsld  of  aeronautic*,  toe. 

At  ths  Institute  far  Guidance  and  Control  of  ths  Technical  University  ef  Braunschweig  different  flight  test 


2-2 


program  far  the  iM  determination  an  beard  af  aircraft  have  fcaan  Initiated  during  the  laaf  tan  years  Ona  af 
thorn  ueed  an  AIRBUS  A300  af  tha  DEUTSCHE  LUFTHANSA  callacting  raw  data  during  taka-aft  and  landing  far  tha 
offline  wind  calculation  (Fig.  I).  Tha  experiments  damanatrated,  that  tha  standard-equipment  of  an  alrllnar  gives 
oufflclant  result  a,  If  thara  la  parformad  an  extensive  and  accurata  data  praparatlan  and  fitting  (1),C21.  A  hlghar 
level  af  accuracy  aa  wall  aa  an  Incraato  In  frequency  ranga  can  ha  ahtalnad  ualng  a  apodal  equipped  raaaarch 
aircraft.  On  board  tha  DO  28  raaaarch  aircraft  (Fig.  2)  af  tha  Technical  Dnlverelty  af  Braunachwalg  an  on-llna 
wind  and  turbulanca  maaauring  ayatam  haa  baan  Implemented,  which  calculataa  all  thraa  components  af  tha  wind 
vaster  In  raal  tima  at  a  campling  ratn  of  23  Hi  [3] .  A  modified  ayatam  haa  baan  davaloppad  In  calaboratton  with 
tha  AERODATA  FtugmeBtechnlk  UmbH  Juat  now  far  tha  aacand  raaaarch  aircraft  of  tha  TU  Braunachwaig,  a  DO 
(28  (Fig.  3),  ualng  mora  powarful  computer  and  aanaer  ayatama.  Tha  DO  128  la  bara  alao  for  tha  flight  toot  of 
the  nawaat  development  of  AERODATA  Flugmaltachnlk  GmbH  In  tha  flald  of  alrborn  wind  end  turbulanca 
maaauring  ayatama,  tha  aocallad  METEOPOD  (Fig,  4).  Thla  ayatam  la  Inatallad  In  a  alandar  body  attachad  undor 
tha  wing  of  tho  raooarch  aircraft  inatead  of  tha  auilllary  fual  tank.  Ita  advantage  la  tha  high  flailblllty  of  tha  hola 
wind  datnrmlnatlan  and  procaaalng  ayatam. 

Tha  raaulta  of  tha  maaauring  projects  ara  uaad  to  gat  mor  .  Information  about  wind  ahaar  and  turbulanca 
phanomana,  and  to  perform  invoatlgationa  of  tha  aircraft  raaponaa.  Furthermore  it  auccaoda  tha  development  of 
flight  control  ayatama  aa  wall  aa  wind  and  turbulanca  engineering  modala  for  flight  almulatlon  and  hatard 
Inveatlgatlon. 


t.  Tha  arfaefafa  a  f  oa-b—rd  wfaaf  and  lartWawa  AdargfMljw 

On  board  of  an  aircraft  wind  can  not  bo  moaaurod  directly.  Only  by  taking  tha  vectordlffaranca  between  tha 
inertial  velocity  a.'d  tho  true  alrapaad  H  tha  windvector  may  be  computed  (Fig.  5).  Fallowing  the  Gorman 
Aviation  Standard  the  wind  compononta  aro  defined  In  tha  aarth-ftaad  coordinate  ayatam  aa  ahown  in  Fig.  6:  A 
wind  blowing  from  tha  aouth  rapraaanta  a  positive  north  component  Uwa  .  a  wind  from  tho  waat  yialda  a  poaitlva 
aaat  component  vw_  ,  and  a  downdraft  indlcataa  a  pcaittve  vertical  wind  component.  Tha  compononta  eipreaaed 
In  earth-fixed  coordinates  can  be  written  In  tha  following  form: 
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The  egu.  3-5  ara  representing  tha  complete  computation  algorithm  for  tho  tranaformatlon  of  tho  true  airapeed 
compononta  In  tho  earth-fixed  coordinate  ayatam  canaldaring  any  kind  of  aircraft  manoeuvre  (31. 

Baaidaa  of  tho  inertial  velocity  compononta  (u  Kg  .  vKfl  .  w  Kg  )  and  tha  magnitude  of  tho  airapaad  IVI  five  anglea 
have  to  ba  maaaurad  on-board  tha  aircraft:  Tne  flow  anglea  of  tha  air  relative  to  tho  aircraft  (a,  g)  and  tha 
attitude  data  true  heading  T,  pitch  angle  G  and  bank  angle  *.  Since  tha  air  velocity  aanaora  may  not  ba  located 
near  to  tha  aanaora  for  tha  Inertial  data,  a  corrective  term  Q.  a  ^  haa  to  ba  added,  whore  Q-  la  tha  vector  of  tha 
angular  velocity  of  tha  aircraft  and  ^  tho  diatanca  from  the  inertial  aanaora  to  tha  air  data  aanaara.  For 
approximate  calculation  thia  affect  can  ba  axpraaaad  by  additional  torma  for  tho  angle  of  attack  and  angle  of 
aide  alip: 


AS 

(8) 

A* 

(7) 

with 


where  ^  I*  (in  ieogitudtnal  blatant*  between  the  Inertial  date  aaraara  and  «ha  air  data  aanaora.  The  esprostiens 
q  aa4  r  ara  representing  tha  pitch  rata  respectively  yaw  rata  at  tha  aircraft,  which  alas  ha va  t«  ba  measured 
en-beard. 
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t.l  St  MMflB  *  *•  AmmM  A90 0 

Tha  Aimus  wind  ahaar  prajact  atartad  In  tttO  C 23 .  Subject  af  thla  raaaarch  program  waa  tha  Investigation  af 
character,  mapnltuda  and  frequency  af  wind  varlatlani  In  tha  terminal  area  of  alrparta.  One  requirement  fer 
describing  statistical  prapartlas  af  wind  shear  is  te  gather  data  from  a  large  number  af  take-off*  and  landings. 
Qua  ta  tha  medium  and  shert  range  sarvica  af  LUFTHANSA  AlflSUS  A 300  In  tha  european  and  tha  north  afrlcan 
region  tha  high  number  af  operations  wars  guaranteed.  Tha  A300  Is  equipped  with  a  medarn  Aircraft  Integrated 
Data  System  for  recording  numerous  signals  and  sensor  outputs  about  the  flight  conditions  of  the  aircraft.  Most 
of  tha  required  parameters  shewn  In  aqu.  3-5  ware  available  en-beard  In  ena  way  or  another.  During  normal 
sarvica  air  datu  from  tha  Air  Oata  Computer  (ADC)  along  with  other  signals  la  sent  to  tha  Flight  Data  Acquisition 
Unit  (FDAU)  which  performs  tha  signal  processing  for  the  Performance  Maintenance  Recorder  (PMR),  see  Fig.  7. 
In  order  ta  overcome  legal  prablama  In  the  field  of  callacting  personal  data  and  to  consider  flight  safety  aspects, 
a  second  FDAU  and  PMR  wars  Installed  In  parallsl  to  tha  primary  sat  used  by  tha  LUFTHANSA  for  their 
maintenance  purposes.  Hsnco,  an  output  saquanca  of  24  parameters  with  a  sampling  ruts  up  tr  4  Hz  could  bo 
recorded  independently  af  airline  requirements. 

One  problem  was  to  make  the  INS-data  available  fer  acquisition.  On-board  no  INS-data  was  processed  by  FDAU 
for  recording.  Thus  a  apodal  airworthy  INS-Intsrface  had  ta  ba  built  and  Installed  which  concerted  and  reduced 
tha  32  bit  serial  INS-data  bus  (ARINC  561)  inta  parallel  12  bit  format  far  the  digital  input  channels  of  the  FDAU. 
This  Interface  has  alaa  was  used  to  hauee  a  timar  switching  the  PMR  an  and  off  and  spodfla  take-off,  landing  and 
ga-areund  conditions. 

The  recorded  raw-data  have  been  converted  by  LUFTHANSA  from  PCM  structure  Into  readable  data  at  the 
University  computer.  Here  an  off-line  data  procasslng  waa  parformad  according  to  equ.  3-5.  But  It  turned  out 
that  an  astanslve  post  processing  of  tha  raw-data  had  te  bs  performed,  like  filtering,  wlldpoint  check  and  error 
modelling,  to  get  acceptable  accuracy  of  the  results.  The  problem  af  Influence  of  sensor  errors  on  the  accuracy 
of  wind  determination  will  be  pointed  or)  later  on. 


3.2  The  maaaur*  i  lasfassiaf  sf  Iks  DO  M/Dm  lit  rsiisrrt  ajrcraff 

It  is  easy  to  realize  that  a  high  accuracy  of  on-beard  wind  determlnotion  requires  a  apadfldy  equipped  aircraft. 
In  this  case  the  sensor  equipment  and  its  location  can  be  chosen  with  reyard  to  the  specific  requirements,  and 
callbratlen  may  be  performed  when  ever  it  is  required. 

In  tSBO  the  TU  Braunschweig  acquired  tha  DO  2S  aircraft  frem  the  BODENSEEWERK  Qerltetechnlk  (BUT)  where 
It  was  formerly  used  as  a  company  raaaarch  aircraft.  Since  than  a  CAROUSEL  IV  Inertial  Navigation  Systom, 
powerful  computer  hardware  and  aavaral  sensors  have  been  installed.  This  was  a  prerequlsita  for  the  realization 
of  the  on-line  wind  measuring  systems. 

Tha  DO  2S  computer  system  consists  of  two  digital  cemputere  (Fig.  S).  In  the  main  cemputer,  a  NORDEN  11/34 
(military  varslan  af  tha  POP  11/34),  tha  wind  vaster  Is  calculated  In  real  time  at  a  sampling  rate  af  23  Hz.  Far 
handling  Input  and  output  af  analog,  digital  and  synchro  signals  a  MUDAS  processor  designed  by  DORNIER  System 
is  used.  Thla  procaasor  also  transfers  those  signals  ta  the  main  cemputer.  Metearelegical  data  can  be  viewed 
during  flight  an  a  CRT  terminal  af  tha  main  computer  an  an  alphanumeric  cockpit  display.  Far  esample  the  average 
of  tho  wind  vector  and  cerrospandlng  aircraft  data  Ilka  peaitien,  altitude,  Heading  and  flight  time  can  ba  displayed 
and  may  alaa  ba  selected  as  output  on  tha  printer.  Data  sterage  is  parformad  by  a  PCM  recording  system,  which 
can  sample  32  channels  at  BS  Hz  on  ana  taps  recorder.  As  another  quick  look  feature  tha  flight  path  and  the 
horizontal  wind  vector  are  platted  in  real  time  on  an  XY-racerder.  Fig.  g  shews  an  asample  with  a  tan  second 
average  af  tha  wind  vretar. 

Tha  lecatlan  af  tha  cemputer  and  tha  PCM  hardware  as  wall  as  the  lecatlen  af  the  essential  sensors  Is  Illustrated 
In  Fig.  10. 

Mast  parameters,  which  ara  required  far  tha  wind  vector  determination  can  be  measured  directly  sscapt  true 
airspeed  and  the  vertical  spaed  af  tha  aircraft.  On-beard  tha  DO  2B  the  vertlel  speed  Is  synthetically  derived  by 
means  af  cemplamentary  filter  using  tha  high  fraquency  infarmatien  af  the'  vertical  acceleration  signal  and  low 
frequency  Information  af  tha  barometric  altimeter.  This  filter  obtains  a  t-s-accuracy  far  tha  vertical  spaed  sf  tha 
aircraft  of  about  0.05  m/s  (33.  Tha  true  airspeed  Is  derived  In  tha  erdinary  way  from  dynamic  pressure,  static 
prasaura  and  total  temperature. 

It  la  planed  a  medium-term  replacement  ef  the  00  2B  by  tha  DO  12B  research  aircraft  which  has  bean  acquired 
by  tha  TU  Braunschweig  In  1005.  The  00  126  has  4n  air  fralm  similar  ta  tha  DO  2B.  But  instead  ef  two  piston 
anginas  it  la  equipped  with  turboprop  anginas,  which  give  improvements  in  flight  performance  and  nelss 
abatement. 

Tha  alrdata  sanssrs  ara  located  at  tha  tip  of  a  new  designed  carbon  fibre  composite  nose  beam  about  3  m  In 
front  sf  tha  aircraft  nosa.  Available  ara  tha  DORNIER  flight  lag  or  s  five-hole  prsbe  whle>  -an  be  completed  by 
flsed  wind  vanes  for  measurement  af  tha  higher  frequency  turbulent  flue  tuitions  (Fig.  11).  The  etgenfrequency  ef 
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the  !!•*•  beam  including  tht  tenser-matt  It  about  IS  Mx.  Tht  cutting-sff  frequency  af  tha  flitX  wins  ,ant«  It 
roughly  100  Hx.  Hance,  III*  lawtat  turbulence  wave  length  which  can  be  measured  by  thlt  system  Ktt  a  valut  af  2 
matrat. 

Tha  Inartlal  data  ara  measured  by  a  ttrap  dawn  Inartlal  navigation  system,  a  rleneywell  LASER-NAV.  Tha  concept 
far  l ha  data  acquisition  and  precasting  tyatam,  which  It  Implemented  at  present,  la  llluatratad  In  Fig.  12.  Tht..  idat 
It  ta  raallza  an  Intafratad  navlfatlan  tvatam.  Far  tha  calculatlan  af  hlfh  accuracy  petition  data  a  camblnatian  af 
INS-data  and  data  from  a  Alebal  Position  System  (warklnf  In  a  differential  mada)  It  projected .  Invaatlgatlan 
thuwt  that  thlt  complementary  tyatam  alta  praducaa  a  flight  path  velocity  calculatlan  af  vary  high  accuracy, 
which  rpducat  tha  arrar  In  wind  determination.  Tha  wind  datarmtnatlan  Including  tha  camplata  caardlnata 
trantfarntatlan  it  parfarmad  In  tha  main  camputar,  probable  a  ruggadltud  Micro-Vat.  Far  data  racerdtng  a 
streamer  tapa  racardar  It  utad. 


£J  Tkm  JMOOATA  AdETEOPOC 

Tha  AEflOOATA  FlugmaBtachnlk  GmbH  haa  devaleppad  ar  aircraft  and  hallcaptar  pod  far  application  In  air 
pollution  and  mataarolaglcal  raaaarch.  It  incorparataa  an  en-boarH  maaaurlng  tyatam  far  wind  turbulanca  and 
athar  mataarolaglcal  paramatara.  Thlt  tyatam  Integrate:,  tha  whale  venter  equipment  (Fig.  13)  Ilka  S-hala  probe, 
temperature-,  pratauro-  and  humidity-prebei,  INERTIAL  NAVIGATION  SYSTEM,  radar  altimeter,  at  well  at  tha 
real  time  data  protecting  tyatam  In  a  alandar  body,  which  may  be  attached  at  a  pylen  under  an  aircraft  wing  ar  at 
a  wlra  beneath  a  hallcaptar.  Lika  tha  aircraft  Integrated  ayitem  af  tha  DO  2S/I2S  tha  METEOPOD  perfermt  real 
time  precettlng  of  tha  wind  vecter  Indutivaly  campantatien  af  tha  aircraft  motion  and  camplata  pracaatlng  af 
navigation  and  aircraft  data.  The  apedflc  advantage!  af  METEOPOD  ara  turbulanca  meaaurement  with  a  tampllng 
ra.a  of  100  Hx,  axtramaly  email  dlatancee  between  tha  different  tantars  and  exchangeability  between  aircraft 
and  hallcaptar.  There  ara  no  maaaurlng  reaultt  available  up  ta  new,  bacauto  tha  In-flight  teat  program  with  the 
DO  t2S  raaaarch  aircraft  etarta  In  actabra  IMf. 


*■  Tba  efYlaa f  of  erne re  fa  meeaeremeef  aw  (ba  accuracy  af  tkm  wind  dhrtamM— Men 

An  Important  factor  In  attimatlng  arrera  In  maaeuremant  and  in  cerracting  aircraft  wind  determination  ara 
In-flight  callbratlonr  Completed  by  apecl.'ic  error  medeia  far  tha  3-dlmenatenal  wind  determination  a 
far-reaching  elimination  of  tanaor  orrora  can  ka  obtained.  Tha  arrar  affacta  arc  demonatratod  by  a  linear  error 
modal,  ualng  tha  north  component  of  tha  wind  vector.  Tha  calculated  value  u  Wg  canaitti  of  tha  true  component  u 
and  lit  error  Au: 


u  *  u  ♦  Au 
w,.c  we  W| 


For  tha  total  error  Auw,  tha  following  linear  modal  it  utad: 
5u„ 
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Tha  total  error  of  tha  wind  component  la  tha  turn  of  all  partial  darlvativea  multiplied  with  tha  correapondlng 
tenter  error  Af  (  .  The  example  In  Fig.  14  illuitrataa  the  arrera  of  tht  horizontal  wind  copmonantt  cautod  by 
errort  In  true  alrtpeed  and  angular  parameter!.  Theta  are  a  functlen  af  the  SINE  and  the  COSINE  af  the  true 
heading  while  of  couraa  tha  error  of  tha  vertical  wind  component  doaa  not  depend  on  truo  heading  C3J.  At  a 
aecend  general  atatement  can  be  mada  that  tha  wind  component  arrera,  which  ara  cautod  by  angular  paramatara, 
ara  directly  proportional  ta  tha  true  alrapeed  af  the  a.  aft.  Thit  alta  matna  an  Incraaaa  in  error  magnitude  of 
wind  speed  and  direction  with  Increaolng  trua  alrtpeed.  rig.  12  damenatratai  thlt  relation  far  the  vertical  wind 
component.  It  can  bo  aeon  fram  thit  flgura,  that  in  tha  cata  of  level  flight  (y  ,  •  •  0>  only  arrera  in  tha  angina  of 
attack  and  pitch  anO  an  arrar  In  tha  vertical  apaad  of  tha  aircraft  contribute  te  tha  total  error  in  tha  vertical  wind 
component. 

Bocauae  of  the  error  behaviour  of  tha  herizontal  wind  companantt  thown  In  Fig,  14,  t  reconstruction  of  senior 
errors  or  at  a  combination  of  tanaor  errors  can  ba  carried  out  by  flying  specific  flight  pattern.  Fig.  IS  Illustrates 
the  calculation  af  the  horizontal  wind  components  respectively  the  magnitude  ana  wind  direction  In  a  standard 
turn  C4I.  Curve  1  thaws  tha  typical  error  behaviour.  In  curve  2  tha  error  affacta  are  eliminated  using  tha  above 
mentioned  linear  error  modal.  Furthermore  thit  flgura  demonstrates,  that  tha  static  tenter  errors  sasentlally 
affect  the  mean  wind  accuracy.  Whareat  tha  accuracy  of  tha  maaaurad  turbulent  fluctuations  ara  influenced  by 
tha  tanaor  dynamic. 


#.  Etwofot  af  WeAf  tool  reoo/to  end  dwbtd  wjtd  mgggjt 

A  lot  af  wind  and  turbulanca  maaauramant  have  boon  carried  out  by  meant  af  tha  AIRBUS  A300  and  tha  DO  IS 
research  aircraft  In  the  field  af  wind  shear  invaatlgatlens,  aircraft  retpenta  on  wind  and  turbulonce  at  wall  at  In 
the  area  of  meteorological  experiments.  Flight  testa  for  the  In-fttght  measurement  of  pollution  trantpart  In  tha 
atmosphere  alee  have  bean  carried  out  In  tha  middle  af  thlt  year.  In  this  cate  tht  simultaneous  mrasurement  af 
the  wind  vector  It  an  important  factor  far  tha  Interpretation  af  the  data.  Samples  of  reaultt  of  the  different 
research  actlvitlae  will  ba  demonstrated  In  this  chapter. 

In  Fig.  IF  reaultt  af  the  AIRBUS  experiment  ara  praaantad  showing  wind  and  temperature  prafilat  in  connection 
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with  temperature  tnvcraiena.  On  the  nnn  bant  than*  ezamploa  Indicate  n  algnlf.tent  Inftuaaca  ef  laveralen  *n  wlnt 
epaod  ant  hrMtlw.  On  thn  other  hand  th«  towvpocotuco  profile  It  effecting  tht  tateaaity  af  turbulence.  Thtrt  it  t 
much  lower  wy»H  leteodty  rtm  tht  laccnloa  Itytr  IMiftrtt  with  tht  region  beneath  tht  Inverelen. 
Obvteudy  t  pampte  rate  ef  4  Hz  In  auffleieat  ta  Identify  thaaa  effecta. 

Tht  Mat  tart  ezarvptea  tea  tart  rtttarch  aatMtlea  ht  tht  araa  tf  law  Itatl  jot  phenomena  which  hart  httn 
tar  far  mat  la  tttptrtttea  with  tht  taatHata  ter  Meteor  tlagy  ant  Climatology  tf  the.  Untveralty  tf  Hannover .  Fig.  it 
ahtara  typical  profWoa  tf  wind  a  peat,  direction  nad  temperature  tf  a  lew-levei-let  moeeeuret  by  mtan  tf  tha 
LUFTHANSA  AMBUS.  Tbit  wind  phonemcaan  hat  hath  ft  ant  in  tha  narthara  part  tf  dermaay  approximately  in  101 
tf  alt  dghte  (81.  Typical  teataraa  arm 


-  a  Jat  Hfet  witf  prafita  It  aaanattlan  with  Intaaalra  greeted  baaed  temperature  Invaralan 

-  law  ttrhalaaaa  tntenelty 

-  hartaantai  hamoganlty  tf  tha  wM  Halt 


Additional  ataaaarlap  aetata  la  tha  dlipram  originate  from  matt  moaoaromaet  whith  hart  baaa  performed  70  ta 
100  km  ttataat  fram  tha  alrpart  Sraman.  TMa  under Haoo  tha  hariiattal  hamapanity  at  thli  phenomenon. 

Fig.  10  tamparaa  maaaaramaala  af  tha  DO  28  raaaarth  airtraft,  a  mataaralaglaal  mrat  ant  tata  af  a 
lew-level- Jot  angtnoorlng  matal  taratappat  at  tha  hnUtata  far  SaMaaaa  ant  C antral  (01.  Tha  aaampla  ahawa  a 
gaat  agreement  af  tha  tlfterant  maaaaramanta  ant  tha  matal.  ThU  aaglnaarlag  matal  haa  haan  uaat  for  tha 
approximation  af  nnmeroaa  maaaarat  lew-level-jeta  ter  haiai  t  iaraatlgatlaaa.  In  a  aim  liar  way  matailng  of  other 
haxartaaa  wlnt  phonomona  Ilka  ta wnhyrat  ant  warm-  ant  caM-franta  haa  haan  porfarmat  [11. 

In  ganaral  a  turbo  lent  wlnt  praftla  can  ha  aeparatat  in  a  largo  acala  moan  wlnt  or  front  ant  a  abort  acala 
turbulent  portion  (Fig.  20  ).  The  cumulation  of  tha  template  wlnt  flaw  flalt  tnay  ha  ayntheetzet  af  thaaa  two 
portiana.  Fig.  21  llhiatratao  tha  analyaaa  of  the  turbulent  fluatuatlana  ant  tha  fitting  hy  tha  Drytan  raapoctlvoly 
tho  *.  Ctarnmn  turbulence  matal.  Subject  af  thla  reaearch  la  the  evaluation  af  ralatlana  botwaon  tha  modal 
paramotora  ant  meteorological  continent. 

Tho  loot  oxamplo  rolatoa  to  DO  ••  maaaaramanta  af  tha  wlnt  altaatlan  influenced  by  a  met' Main  ridge  near  to 
Stuttgart  airport.  Invnatigatiena  by  mean  a  af  vary  aimpie  wind  matal  aancapta  have  temenatrated  that  flight 
aafaty  may  bo  effactot  taring  taku-ef  I  ant  go-araunt  by  loo-affaata  af  tha  hW  In  ona-onglno-out  operation  (71. 
Tho  moaaurot  wlnt  tata  Intlcala  a  much  higher  haxard  lovol  ant  faraa  ta  aarract  the  matal  concapta.  Fig.  22 
(lluatrataa  tho  raaulta  ef  a  taka-off  aimulatlan  with  an  angina  failure  at  *  I  ,  uelng  moaaurot  wind  data.  Tha 
aircraft  cnnrui  dear  tho  ebatadoa  on  thu  hill  aa  required  untar  thla  ceatltleaa. 


Thla  paper  proaanta  tho  principloa  of  in-flight  wlnt  ant  turbulanco  tatarmlnatlon  and  tho  hardwar*  realization  af 
on-board  maaaun.-g  ant  tata  procaaaing  ayatoma.  Tha  raaulta  of  aavoral  maaauring  project  a  Including  tower 
fly-by  and  compariaon  with  other  raforanco  data  a  how  a  ralatiao  high  accuracy  i.i  on-beard  wind  computation.  Tha 
horizontal  wind  apeod  la  calculated  with  a  pr-claien  botwoon  0.5  m/a  ant  0.7  m/a  (ta-ualuol.  For  tho  vortical 
wind  campanont  a  prodaion  of  0.3  m/a  can  bo  atatot.Sut  high  mcaaurlng  accuracy  procumoo  oxtonalvo  In-flight 
calibration  aa  well  aa  the  uaa  of  apadflc  error  motela  for  tha  campanaatlan  of  aonaor  orrora.  Further  Incraaaa  li 
prodaion  may  bo  obtained  by  modification  of  tho  maaauring  equipment,  eapeciallly  by  future  Inatallatlon  of  OPS  aa 
an  attitionr.l  aonaor  ayatom  tn  campanaata  tha  bad  long  term  accuracy  af  tho  INS-ayatem.  Firat  flight  taata  with 
different  OPS  racoivara  have  boon  performed  Juat  now  ant  a  raaoarch  project  in  thia  flalt  will  bo  atartod  In  1080. 
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Fig.  t :  Nwsurlng  tlrcnft  AIRBUS  A  300 


Fig. 4:  Th*  AEROOATA  METEOPOO 


Fig.  S:  Ddtorwlnition  of  tlio  Hind  Voctor 


Fi(.  6 1  Definition  of  Hind  Component* 
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Flj.13:  Herd**  re  locetion  in  the  HETEOPOD 
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«rror  of  the  TerUcol  wind  Component  due 
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Fij.lS:  Heesureaent  correction  by  aeons  of 
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Fig. 19:  Comparison  of  aircraft  and  mast  measurement  with  a  low  level  jet  model 
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MATCHING  P.S.D.  -  DESIGN  LOADS 

>>y 

R.  Nohack 

National  Aerospace  Laboratory  NLR 
P.0.  Box  153.  8300  AD  Rmmeloord, 
Tha  Natharlanda 


SUMMARY 

A  method  to  match  loads  obtalnad  with  tha  Design  Envelope  crltarlon  of  tha  P.S.D. -method  la  presented. 
Con a 1st ant  aats  of  design  load  conditions  can  ba  ganaratad  using  tha  corralatlon  coafficlanta  batvaan  tha 
loada. 

Two  of  thaaa  acts  ara  propoaad  for  practical  uaa. 


1 .  INTRODUCTION 


Tha  Power-Spectral-Density  (P.S.D.) -method  for  tha  calculation  of  airplane  loada  is  baaed  on  the 
assumption  that  atmospheric  turbulence  la  a  random  quasi-  stationary  Gaussian  process*  acting  as  Input  to  a 
linear  system*  tha  airplane.  The  output*  i.a.  loads*  accelerations*  stresses,  etc.  also  are  quasl- 
statlonary  Gaussian  processes. 

The  calculation  of  design  loads  y..  with  the  P. S.D. -method  can  ba  based  or  two  criteria  as  described 
in  references  1  and  2. 

The  first  one  is  the  Design  Envelope  criterion.  The  ratio  of  the  standard  deviations  of  the  load  and 
turbulence  is  multiplied  with  a  design  value  U  to  obtain  the  design  load  y^. 
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(i) 


it 

to 
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should  be  noted  that  here  and  in  the  following  the  design  load  or  stress  is  the  load 
turbulence.  For  the  real  design  value  the  lg-load  has  to  be  added. 

The  second  criterion  is  the  Mission  Analysis  criterion.  The  number  of  exceedances  of 
a  certain  flight-segment  (k)  is  calculated  with 


‘k 


(2) 


or  stress  due 
load  level 


Alk  is  the  ratio  Oj/o^  and  NQik  is  the  number  of  aero  crossings  of  load  y±  pertaining  to  segment  k. 

PLk*  P2k*  bik  aD<*  b2k  ,re  constants  describing  the  atmospheric  turbulence  in  segment  k. 
yglk  is  the  lg-load  of  load  y ^  in  segment  k. 

The  total  number  of  exceedances  of  load  level  y.  is  obtained  by  summing  Eq.  (2)  over  all  flight 
segments.  The  design  load  is  defined  as  the  load  for  which  the  number  of  exceedances  is  equal  to  the  design 
value  N.  . 

The  ^.S.D. -method  produces  the  design  loads*  however  not  the  mutual  relation  or  phasing.  If  the 
positive  and  negative  values  of  these  design  loads  are  used  to  calculate  a  stress,  depending  on  more  than 
one  load,  the  design  stress  usually  will  be  overestimated.  Besides  that  the  design  loads  will  not  be  in 
equilibrium. 

Hence  a  method  is  needed  to  combine  or  match  loads  such  that  stresses  as  calculated  with  such  a 
combination  of  loads  or  design  load  condition  will  give  a  good  estimate  of  the  correct  values  of  the 


In  the  discrete  gust  case  design  load  conditions  usually  will  be  defined  au  the  loads*  occurring  at 
the  same  time,  usually  at  the  time  that  one  of  the  loads  is  at  its  maximum  or  minimum  value.  Estimates 
for  the  stress  in  a  part  of  the  structure  can  then  be  calculated  for  each  one  of  the  design  load 
conditions. 

The  maximum  positive  or  negative  value  of  the  stresses  thus  obtained  will  be  the  dealgn  stress.  It 
will  be  clear  that  this  calculated  value  of  the  stress  wi]l  be  lower  than  or  just  equal  to  the  maximum 
value  of  the  stress.  The  stress  generally  will  have  its  maximum  positive  or  negative  value  not  exactly  at 
the  time  that  one  of  the  loads  reaches  ito  design  value  (see  Fig.  1). 

A  method  to  generate  dealgn  load  conditions  or  to  match  design  loads,  obtained  with  the  Design 
Envelope  criterion  of  the  P.S.D. -method  will  be  described  in  this  paper.  Only  the  case  with  two  loads 
will  be  treated.  The  derivation  for  the  general  case  with  N  loads  is  given  in  reference  3,  only  the  results 
will  be  given  here. 


2.  CORRELATION  AND  EQUAL  PROBABILITY 


In  the  P.S.D.^ethod  it  is  assumed  that  atmospheric  turbulence  is  a  quasl-atatlonary  random  process, 
with  Gaussian  probability-density  function  (p.d.f.)  and  with  normalised  power  spectrum  and  standard 
deviation  o  .  This  random  process  acts  as  input  to  a  linear  system,  the  aircraft.  The  outputs,  loads  and 
stresses,  have  the  same  properties  as  the  input. 


The  power  spectrum  of  the  output  y±  is 


v  -Mr 


(3) 
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The  ratio  of  tha  standard  deviations  of  output  y^  and  input  w  is 


Ih1w(«)I2 


(4) 


Hiw(«)  is  tha  transfer  function  of  output  y ^ 

According  to  tha  raqulraaants  (Refs*  1  and  2)  tha  design  load  has  to  be  calculated  with 

■Ii“.-3ir  <5> 


Th*  p.d.f.  of  tha  load  y^  la 
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and  it  follows  that  tha  probability  that  load  y  is  larger  then  tha  design  load  la 
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Tha  outputa  and  Input:  w  ar*  corralatad.  The  correlation  coefficient  between  outputs  and  la  (see 
Ref.  3) 
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J  Uiv(»H*v(jw)  +  H*w(jw)HJv{j.i1)}*w(w)dU 


(8) 


The  Joint  p.d.f.  of  two  loads  y(  and  y 2  with  correlation  coefficient  p;2  le 
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has  the  same  probability  density,  namely 
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Any  coabinatlon  of  loads  y^  «*nd  y2  such  thst 
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It  should  be  noted  that  Eq.  (10)  represent  an  alllpae. 


3.  DESIGN  STRESS  AS  FUNCTION  OF  TWO  LOADS 


(9) 


(10) 


(11) 


It  will  now  be  aaauaed  that  scresa  q  la  a  linear  function  of  loads  and  y^. 


(12) 


Tha  coefficients  Sj  and  depend  on  tha  dlaanalona  of  the  structure. 
The  design  value  for  the  etress  q  can  ba  calculated  with 


"d  '  *q  U0 

in  which  Aq  can  bs  expressed  as 


(13) 
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From  Eq.  (12)  follows 

yj«)  -  «!  +  «2  H^O#) 


(15) 


This  glvti  with  Eq..  (8)  nd  114) 

2 

Xq  •  ^  '  /  {tl  Hlw(J“)HIw(J“>  +  *1  *2  + 


+  •l  »2  Hlw(Ju)H2l((3u)  +  11*  b2w(5“)b2»<-I“)J*»<“)<1*‘ 
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It  should  b«  Dotsd  that  stress  q  is  a  linear  function  of  the  Geuselan  processes  y. .  This  implies  that 
q  also  la  a  Gt.ua  si  an  process  with  p.d.f. 
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The  probability  that  q  Is  larger  then  q.  -  A  U  Is  equal  to  the  probability  that  load  y  la  large*- 
yid  ■  \  Ua  <*<•  7> 

P(q  >  qd)  •  J  |l  -  «r£  (  ^ —)j  (18) 
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Ji(c)»  y2(t)  end  q(t)  '•  a  y.(r)  +  a2  y2(t)  ere  outputs  of  a  system  having  as  Input  the 
>se  v(t).  It  will  be  clear  thAc  if  q  ie  equal  to  Its  design  value  q  ,  that  then  y.  end  y9 

i - 2  _r ..  a  1  *- 


The  proce 
Geuesien  procesi 
can  have  all  values  that  satisfy 
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It  can  be  shown  (see  Ref.  3)  that  the  p.d.f.  of  y.  under  the  condition  that  q  -  q  ,  Is  a  Gaussian  p.d.f. 
with  mean  d 
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and  standard  deviation 
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The  design  lose  condition,  having  the  highest  probability  under  the  condition  that  q  is  equal  to  q.,  Is 
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The  correlation  coaftlclants  plt^  and  are  (Ref.  3) 


*1  °1  +  p12  *2  °2 


Plq 


P2q 


p12  *1  C1  +  *2  a2 
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q  q 

The  deelgn  load  condition  aj  glvan  In  Eq.  (22)  dafincs  tha  "optimal'Meaign  load  condition  for  atraaa  q 
with  coefficients  and  a2> 

Tne  locua  of  the  loads  y.  and  y~  will  now  ba  determined.  The  loads  y.  and  y2  will  ba  expressed  in 
tha  non-dimen  sJcnal  form 
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Tha  oon-dlrntnalonal  dealgn  loads  are 
yl 


X.  ■  - —  -  p . 

1  7ld  lq 


y2 

y2d  ’ 


P12  °2  +  °2 


(25) 


In  which  n^  “  a^ 
and  (aaa  Eq.  16) 


„i  +  2  p12  n,  »2  +  n2 


(26) 
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Eq.  (19)  btcowa  with  loads  and  y 2 


(27) 


Solving  Eq.  (25)  for  n  and  n.  and  insartlng  tha  result  In  Eq.  (27)  gives  the  relation  between  x.  and  5L. 
The  locus  of  the  points  (Xj,  x2)  is  the  ellipse 

+  2  *2  +  *2  "  1  "  p12  (2®> 

This  equation  represents  the  mb*  ellipse  as  Eq.  (10). 

An  example  la  given  In  figure  2. 

Eq.  (27)  represents  the  tangent  to  the  ellipse  In  point  (x.»  3L /*  This  poln>'  represents  the  optimal 
design  load  condition  (Eq.  22). 

The  distance  from  the  centre  of  the  ellipse  to  this  tangent  Is 


D  (n1#  n2)  - 


(29) 


The  coefficients  a.  and  a2»  and  thus  n.  and  n2  are  not  known  In  the  design  stage.  Besides  that  for  various 
parts  of  the  structure  different  values  of  a^  and  *  are  valid* 

All  points  on  the  ellipse  have  the  same  probability  as  the  optimal  design  load  condition.  Suppose  that 
another  design  load  condition  is  chosen  for  the  calculation  of  stress  q,  for  example  with  parameters  k.  and 
k2  instead  of  n^  and  n2  (see  Fig.  2). 


kl  +  P12  k2* 


or  y. 


kl+p12  k2 


5 12  kl  +  k2 


or  y  - 


p12  kl  +  k2 


U 

O 

°2  o 


(30) 


with 


°k  ’  kl  +  2  P12  kl  k2  +  k2 


(31) 


The  (Cress  that  will  b.  obtained  with  thesa  values  of  Che  design  load  condition  Is 


0. 


•i  h  +  *2  h 

°l(kl+  p  12  V  *  n2(p12  kl,f  k2}  JJo 

°k 

o  i 

qe  0 

w 


(32) 


The  line  through  the  pr;,'  *.tr  (r..,  x2)  18  defined  with  the  parameter  values  k.  and  k?,  and  parallel  to 
the  tangent  as  given  *n  Eq.  (2/)  Is  (aee  Fig.  2)  z 


"l  +  n2 


(33) 


The  distance  from  the  centre  of  the  ellipse  to  this  line  Is 

0  <V  V  -  (3A) 

/nl  +  n2 

It  can  b.  shown  for  N  »  3  (Ref.  3)  and  It  Is  aaslly  visible  in  figure  2,  thst  the  estimate  o  la  slwaya 
saall.r  than  or  equal  to  o  ,  or  q,  <  q  .  *** 

From  tbs  foregoing  can  bs  concluded  that  combinations  of  loads  y  and  y  can  bs  defined  that  have  an 
equal  probability  density.  Tha  locus  of  these  design  load  conditions  is  an  ellipse,  that  has  aa  a  tangent 
the  line  r.preeantiog  tha  relation  between  tha  design  stress  qd  and  tha  loads  y^  and  y^.  Each  combination 
of  paraaacers  k,  and  k  dsfinas  a  point  on  tha  allipss  and  each  point  on  tha  ellipse  daflnas  an  sstlaate 
q(  <  q^*  One  point  on  the  allipss  (the  point  of  contact  with  tha  tangent)  gives  the  exact  value  of  q^. 

Tha  locua  is  tha  case  of  N  loada  in  an  N-diaenalonal  aacond  order  surface,  for  H  ■  3  in  ellipsoid. 


A.  DESIGN  LOAD  CONDITIONS 

Aa  shown  in  tha  previous  paragraph  any  point  on  the  ellipse  can  be  chosen  to  represent  an  equal 
probability  daalgn  load  condition.  Ar  in  tha  case  of  tha  diacrats  gust  design  load  conditions,  cars  Bust  ba 
taken  to  chooaa  meaningful  condltlona  (for  axaaple  in  tbs  dlacr.ta  gust  case,  those  conditions  for  which 
one  of  the  loads  attains  its  mexiaua  or  alnlaua  value) . 

Two  eats  of  d'sign  load  condltlona  will  ba  proposed  for  practical  use.  Tha  first  one,  tha  correlated 
design  load  condition  is  coaparabla  to  tha  dlscrata  gust  case,  where  aaxlaua  loads  are  combined  with  tha 
loads  occurring  at  the  saaa  tiaa.  Tha  daalgn  load  conditions  in  tha  dlscrata  gust  case  will  produce 
relatively  low  estiaatts  of  tha  daalgn  stress  if  tha  phase  dlffarancsa  bstwsan  tha  loads  ara  large.  Tha 
saaa  will  ba  true  for  the  correlated  design  load  conditions  if  tha  correlation  cosfficianta  ara  saall. 
Therefore  a  second  sat  of  design  load  coodltlons,  tha  "aigan-vactor"  loads  is  proposed. 
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1*  Correlated  daalgn  load  conditions 

lech  dealga  load  condition  conalsta  of  one  daalgn  load  plus  the  correlated  value*  of  t'»  other  loads*  Thla 
la  analogous  to  the  dlacrete  gust  case*  where  each  design  load  condition  is  composed  of  one  design  load 
plus  the  values  of  the  other  loeda  at  the  ease  tine  that  the  first  load  reaches  its  maximum  value  (see 
Fig.  1). 

The  two  dealgo  load  conditions  can  be  generated  with: 


condition  1 

*1  1 

*2  0 

condition  2 
0 

1 

(35) 

Th.  ra.ult  1. 

*11  " 

yu 

*21  *  P12  yld 
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p12  y2d 

y22  ”  y2d 
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The  correlated  design 
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conditions  for 

H  loads  ara 
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2  3 
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yld 

P12  yid  p13 

yid . 

P12 

y2d 

y2d  P23 

y2d  . 

y  - 

p13 

y3d 

P23  y3d 

y3d . 

(37) 

y4d  ' 

P1N 

yNd 

P2N  yNd  . 

2.  Eigen-vector  design  load  conditions 

The  second  set  of  design  load  conditions  will  be  defined  as  the  loads  that  are  represented  by  the  end 
points  of  the  wain  axes  of  the  ellipse. 

These  points  can  be  determined  using  eigen-values  and  el  gen-vectors  of  the  matrix  R  with  the 
correlation  coefficients.  This  will  be  shown  for  the  2-dlmenslonal  esse. 

The  eigen-vector  ic^  la  definad  with  the  set  of  equations 


*1  +P12  *2+  •• 

.  .  •  iKj 

(38.) 

p12  *  *2  *  * ' 

..  *  1»2 

(38b) 

The  values  of  k  define  the  direction  of  the  main  axes. 

X  la  a  scale  factor. 

This  set  of  equatlona  hat  a  solution  only  If  the  determinant  la  equal  to  sero. 


I -A 
P12 


12 

1-i 


(39) 


utb 


Thi.  1.  an  N  ord.r  .quit ion  in  X,  giving  N  root,  or  .lg.n-v.lu*.  For  N  *  2  follow. 


*1  ‘  1  +  p12  *2  '  1  -  p 


12 


Ina.rtlng  th...  .olution.  in  Eq.  (38)  giv*. 
for  X, 


1  • 


'll 


for  1 


2  -  *21 


12 


22 


(40) 


(41) 


Th.  nor»»Xi..d  .lg.n-v.ctor.  will  now  b.  u«*d  to  d.fln.  th.  d.ilgn  lc.d 
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(42) 


*11*  p12  *12 


*21  *  p12  *22 


Th.  d..lgn  load  condition,  b.coa.  (...  Eq.  30) 
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It  cm  ba  proven  that  ofcl  •  A^  (Raf.  3). 
for  the  2  dimensional  casat 
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from  equations 

(38.  ud  b) 

follow. 
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These  resulte  inserted  in  I 
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gives 
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The  eigen-vector  deelgn  load  conditions  for  N  loads  are  with  k  as  the  normaliaad  eigen-vector: 
condition 
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The  design  load  conditionr  are  presented  in  figure  3.  The  point  opposite  to  the  ones  defined  above  produce 
the  same  design  load  conditions,  however  with  opposite  sign. 

The  deelgn  load  conditions  as  defined  have  special  properties.  When  the  design  is  finished  and  the 
dimensions  have  been  defined  it  is  then  possible  to  calculate  the  correct  value  of  the  stress  q  “  q ..  Of 
course  the  correct  value  can  be  calculated  also  with  Eq.  (16),  using  the  o.  end  p , .  values.  Assuming, 
however,  that  only  the  design  load  conditions  are  available  to  the  stress  office,  the  design-stresses  can  be 
calculated  with  the  following  rules. 


a.  Using  the  astlmatas  q  .  of  the  stress  as  calculated  with  the  correlated  design  load  conditions. 
The  square  of  the  design*itress  q.  is  equal  to  the  sub  of  the  products  of  the  i-th  estimate  and  the 
stress  due  to  the  i-th  design-load. 

S  ■ 1  yu  Si  <«> 


The  estimates  q#i  in  the  2-dimei.alonal  case  are 

51  ■  *1  yld  +  «2  p12  y2d 

52  ‘  *1  p12  yld  +  *2  y2ci 


(50) 


Inserting  this  in  Eq.  (A9)  and  using  Eqs.  (5),  (13)  and  (16)  proves  the  rule  for  the  2-dlmenslonal  case, 
b.  Using  the  estimates  of  the  stress  as  calculated  with  the  eigen-vector  design  load  conditions 

The  square  of  the  etress  la  aqual  to  tha  sum  of  ths  square*  of  the  estimates 


(51) 


The  iiriMCii  in  the  2 -dimensional  cant  art 


Inserting  thla  In  Kq.  (51)  and  using  Eqe.  (5).  (13)  and  (16)  provaa  tha  rula  for  Cha  2 -dimensional  caaa. 
Proof  a  for  these  rulaa  for  cha  N-dlmenaional  caaa  ara  given  In  rafaranca  3. 

It  alao  la  poaalbla  to  daflna  daalgn  load  conditions  that  provlda  a  lower  and  an  uppar  linlt  for 
atraaa  qrf.  Tha a*  condition*  ara  daacrlbad  In  Appaodlu  A. 


5.  DISCUSSION  AND  CONCLUSIONS 

It  haa  baan  ah  own  chat  it  ia  poaalbla  to  generate  equal  probability  daalgn  load  condltlona  ualng 
P.S.D. -daalgn  loada  obtalnad  with  tha  Daalgn  Envelope  criterion,  and  tha  corralatlon  coafflclanta  batvaan 
thaaa  loada.  Tha  corralatlon  coafflclanta  can  ba  calculatad  aaally  togathar  with  tha  X- values. 

Tha  Matrix  of  corralatlon  coafflclanta  can  ba  uaad  to  daflna  an  N-dlmanelonal  aurfaca.  For  tha  2- 
dimenslonal  caaa  thla  aurfaca  raducaa  to  an  ellipse*  Each  point  on  tha  aurfaca  daflnaa  an  aqual  probability 
daalgn  load  condition.  It  can  ba  shown  (Ref.  3)  that  aoch  a  condition  la  In  equilibrium. 

Each  daalgn  load  condition  can  ba  uaad  to  calculata  an  estimate  for  tha  atraaa  In  a  point  of  tha 
atructura.  Each  aatlmate  la  lovar  than  tha  corract  valua.  Only  ona  point  on  tha  N-dlaanalonal  aurfaca 
represents  a  daalgn  load  condition,  that  will  give  tha  corract  valua  of  tha  atraaa. 

Two  aata  of  N  daalgn  load  condltlona  ara  propoaad  for  practical  uaa.  Thay  Hava  baan  choaan  auch  that 
It  can  ba  axpactad  that  at  laaat  ona  of  tha  estimates  will  davlata  not  too  Much  from  tha  correct  valua. 
Thla  howavar  can  not  ba  guarantaad.  Fortunataly  tha  choaan  aata  of  daalgn  load  condltlona  both  hava  tha 
proparty  that  tha  corract  valua  of  tha  atraaa  can  ba  calculatad,  ualng  tha  estimates.  Thla  knowledge 
can  than  ba  uaad  to  radafina  tha  dinanaiona.  Nota  that  with  tha  dlacrata  guat  mat  hod  auch  a  check  la 
not  poaalbla. 

Tha  nathod  also  can  ba  uaad  to  ganarata  daalgn  load  condltlona  for  tha  coMplata  atructura,  for  exanple 
tha  wing.  A  probleu  arlaaa  If  atraaaaa  hava  to  ba  calculatad  for  a  auction  of  tha  atructura  for  which 
tha  daalgn  valuaa  of  the  loada  and  tha  corralatlon  coafflclanta  hav*  not  baan  calculatad.  It  than  aeaua 
a  logical  approach  to  interpolate  within  a  daalgn  load  condition.  Thla  however  laada  to  inconalatant 
valuaa  of  both  loada  and  atraaaaa  In  that  aectlon. 

Tha  mm  problew  arlaaa  In  a  dlacrata  guat  analyala  if  tha  daalgn  load  condltlona  are  defined  aa  tha 
loada  occurring  at  tha  mm  tine.  Thla  problan  ia  dlacuaaad  In  wore  detail  In  reference  3.  A  poaalbla 
aolutlon  la  to  lnterpolat  batvaan  tha  correapondlng  daalgn  loada  and  correlation  coafflclanta  of  tha 
adjoining  taction*. 

A  n unbar  of  hlnta,  that  nay  ba  uaaful  In  tha  application  of  tha  Method  la  given  in  rafaranca  4. 

Tha  methods  for  tha  determination  of  aqual  probability  daalgn  load  condltlona  can  not  ba  uaad  if  tha  loada 
are  obtained  with  the  Mission  Analysis.  Tha  correlation  coefficient*  in  that  caaa  ara  not  defined.  However 
It  ia  poaalbla  to  approximate  tha  aqual  probability  desljn  load  condltlona.  Tha  correlated  daalgn  load 
condltlona  In  that  caaa  conslat  of  tha  Median  value  of  tha  loada  under  tha  condition  that  ona  of  tha  loada 
exceeda  ita  daalgn  valua.  Tha  derivation  and  also  tha  application  ia  rather  involved.  A  description  la 
given  In  rafaranca  3. 

Tha  support  of  tha  Netherlands  Agency  for  Aerospace  Programs  (NIVR)  for  thla  Investigation  la 
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AfftMBIX  A  Upper  and  lower  Halt  design  load  condition? 

One  poaalble  aat  of  design  load  conditions  that  will  produce  conservative  estimates  of  the  atreaa, 
consists  of  the  combinations  of  positive  and  negative  values  of  the  design  loads.  The  estimate  of  the 
stress  calculated  with  these  design  load  conditions  (d.l.c.)  will  he  (nuch)  higher  than  the  correct  value 
oi  the  stress  g^. 

In  the  following  will  bo  shown  that  It  is  possible  to  define  d.l.c.  that  produce  at  least  one  value 
for  tbs  stress  larger  then  g.,  but  generally  much  lower  than  those  obtained  with  the  design  loads.  The 
nsalaua  possible  error  cen  be  eetsbllshed.  Using  this  result,  d.l.c.  will  be  defined  that  give  a  lower 
Halt  for  stress  q,. 

The  proposed  conservative  d.l.c.  are  based  on  the  eigen-vector  d.l.c.  In  the  cess  of  N  losds  also  N 
elgen-vsctor  d.l.c.  are  defined.  The  conservative  d.l.c.  In  this  case  consist  of  N  sets,  each  consisting  of 
2  to  the  power  H-l  d.l.c. 

The  a-th  set  la  equal  to  the  a-th  algen-vector  d.l.c.  plus  or  alnua  a  fraction  c  of  each  of  the 
other  algen-vector  d.l.c.  The  1-th  d.l.c.  of  this  set  esn  be  presented  In  vector  notation  with  y  as  the 
3-th  aigen-viictor  d.l.c..  as 
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Th«  A-th  d.l.c.  Is  defined  vith  one  of  the  poaalble  2  combination*  of  plue  and  alnua  eigne. 
The  atreaa  due  to  an  elgen-vector  d.l.c.  la 
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and  It  follows  that  the  atreaa  due  to  the  d.l.c.  of  Eq.  (Al)  Is 
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The  maximum  estimate  for  the  atreaa,  as  produced  by  the  a-th  set  is  equal  to 
a-1  N 
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If  c  <  1 

One  of  the  estlaatea  9  will  be  the  largest  one.  Without  loss  of  generality  It  can  be  assumed 
that  It  le  the  tatlaatS’of  tho  first  set,  hence 

*.l  *  «  *  2.  N  <A5> 

It  follows  that 
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(A6) 

It  will  now  be  shown  that  the  eetlaate  1),]  1*  *  conservative  estimate  for  certain  values  of  c. 

The  square  of  Is 
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Th«  sum  of  tha  aquaraa  of  tha  astlmataa  lq#J  of  th«  algm-vactor  d.l.c.  la  tqual  to  tha  aquara  of 
atraaa  and  it  follova  that 
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TIm  ttlatlva  error  lo 

•  >Jl  +  !j  -  l  -  t/l  *  T  -  1 

with  r  »  -5*  *  -V 
^  1  '.i 

mo  raluo  should  be  positive,  but  so  email  as  possible, 
eeeood  car*  lw  positive  (or  c  »  0.  bscauao  |a  , |  <  |a  , | . 

Ths  first  tan  la  ssro  or  positive  l( 

le  -  l  ♦  e1  h  0  or  c  8  /l  -  1 


(A10) 

(All) 


The  last  Cara  of  Eq.  (A9)  lo  pooltlvo.  Tho 


(A12) 


Thla  aloo  turoo  out  to  be  th«  loveet  allowable  value  of  c.  Aaeume  that  q  .  ■  0*  except  for  J  »  1  and  2. 
Tho  third  (pooltlvo)  com  of  Kq.  (19)  la  oquol  to  aoro  aod  thla  aquation  bomooa 

I  -  (2c  -  1  ♦  c1)  q^  ♦  2c  (Iq^l  Iq^l  -  ^,) 

-  -  (1  -  cJ)  q*2  +  2c  |,#ll  ,tl||  (A13) 

Tho  ratio  |q  _|/|q  ,|  rangaa  from  taro  to  ona  and  It  follows  that  E  and  thua  F  ara  pooltlvo  If  c  lo  aqual 

to  St  -  1. 

Tha  aaxlaua  valua  for  F  la  obtalnad  aa  follow. 

Assume  that  tho  darlvatlvoa  8F/8q  .  (J  ■  2,  N)  ara  aatobllahad.  Than  for  symmetry  raaaona  tha  extreme 
valua  If  it  exists,  will  occur  at  q>2  -  q#3  •  ...  q#JJ. 

Inter ting  thla  raault  In  Eq.  (A9)  and  (A14)  and  taklnt  tha  dorlvativa  IF/Iq  2 ,  taking  Into  account 
that  c  ■  /2-1  It  la  found  chat  * 

q,2  -  q.3  -  •••  t.h  -  c  «#l  <*u> 

Thla  raault  lnaartad  lo  Eq.  (A9)  and  (All)  glvaa 

r  -  (N-l)  c*  (A15) 

and  It  follows  that 

0  S  r  S  (N-l)  c2  (Alb) 


In  tha  table  tha  masher  of  estimates  and  tha  maximum  possible  relative  errors  ara  given  for  noma 
values  of  N 


N 

number  of 

estimates 

Miimm 

relative 

error 

/I+F-l 

correlated  + 

•lg.n-valuea 

conservative 

2 

4 

4 

0.0824 

3 

6 

12 

0.1389 

A 

8 

32 

0.2302 

5 

10 

ao 

0.2986 

6 

12 

192 

0.3630 

Tha  conservative  design  load  conditions  aa  defined  In  Eq.  (Al)  with  c  ■  *'2-1  produce  an  upper  limit  for 
etreea  qd-  It  la  now  also  poeelble  to  define  a  lower  Halt  for  atraas  qrf ■ 

from  lq.  (A10)  and  (Alt)  follow. 


9d  M.  *  ,d  Vl+F 


(A17) 


^  lo  tho  aixlra  eatlmste  obtalnad  with  tho  uppor  limit  d.l.c.. 
E$»  (A17)  con  be  written  a« 


(A18) 


and 

Tha 


thua 
9*  1 


~  ‘  qd  * 

/T+r  a  * 

lovor  limit  4  for  q^ 

y  .  - 1 -  0 

*l  /Hi 


can  ha  obtalnad  with  tha  lower  Unit  d.l.c. 


(Al») 

(sea  Eq.  Al) 

(A20) 


It  can  be  shown  that  tha  points  rapraaantlng  tha  lower  limit  d.l.c.  ara  located  on  tha  M-dlmenslonel  second 
order  surface  rapraaantlng  tha  locus  of  tha  aqual  probability  d.l.c..  The  upper  and  lower  limit  d.l.c.  for 
tha  case  M  •  2  ara  shown  In  figure  Al. 

Tha  correct  valua  for  qd  can  ha  expressed  aa  a  function  of  tha  estimates.  It  can  ba  shown  that 


(All) 


J-12 


1  ’!k 

k-i 

n ei— 
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with 


•■>4  F  -  e1  (H-1) 

The  upper  and  lower  limit  d.l.c.  hara  the  adviBtaia  that  a  range  tor  tha  correct  valua  of  atraaa  can 
ha  eetabllahed.  Tha  drawback  ia  that  a  much  larger  nuaber  of  eatlaatea  have  to  ha  calculated. 


Pig.  A1  Upper  and  lower  Halt  deal go  load  condition 


i 


SimMiri aatlWant-Uc.  - 


A  SUMMARY  OF  METHODS  FOR  ESTABLISHING  AIRFRAME  DESIGN 
LOADS  FROM  CONTINUOUS  GUST  DESIGN  CRITERIA 


4-1 


Richard  If.  Moon 
Design  Specialist  Senior 
Dynamic  Loads 

Lockheed  Aeronautical  Systems  Company 
P.O.  Box  551,  Burbank,  CA  9-1520  USA 


SUMMARY 

Continuous  gust  design  criteria  for  airframe  design  are  specified  in  FAR  25,  JAR-25  and 
various  United  States  military  specifications.  Two  forms  of  criterion,  the  design  envelope 
approach  and  the  mission  analysis,  are  usually  referenced  as  “an  acceptable  means  of  compliance." 
However,  these  criteria  do  not  provide  methods  of  applying  the  statistical  results  to  the  design  of 
the  structure.  Development  of  such  methods  is  left  to  the  imagination  of  the  airframe 
manufacturer,  subject  to  the  approval  of  the  certifying  agency.  Some  of  the  methods  that  arc 
currently  used  by  United  States  airframe  manufacturers  arc  summarized  here.  Continuous  gust 
design  requirements  from  various  certifying  agencies  are  reviewed.  A  brief  discussion  is  also 
provided  on  the  methods  employed  to  include  the  effect  of  the  L-1011  Tristar  active-controls 
wing  load  alleviation  system  on  the  loads  due  to  corrective  roll  control  in  turbulence. 


INTRODUCTION 

This  paper  is  composed  of  three  related  sections. 

First,  a  brief  review  is  provided  of  continuous  turbulence  requirements  specified  by  various 
certifying  agencies.  It  is  the  purpose  of  this  review  to  provide  background  material  for  the 
reader.  Design  gust  load  continuous  turbulence  psd  requirements  are  quite  consistent  in  their 
essential  aspects  among  all  certifying  agencies.  Discussion  of  some  of  the  underlying  concepts 
that  contributed  to  current  design  gust  criteria  is  includea.  Much  of  the  discussion  reflects  the 
work  of  Frederic  M.  Hoblit  as  presented  in  a  pre-publication  version  of  Reference  1.  Mr.  Hoblit 
was  the  Lockheed-California  Company  lead  engineer  on  the  FAA  contract  during  1964  -  1966 
that  resulted  in  repjort  FAA-ADS-53,  which  provides  the  basis  for  current  continuous  turbulence 
gust  loads  requirements. 

The  second  part  of  this  paper  deals  with  the  problem  of  applying  continuous  turbulence  design 
criteria  to  the  design  and  sizing  of  structure.  In  short,  what  do  you  do  with  psd  statistical 
parameters  after  you  have  them?  There  has  been  very  little  written  with  respect  to  practical 
application.  A  summary  of  some  of  the  methods  used  by  United  States  airframe  manufacturers 

is  presented.  These  methods  are  dependent  on  the  nature  of  the  airplane  to  be  analyzed 

(similarity  to  prior  designs  and  anticipated  operation),  the  criticality  of  the  structure  to  gust 
loading,  and  the  complexity  of  the  structure.  Different  methods  arc  often  applied  by  the  same 
manufacturer  in  the  analysis  of  different  airplane  configurations,  or  for  that  matter,  for  different 
components  of  the  same  airplane.  Because  a  mixture  of  procedures  is  likel;  be  used  by  all 

companies,  these  methods  are  not  associated  with  a  single  company,  but  are  presented  as  a 

composite. 


As  the  use  of  active  controls  becomes  more  prominent  in  modern  aircraft  design,  the  adequacy 
of  the  continuous  gust  design  criteria  and  the  associated  methods  of  assuring  adequate  strength 
become  of  concern.  The  last  section  of  this  paper  provides  a  summary  of  the  studies  that  were 
performed  to  include  the  effects  of  the  Lockheed  L-1011  Trisur  active  control  system  on  the 
determination  of  gust  loads.  The  active  control  system  developed  for  the  L-1011  Tristar  included 
a  wing  load  alleviation  system.  Although  wing  loads  due  to  gust  were  significantly  reduced, 
loads  due  to  corrective  roll  control  in  turbulence  were  unaffected  and,  as  a  result,  became 
significant  Considerable  modification  of  the  standard  psd  design  gust  procedures  was  required  to 
account  for  this  effect. 

REVIEW  OF  CRITERIA 

The  FAA  incorporated  explicit  continuous  turbulence  power-spectral  gust  loads  criteria  into 
Appendix  G  of  FAR  25,  Reference  2,  in  September  of  1980.  These  criteria  were  the  result  of 
nearly  20  years  of  study  involving  close  coordination  between  the  FAA  and  the  manufacturers. 

The  first  attempt  to  develop  a  comprehensive  set  of  power-spectral  gust  loads  criteria  that  had 
general  applications  was  the  study  conducted  by  Lockhccd-California  Company  in  1964  -  1966, 
under  contract  to  the  FAA.  This  study  produced  two  reports;  FAA-ADS-53,  “Development  of  a 
Power-Spectral  Gust  Design  Procedure  for  Civil  Aircraft,”  Reference  3,  and  the  companion  report 
from  Boeing,  FAA-ADS-54,  Reference  4.  The  criteria  formulated  in  FAA-ADS-53  were 
immediately  recognized  by  the  FAA  as  an  “acceptable  means  of  compliance”  with  the 
requirements  of  FAR  25.305(d).  At  that  time  the  only  requirement  was  that  “The  dynamic 
response  of  the  airplane  to  vertical  and  lateral  continuous  turbulence  must  be  taken  into  account.” 

FAA-ADS-53  provided  the  basis  for  the  current  Appendix  G  to  FAR  25,  and  for  that  matter, 
it  provided  the  basis  for  all  current  continuous  turbulence  criteria  regardless  of  the  certifying 
agency.  The  primary  difference  between  the  criteria  prescribed  in  FAA-ADS-53  and  current 
criteria  ate  in  the  specified  design  gust  velocities  and  their  variation  with  altitude. 

European  civil  regulations  as  of  January  1987  are  specified  in  JAR-25,  Reference  5.  The 
continuous  gust  psd  design  criteria  given  here  are  identical  to  those  given  in  Appendix  G  of 
FAR  25  with  one  major  exception.  JAR  25  makes  no  reference  to  reduced  design  gust  velocities 
for  airplanes  similar  to  those  having  extensive  satisfactory  service  experience  despite  a  lower  gust 
velocity  capability.  This  will  be  discussed  later. 

The  U.S.  Air  Force  and  Navy  are  in  the  process  of  updating  (simplifying)  some  of  the 
military  specifications.  The  Air  Force  Document,  for  example,  is  MIL-A-87221(USAF). 

However,  U.S.  Air  Force  requirements  for  power-spectral  determination  of  limit  design  gust 
loads  are  still  as  given  in  MIL-A-008861A(USAF),  Reference  6.  This  document  was  first  issued 
in  1971.  The  specified  power-spectral  gust  criteria  are  essentially  those  of  FAA-AD5-53  for  a 
mission  analysis  with  design  envelope  floor.  The  scale  of  turbulence  values,  L,  have  been 
reduced  at  altitudes  below  2500  feet  and  the  gust  intensity  parameter  values,  b’s,  have  been 
correspondingly  adjusted.  Provision  is  also  made  for  evaluating  mission  analysis  loads  on  an 
ultimate  as  well  as  on  a  limit  basis.  The  same  basic  criteria  were  published  in  SEG-TR-67-28  in 
1967,  Reference  7.  The  normalized  power  spectrum  specified  for  the  criteria  is  the  Von  Karman 
representation. 

Prior  to  1986  the  US.  Navy  requirements  did  not  require  power-spectral  determination  of  limit 
design  gust  loads.  However,  MIL-A-8861B(A3),  Reference  8,  issued  February  1986  supersedes  MIL- 
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A-8861CASG)  and  specifies  essentially  the  same  requirements  for  power-spectral  determination  of 
limit  design  gust  loads  as  MIL-A -008861 A(USAF). 

Gust  fatigue  and  design  load  psd  requirements  for  both  services  utilize  the  missi  n  analysis 
approach.  Discussion  in  MIL-A-87221(USAF)  does  provide  for  the  use  of  “envelope  imit  gust” 
which  is  similar  to  the  Design  Envelope  criterion  provided  in  FAA-ADS-53. 

Banc  Forms  of  Criteiion 

As  specified  in  Appendix  G  of  FAR  25,  powcr-spcctra!  gust  loads  criteria  are  presented  in  two 
basic  forms;  the  design  envelope  analysis  and  the  mission  (flight  profile)  analysis.  Provision  is 
also  made  for  a  modification  of  the  design  envelope  analysis  (i.c.,  reduced  design  gust  velocities) 
that  considers  the  service  experience  of  cxistiitg  airplanes. 

In  the  past,  two  shapes  of  gust  velocity  psd  were  commonly  used,  the  Von  Karman  and  the 
Dryden.  Currently,  however,  the  Von  Karman  psd  representation  is  specified  by  FAR  25,  JAR- 
25  and  the  U.S.  military  specifications  for  use  in  both  the  design  envelope  analysis  and  the 
mission  analysis. 


The  response  of  the  airplane  in  vertical  and  lateral  turbulence  is  characterized  by  the  response 
parameters  A  and  Nq  .  A  is  defined  as  the  ratio  of  root-mean-square,  rms,  incremental  load  to 
root-mean-square  gust  velocity,  expressed  as: 


A 


o±  at/0|Hyi\(Q)  dQ 

"  WQ>  « 


Eq.  (1) 


Nq  is  defined  as  the  characteristic  frequency  of  response  (the  average  number  of  times  per 
second  that  the  response  crosses  the  value  zero  with  positive  slope)  and  is  expressed  as: 

1 

N  -JLl.JL 
0  2n  Oy  2n 

where, 

0  »  spatial  frequency,  radians/ft. 

0  -  cut-off  frequency,  an  upper  limit  of  integration  used  in  calculations,  chosen  such  that 

the  calculated  integral  adequately  approximates  the  integral  from  zero  to  infinity. 

*w(0)  is  the  Von  Karman  power  spectrum  for  vertical-lateral  gust  pod’s,  plotted  in  Figure  1. 


Qc/0g2|Hy)2+w(g)  dQ 

Qc/0IHyl2^(Q)  dQ 


Eq.  (2) 
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Eq.  (3) 


L  -  scale  of  turbulence 
ow  -  rms  gust  velocity 

Design  Envelope  Criterion 

The  design  envelope  criterion  is  similar  to  past  discrete  gust  criteria  as  well  as  to  current 
limit  design  maneuver  loads  criteria.  Operational  usage  of  the  aircraft  is  ignored.  Instead  the 
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Figure  1.  Von  Karmon  Vertical— Lateral  Gust  PSD"s 

aircraft  response  is  evaluated  for  a  specified  design  envelope  of  speed,  altitude,  gross  weight,  fuel 
weight,  and  center  of  gravity,  c.g.,  position.  For  the  load  (response)  quantities,  y,  that  are  of 
interest,  A’s  are  obtained  by  dynamic  analysis.  The  limit  design  value  of  y  is  given  as: 

Y  design  ’  Vd  '  °")  '  (*  0w>  Id  *  *  (ow  V  *  *  Uo  Eq-  (4) 

where, 

ow  -  design  rms  gust  velocity 

Pj  -  design  ratio  of  peak  to  rms  values 

A  is  from  the  dynamic  analysis 


4-5 


The  quantity  Ufl  is  the  product  of  a  design  rms  gust  velocity  ow  and  a  design  ratio  of  peak 
to  rms  values,  jjd.  U0  has  the  units  of  velocity  and  can  be  thought  of  as  a  continuous 
turbulence  design  gust  velocity.  The  breakdown  between  ow  and  is  shown  here  to  help  in 
visualizing  the  criterion;  only  the  product  U0  is  specified. 

U0  can  also  be  expressed  as  a  design  value  of  y/A;  that  is. 


Design  values  of  U0  arc  specified  as  a  function  of  altitude,  much  like  the  U{jc  values  of 
discrete  gust  velocity.  Ug  ,  however,  is  a  true  gust  velocity  and  is  an  equivalent  gust 

velocity.  Values  of  UQ  ,  at  speed  Vq  ,  are  defined  as  85  fps  true  gust  velocity  from  an 
altitude  of  0  to  30,000  ft  with  a  linear  reduction  to  30  fps  at  an  altitude  of  80,000  ft.  The 

variation  of  design  UQ  with  altitude  at  is  included  in  the  comparison  of  design  velocities 

given  in  Figure  2.  At  speed  Vg,  U0  is  taken  as  1.32  times  the  V^-,  values  and  at  speed  V^, 
0.5  times  the  Vq  values. 

Reduced  U0  Requirement 

Appendix  G  of  FAR  25,  Item  (b)(3)(i),  provides  for  reduced  design  values  of  U0  . 
Specifically,  “Where  the  Administrator  finds  that  a  design  is  comparable  to  a  similar  design  with 
extensive  satisfactory  service  experience,  it  will  be  acceptable  to  select  at  less  than  85 
fps,  but  not  less  than  75  fps,  with  linear  decrease  from  that  value  at  20,000  feet  to  30  fps  at 

80,000  feet"  A  plot  of  the  variation  of  UQ  with  altitude  is  included  in  Figure  2.  To  apply  the 

reduced  UQ  values  requires  that: 

1)  Transfer  functions  of  the  new  design  are  similar  to  the  prior  designs. 

2)  Typical  missions  of  the  new  airplane  are  substantially  equivalent  vO  that  of  the  similar 
design. 

3)  The  similar  design  should  demonstrate  the  adequacy  of  the  U0  selected. 

This  modification  to  the  design  envelope  criterion  came  about  from  an  AIA  proposal  to  the 
FAA  after  extensive  studies  of  mid-range  to  long-range  transports,  such  as  the  L-1011,  DC-9  and 
DC-10,  and  the  Boeing  727,  737,  747,  757,  and  767,  that  showed  UQ  of  75  fps  at  V^,  was 
adequate  under  FAR  25  Appendix  G  for  this  type  of  transport.  The  higher  values  specified 
by  the  basic  design  envelope  criterion  are  more  appropriate  for  the  lower  cruise  altitude  more- 
severe  types  of  operation.  The  more-severe  types  of  operation  are  represented  by  short  range  or 
commuter  operations  where  cruise  altitudes  of  20,000  to  30,000  feet  are  typical.  The  mid-range 
to  long-range  airplanes  normally  have  cruise  altitudes  in  the  vicinity  of  35,000  feet 

The  requirement  of  similar  transfer  functions  to  qualify  for  use  of  the  modified  design 
velocities  does  not  seem  to  be  particularly  relevant.  Operation  of  the  airplane  is  a  much  more 
significant  consideration.  It  is  logical  to  assume  that  the  mid-  to  long-range  transports  could  be 
classified  a  single  type  of  airplane  having  satifactory  service  experience;  new  airplanes  in  this 
catagory  should  then  qualify  for  use  of  the  reduced  U0  design  velocities. 

It  was  noted  earlier  that  JAR-25  makes  no  reference  to  modification  of  the  design 
velocities,  nor  does  FAA-ADS-53,  which  is  referenced  by  JAR-25.  Depending  on  the  methods 
used  by  an  individual  manufacturer  and  the  service  history  of  their  past  airplanes  this  may  or 
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Figure  2.  Variation  of  Design  Ug  with  /iiiifude,  Criferia  Comparison 

may  not  cause  some  hardship.  If  the  design  gust  loads  for  a  specific  airplane  are  obtained  using 
the  mission  analysis  criterion,  later  studies  of  derivative  aircaft  can  easily  be  related  to  the 
original  mission  analysis  in  much  the  same  manner  as  a  design  envelope  ana.ys.s  U.e  by 
comparison  of  predominant  mission  analysis  flight  conditions).  It.  t.ns  nsarce,  lace  trf  the 
reduced  U  design  envelope  criterion  should  not  cause  significant  nobler  a  v.tthods  are  jeing 
proposed  for  JAR-25  evaluation  to  define  a  reduction  or  PSD  gust  intensity  w  lets  '-an 

ft/sec. 


Mission  Analysis  Criterion 

As  originally  developed  in  FAA-ADS-53,  the  mission  analysis  approach  was  a  “stand  alone” 
method.  It  was,  however,  suggested  that  the  most  appropriate  criterion  would  be  a  combination 
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of  the  design  envelope  approach  and  the  mission  analysis  approach.  Combining  these  two 
approaches  now  constitutes  the  Mission  Analysis  Criterion  specified  in  FAR  25,  JAR-25  and  the 
various  U5.  military  specifications. 

In  addition  to  a  mission  profile  analysis  the  Mission  Analysis  Criterion  requires  that  a  design 

envelope  analysis  be  performed  similar  to  the  design  envelope  criterion,  but  with  reduced  U0 

values  to  provide  a  design  envelope  floor.  The  U0  values  at  are  specified  as  60  fps  from 
0  to  ?0,000  feet  with  a  linear  reduction  to  25  fps  between  30,000  feet  and  80,000  feet.  The 
Vg  ana  Vp  values  arc  still  1.32  and  0.5  times  the  V^.  values,  respectively.  The  variation  of 

UQ  with  altitude  for  use  with  the  mission  analysis  is  also  shown  in  Figure  2. 

The  mission  analysis  approach  is  based  on  “Rice’s  Equation”,  first  published  in  Reference  9. 
Rice’s  equation  is: 

N(Y)  -  N0»  *  W  Eq-  (6) 


This  equation  yields  N(y),  the  number  of  crossings  of  a  given  y,  per  unit  time,  with  positive 
slope.  N0  is  the  number  of  zero  crossings  per  unit  time  with  positive  slope. 

Application  of  Rice’s  equation  to  determine  frequency  of  exceedance  as  a  function  of  load 
level  requires  that  the  equation  be  modified  to  accommodate  one  or  more  mission  profiles.  The 
mission  profiles  represent  the  expected  utilization  of  the  airplane.  Each  profile  is  divided  into  a 
number  of  mission  segments  to  account  for  variations  over  the  flight  profile  of  the  various 
parameters  affecting  A  and  Nq  and  the  variation  with  altitude  of  the  expected  ow  exposure. 
The  required  modifications  to  Rice’s  equation  result  in  the  following  expression. 

NIY)  -  X  tN„  [plWp P2«P  (-^r1')  E9-  (7) 

where, 

fraction  of  total  mission  time  in  each  segment 
net  value  of  response  quantity 
value  of  response  quantity  in  one-g  level  flight 
denotes  summation  over  all  mission  segments 
parameters  determined  by  psd  dynamic  analysis 

parameters  defining  the  probability  distributions  of  rms  gust  velocity. 
These  values  are  shown  in  Figures  3  and  4. 

This  expression  provides  frequency  of  exceedance  curves  from  which  the  limit  gust  loads  are 
read  at  a  frequency  of  exceedance  of  2  x  10  ^  exceedances  per  hour.  The  parameters,  P,,  P,, 
b„  and  b  ,  depend  only  on  altitude.  By  setting  y^  equal  to  zero,  the  variable  N(y)/Ng  can  be 
plotted  versus  y/A  to  produce  the  generalized  exceedance  curves  shown  in  Figure  5.  Actual 
exceedance  curves,  to  be  applied  about  the  one-g  flight  load,  can  be  obtained  by  multiplying  the 
ordinates  by  N0,  and  the  abscissas  by  A,  for  the  specific  load  quantity  and  flight  case  desired. 


ylg 

X 

A,  N0 

b»»  bj  — 


ALTITUDE,  1000  FT. 


Figure  3.  Parameters,  Pj  and  P^,  for  Defining  the  Probability 
Distributions  of  RMS  Gust  Velocity 
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Figure  5.  Generalized  Exceedance  Curves 
Mission  Profile  Generation 

The  generation  of  realistic  profiles  is  not  a  particularly  straight  forward  procedure.  It  requires 
ingenuity  and  judgment  to  apply  information  from  a  number  of  different  areas  that  includes 
outside  sources.  For  example,  development  of  realistic  profiles  for  a  commercial  transport  requires 
knowledge  of  the  target  airlines’  current  and  anticipated  route  structures,  probable  passenger  load 
factors,  and  cargo  loading  practices. 

Mission  profiles  should  reflect  anticipated  operational  usage  for  flight  parameters  such  as  speed, 
payload,  flight  duration,  c.g.  location,  passenger  load  factor,  etc.  A  number  of  different  profiles 
can  be  generated  to  reflect  significant  variations  in  anticipated  usa'ge  from  one  operator  to  the 
next.  However,  the  profiles  should  still  reflect  the  composite  of  all  operations,  rather  than  the 
most  severe. 

Design  loads  obtained  from  a  mission  profile  analysis  are  not  normally  increased  simply 
because  utilization  by  a  new  operator  is  moderately  more  severe  that  anticipated  in  the  original 
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minion  analysis.  In  most  instances,  the  effect  of  a  single  operator  on  the  frequency  of 
exceedance  of  limit  or  ultimate  load  for  the  fleet  as  a  whole  (all  operations)  is  small  enough 
that  no  change  in  toads  is  required.  It  is  recommended,  however,  that  such  changes  be 
qualitatively  evaluated. 

A  number  of  profiles  are  normally  generated  to  adequately  represent  variations  in  profile 
distance  and  duration,  cruise  altitude  and  Mr  ih  number,  payload  (both  passenger  load  factor  and 
cargo),  fuel,  and  pilot  training  or  check  flights.  Plots  of  a  representative  mission  profile, 
developed  for  the  L-1011-1,  are  shown  in  Figures  6  and  7.  Typical  segmentation  of  that  profile 
is  shown  in  Figure  8. 

Flape-ex tended  segments  contribute  to  wing  fatigue,  flap  fatigue,  and  may  affect  limit 
horizontal  tail  loads,  but  they  do  not  contribute  significantly  to  wing  limit  load  exceedances.  The 
flaps-ex tended  segments  are,  therefore,  not  shown  in  Figures  6  through  8. 

The  parameter  values  to  be  used  in  the  definition  cf  the  mission  profiles  must  be  selected 
with  care.  This  is  illustrated  in  Figure  9.  Two  samples  of  airspeed  data  for  the  L-1011-1  are 
shown.  In  the  upper  sketch  an  airspeed  of  290  knots  based  on  a  simple  average  is  reasonable. 
However,  the  airspeed  range  is  so  large  in  the  lower  sketch  that  a  simple  average  is  not  realistic. 
This  is  demonstrated  by  considering  the  airspeed  to  be  represented  by  high  and  low  speed 
segments  with  the  average  speed  in  each  segment  used  for  analysis.  In  this  representation  the 
contribution  of  the  low  speed  segment  to  the  exceedance  curve  is  negligible  and  the  high  speed 
portion  contributes  half  as  many  cycles  as  the  total  distribution  but  with  a  conridcrably  higher 
average  speed.  The  increase  in  loads  due  to  the  higher  speed  has  a  greater  effect  on  the 
exceedance  curve  than  the  reduction  in  cycles.  In  this  case  a  weighted  average  should  be  used. 


Mission  Analysis  nr  Design  Envelope? 

The  question  arises  as  to  which  form  of  criterion,  the  mission  analysis  or  design  envelope,  is 
most  applicable  for  determination  of  limit  design  gust  loads  of  a  specific  airplane.  The  answer 
is  not  always  dear;  both  have  advantages  and  disadvantages. 

Historically,  before  the  introduction  of  pover-spectral  methods,  the  design  envelope  type  of 
analysis  was  the  most  common.  However,  even  then,  a  mission  analysis  type  of  evaluation  was 

performed  if  there  were  doubts  concerning  either  the  ability  of  a  given  airplane  to  withstand  the 

gust  loads  to  which  it  might  be  exposed  or  the  applicability  of  the  existing  criteria  to  new 

aircraft  with  mission  profiles  that  were  considerably  different  than  prior  aircraft 

Moat  airplanes  operate  well  within  their  placard  speeds.  NASA  VGH  data  on  actual 
operational  usage  of  similar  airplanes  shows  that  the  spread  between  actual  and  placard  speeds  is 
never  fas  than  10  to  IS  knots;  the  gap  is  often  greater.  The  mission  analysis  approach  will 
provide  adequate  loads  for  airplanes  that  operate  dose  to  their  design  envelopes  most  of  the  time 
as  well  as  for  airplanes  that  operate  well  within  their  design  envelopes.  The  design  envelope 
approach  will  either  over  estimate  the  loads  for  airplanes  that  operate  well  within  their  design 
envelopes  or  under  estimate  loads  for  airplanes  that  operate  dose  to  their  design  envelopes. 

On  the  other  hand,  the  antidpated  profiles  that  must  be  defined  to  perform  a  mission 
analysis  evaluation  require  considerable  judgment  This  leads  to  differences  of  opinion  that  can  < 

be  quite  difficult  to  reconcile.  The  profiles  are  then  a  compromise  and  the  loads  obtained  are  j 

never  exactly  “right"  Ufc  of  the  design  envelope  criterion  eliminates  this  type  of  uncertainty.  j 
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Typical  Flight  Profile 
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IN  NAUTICAL  MK.E  FLIGHT  PROFILE 


Vt  -  EQUIVALENT  AIRSPEED  ~  KNOTS 


Ftgun  7.  Typical  FUght  Profile  on  Spmd-altitude  Coordinates 


In  addition,  the  mission  analysis  criterion  requires  that  all  flight  conditions  be  analyzed  and 
the  frequency  of  exceedance  data  generated  before  a  single  limit  design  gust  load  is  obtained. 
Using  the  design  envelope  approach,  loads  can  easily  be  obtained  even  for  preliminary  design  by 
analyzing  a  selected  number  of  flight  conditions  that  are  expected  to  be  critical.  Additional 
conditions  can  be  added  throughout  the  design  process. 

Because  a  mission  analysis  is  normally  required  to  obtain  repeated  load  spectra  for  fatigue 
analysis,  use  of  the  design  envelope  does  not  eliminate  the  need  to 'generate  mission  profiles.  In 
addition,  use  of  the  design  envelope  does  NOT  guarantee  a  conservative  design.  Nor  does  it 
guarantee  that  the  airplane  can  be  safely  operated  at  any  point  within  the  design  envelope. 
Design  gust  velocities  are  based  on  the  satisfactory  performance  of  past  airplanes.  The  design 
gust  velocities  are  adequate  for  a  new  airplane  only  if  it  operates  in  a  manner  similar  to  these 
past  airplanes,  for  example  -  well  within  its  design  envelope. 
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Ftgur *  9.  Typical  Alrspmd  Distributions 


If  the  requirements  of  the  certifying  agency  are  not  the  determining  factor,  the  choice 
between  the  mission  analysis  approach  or  the  design  envelope  approach  is  pretty  much  determined 
by  the  anticipated  operation  of  the  airplane  relative  to  prior  airplanes.  If  it  is  basically  an 
existing  airplane  design,  intended  to  fly  in  the  same  manner  ts  past  airplanes,  the  design 
envelope  criterion  should  be  adequate.  Howevei,  for  all  other  designs,  including  the  use  of  active 
controls,  the  mission  analysis  approach  should  be  seriously  considered. 

In  spite  of  the  difficulties  in  defining  mission  profiles,  the  Mission  Analysis  Criterion  has 
always  been  preferred  at  the  author’s  company.  The  other  U.S.  companies  surveyed  are  applying 
the  Mission  Analysis  Criterion  with  increasing  frequency.  Military  specifications  specifically 
require  use  of  the  mission  analysis  approach  for  gust  critics!  airplanes.  The  Design  Envelope 
Criterion  is  used  primarily  for  designing  derivatives  of  existing  airplanes  or  in  the  design  of 


4-15 


airplanes  that  are  similar  to  past  airplanes  in  both  design  technology  and  anticipated  operation. 
It  is  also  used  extensively  in  preliminary  design  studies. 

DESIGN  APPLICATIONS 

In  compliance  with  either  the  design  envelope  or  mission  analysis  criterion,  limit  design  levels 
of  a  variety  of  airplane  response  quantities  can  be  established.  These  response  quantities  may 
include  not  only  external  loads,  such  as  shears,  bending  moments,  and  torsions,  and  translational 
and  rotational  accelerations,  but  also  internal  loads  and  stresses  acting  on  the  various  structural 
elements.  Design  limit  values  of  stress  in  the  various  structural  elements  may  be  obtained  from 
the  statistical  parameters  either  by  direct  computation  of  the  internal  stresses  or  by  generating 
design  load  conditions  that  produce  the  design  stress  levels  when  applied  to  the  structure. 

In  effect,  the  direct  computation  of  internal  stresses  involves  determining  separate  power 
spectra  for  loads  in  every  element  of  the  structure.  For  an  entire  airplane  the  number  of 
structural  elements  could  number  in  the  thousands.  In  addition,  where  the  strength  of  an 
element  such  as  a  wing  surface  panel  involves  the  interaction  of  two  stresses,  for  example 
oompiessior;  and  shear,  these  stresses  must  be  properly  combined  (phased)  to  provide  the  necessary 
stress  information  for  design.  Because  rms  stresses  obtained  in  this  manner  are  directly  tied  to 
the  structure,  changes  in  the  stress  model  require  that  the  internal  load  rms  values  be 
recomputed.  Application  of  this  approach  has  then  been  limited  to  local  areas  of  structure  with 
relatively  complicated  loading  patterns,  where  the  gust  loads  tend  to  be  critical. 

The  more  common  approach  is  to  generate  design  load  conditions.  A  design  condition  consists 
of  a  set  of  external  forces  in  equilibrium  that  represents  the  statistically  defined  parameters  from 
the  psd  analysis  over  specific  regions  of  the  airplane.  By  applying  such  a  set  of  forces,  the 
stresses  in  every  element  of  the  structure  can  be  determined  and  the  same  set  of  forces  can  be 
applied  in  static  tests.  These  conditions  are  analyzed  by  the  Stress  Department  in  the  same 
manner  as  conditions  obtained  from  t;me  history  or  static  loads  analyses  to  determine  the  internal 
loads  and  stresses  acting  on  the  structure.  Secondary  structured  changes  to  the  stress  model  have 
little  effect  on  the  externally  applied  loads,  and  such  changes  can  be  easily  evaluated. 

The  design  procedures  that  are  discussed  here  are  summarized  below  relative  to  the  approach 
used  -  that  is,  computation  of  internal  stresses  or  development  of  design  conditions. 


DESIGN  PROCEDURE  APPROACH 

Interned  Design 

Stress  Condition 

Matching  Conditions  x 

Conditioned  Probability  Method  x 

Equivalent  Discrete  Gust  x 

Internal  Load  Method  x 


Some  companies  apply  more  than  one  method.  This  is  done  because  some  procedures  require 
modeling  of  the  airplane  that  is  not  available  for  older  derivatives.  In  addition,  some  of  the 
more  sophisticated  methods  are  applied  only  to  the  most  complex  areas  of  the  structure  or  to 
areas  that  are  gust  critical  or  sensitive  to  gust  loading.  Conversely,  the  more  easily  applied  but 
potentially  less  accurate  procedures  are  applied  in  a  conservative  manner  to  non-gust-critical 
structure  with  relatively  simple  loading  patterns. 

All  of  the  airframe  manufacturers  contributing  to  this  study  emphasized  that  the  psd  design 
procedures  require  close  coordination  and  a  high  degree  of  cooperation  between  the  Loads 
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Department  (responsible  for  providing  the  loading  or  stress  environment)  and  the  Stress 
Department  (responsible  for  sizing  the  various  structural  elements).  In  many  respects  it  is  an 
interdisciplinary  analysis. 

With  the  exception  of  the  Equivalent  Discrete  Gust  Approach  all  of  the  above  procedures 
explicitly  consider  the  problem  of  phasing. 

Phasing 

Values  of  the  response  quantities  obtained  by  pad  analysis  are  inherently  unsigned.  Because 
positive  and  negative  values  are  equally  likely  to  occur,  both  must  be  considered  in  establishing 
the  design  loads.  For  a  mission  analysis,  separate  exceedance  curves  are  obtained  for  positive  and 
negative  net  loads  (response  quantities).  For  the  design  envelope  analysis,  the  design  net  values 
are  computed  as: 


Net  load  -  +  U0  A  +  Lj^  Eq.  (8) 

Because  the  design  values  of  the  various  response  quantities  generally  occur  at  different  times, 
the  above  design  values  are  “unphased.”  Design  of  the  structure  cannot  be  determined  until 
proper  combinations  of  these  response  quantities  sue  defined.  This  is  true  even  if  internal  stress 
quantities  are  directly  computed  in  the  psd  analysis. 

The  phase  relationship  for  any  two  response  quantities  is  completely  defined  by  their 
covariance  and  respective  variances,  usually  represented  as  con  elation  coefficients.  The  method  of 
fictitious  structural  elements  also  provides  phasing  information.  The  use  of  correlation  coefficients 
is  the  more  popular  procedure;  however,  fictitious  structural  elements  are  occasionally  used  for 
special  considerations. 

Corrciadon  Coefficients 


In  applying  the  pad  method  it  is  assumed  that  atmospheric  turbulence  is  a  random  Gaussian 
process  which  acts  as  an  input  to  the  airplane.  The  airplane  is  represented  as  a  linear  system. 
The  resulting  outputs  from  the  psd  analysis  are  then  alro  a  Gaussian  process.  For  a  Gaussian 
process  the  statistical  dependency  of  any  two  response  quantities  (outputs)  can  be  represented  by  a 
correlation  coefficient,  p  ,  given  by  the  following  expression,  which  follows  from  the  equation? 
given  in  Appendix  B  of  Reference  4. 
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Eq.  (9> 


-  forcing  frequencies,  the  array  of  frequencies  for  which  the  psd  analysis  was 
performed 

-  gust  spectrum  at  each  forcing  frequency 

-  Ratio  of  the  rms  of  the  response  quantity  to  the  rms  value  of  the  gust 
velocity,  obtained  from  the  psa  analysis 

-  characteristic  frequency  of  the  response  quantity,  obtained  from  the  psd 
analysis 

-  complex  transfer  function  of  the  response  quantity  at  each  forcing  frequency 


Eq.  (10) 


TF(w)  -  H(w)  ♦  i  H(ti>) 

rest  imsg 


4-17 


The  phasing  of  any  two  response  quantities  is  then  determined  ty  using  the  correlation 
coefficients  in  the  expression  for  a  joint  probability  density  function  given  in  Reference  4,  that 


Eq.  (11) 


By  assigning  various  constant  values  to  P(x,y),  for  a  given  pv„  ,  contours  of  constant  joint 
probability  density  are  defined.  These  are  ellipses,  collapsing  to  a  straight  line  at  p  -  1  or  a 
circle  at  p  -  0.  An  ellipse  can  be  defined  that  is  tangent  to  the  design  level  values  of  any 
two  specified  response  quantities.  This  is  referred  to  as  a  design  ellipse  and  gives  a  complete 
representation  of  the  phased  design  loads  for  these  two  response  quantities,  Figure  10.  By 
considering  only  the  incremental  loads  and  normalizing  these  loads  to  their  design  values,  the 
variation  of  the  probability  density  function  (ie,  different  values  of  px y)  at  constant  stress  levels 
can  be  illustrated  as  shown  in  Figure  11. 


Fictitious  Structural  Element 


Before  it  was  recognized  that  proper  load  combinations  could  be  easily  obtained  through  the 
use  of  correlation  coefficients  to  define  an  equal-probability  ellipse,  load  phasing  was  accomplished 
using  the  concept  of  the  fictitious  structural  element,  Reference  3.  This  method  still  has  some 
potential  advantages  relative  to  the  correlation  coefficient  approach. 

The  fictitious  structural  element  can  directly  provide  design  level  values  of  combined  internal 
loads  and  stresses.  For  example,  by  expressing  front  or  rear  beam  shear  flow  as 

FE  «  a,  S  +  a,  M  +  a,  T  Eq.  (12) 
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Figure  11.  Variation  of  Probability  Density  Function,  P(xy), 
Maintaining  Constant  Stress  Levels 


the  frequency  response  function  of  " FE "  can  be  determined  as  a  linear  function  of  the 
frequency  response  functions  of  S,  M,  and  T.  Its  psd,  A,  N0  ,  and  design  value  can  then  be 
determined.  To  further  illustrate,  assume  that  shear  flow  is  affected  only  by  shear  and  torsion. 
The  above  equation  then  becomes 


FE  -  a,  S  +  a,  T  Eq.  (13) 

which  defines  a  diagonal  straight  line  on  shear-torsion  coordinates,  as  shown  in  Figure  12.  This 
line  represents  the  shear  and  torsion  that  result  in  the  design  magnitude  of  the  front  beam  shear 
flow.  No  valid  combination  of  shear  and  torsion  can  exceed  this  line.  Therefore,  in  Figure  12, 
points  2  and  3  could  be  considered  realistic  design  load  combinations,  but  point  1  is  obviously 
conservative. 

The  procedure  can  be  applied  to  any  two  response  quantities,  not  just  front  beam  shear  flow. 
In  addition,  arbitrary  values  of  a,  and  a,  can  be  selected.  This  defines  a  family  of  diagonal 
lines  each  representing  the  design  load  or  stiers  in  a  fictitious  structural  element.  For  a  design 
envelope  analysis,  where  the  design  load  is  defined  as  a  constant  times  the  rms  value,  this 
family  of  diagonal  lines  produces  the  same  design  ellipse  that  is  obtained  using  the  correlation 
coefficient  approach. 

Current  psd  gust  procedures  normally  apply  the  correlation  coefficient  method  to  establish 
phasing.  However,  one  potential  advantage  in  using  fictitious  structural  elements  is  that  they  can 
be  carried  through  an  entire  mission  analysis,  thereby  providing  direct  load  phasing  information 
that  includes  the  effects  of  differences  from  segment  to  segment  in  one-g  flight  loads,  correlation 
coefficients,  and  ratios  of  (say)  design  shear  to  torsion.  In  general,  the  result  is  not  an  ellipse. 
This  approach  is  also  easily  applied  to  the  time  history  determination  of  loads  and  is  currently 
used  by  at  the  author’s  company  in  the  determination  of  design  taxi  loads. 
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Figure  12.  Fictitious  Element  Shear—  Tordon  Representation  of  Front  Beam  Shear  Flow 


Matching  Condition  -  Concept 

The  matching  condition  procedure  or  a  variation  is  used  by  all  companies  surveyed  as  one 
method  of  generating  design  conditions.  It  was  first  proposed  in  Reference  3  and  further 
developed  in  References  10  and  11.  In  Reference  12,  linear  optimization  techniques  were  applied 
to  the  procedure,  which  resulted  in  a  batch  processing  computer  program  that  provided  solutions 
to  the  matching  operation.  The  matching  condition  concept  is  described  by  the  following  excerpt 
from  Reference  3. 

"The  basic  concept  employed  in  matching-condition  generation  is  suggested  by  the  fact  that,  in 
flying  through  turbulent  air,  an  airplane  responds  statically  to  the  low  frequency  components  of 
the  turbulence  (long  gradient  gusts)  and  it  responds  dynamically  in  its  various  elastic  modes  to 
the  higher  frequency  components  of  the  turbulence.  The  two  types  of  responses  -  the  static  and 
the  dynamic  -  generally  have  quite  different  distributions  of  load  throughout  the  structure. 
Moreover,  each  elastic  mode  will  have  its  own  distinctive  load  distribution.  In  flight  through 
typical  turbulence  there  is  a  random  interplay  among  these  various  distributions.  As  a  result,  no 
single  distribution  can  be  expected  to  reproduce  simultaneously  the  correct  stress  histories  at  all 
points  in  the  structure.” 

"Accordingly,  in  generating  matching  conditions,  the  approach  is  to  start  with  a  number  of 
'elementary  distributions.’  Each  of  these  consists  of  a  set  of  forces  in  equilibrium,  representing  the 
static  response  of  the  airplane  or  the  dynamic  response  of  a  particular  elastic  mode.  The 
elementary  distributions,  as  building  blocks,  are  superimposed  in  various  proportions  to  give  a 
number  of  design  conditions,  which,  collectively,  envelope  the  statistically  defined  shears,  bending 
moments  and  torsions.” 
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Application  of  this  concept  requires  a  procedure  to  define  each  of  the  following  items: 

•  design  load  combinations  to  be  matched 

•  elementary  distributions 

•  contribution  of  each  elementary  distribution  to  a  design  condition 
Design  Load  Combinations 

A  design  ellipse,  obtained  by  use  of  either  fictitious  elements  or  correlation  coefficients, 
establishes  the  phase  relatiorship  for  any  two  response  quantities.  To  the  extent  that  the  critical 
internal  stress  required  for  design  of  a  structural  element  is  a  linear  function  of  one  or  two 
response  quantities,  the  design  ellipse  provides  sufficient  information  to  define  the  critical  design 
load  combinations. 

An  infinite  number  of  load  combinations  is  required  to  define  each  point  on  a  design  ellipse. 
However,  linear  combinations  of  loads  can  be  easily  defined  that  circumscribe  the  design  ellipse 
and  thereby  provide  a  limited  number  of  load  combinations  that  produce  a  conservative  value  for 
the  design  stress  that  is  a  linear  function  of  two  inputs.  For  the  matching  condition  procedure, 
a  design  octagon  is  defined,  that  circumscribes  the  design  ellipse,  illustrated  in  Figure  10  for  shear 
and  torsion.  Points  A  through  H  are  defined  for  the  incremental  design  values  of  load  from 
the  normalized  equation  that  produced  Figure  11.  The  coordinates  of  points  A  through  H,  Figure 
10,  are  then  dependent  only  on  the  value  of  the  correlation  coefficient.  The  normalized 
coordinates  of  points  A  through  H  are: 
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A  set  of  design  load  conditions  is  obtained  by  matching  the  phased  loads  defined  at  each  of 
the  eight  points  for  a  number  of  design  octagons  defined  throughout  the  structure. 

Prior  discussions  of  this  procedure  given  in  References  3  and  10  thru  12  have  emphasized  the 

use  of  wing  shear-torsion  and  bending-torsion  load  combinations.  Shear  and  bending  at  the 

various  wing  locations  are  in  general  highly  correlated.  The  elementary  distributions  inherently 

reflect  this  correlation  and  the  design  conditions  that  are  developed  based  on  the  shear-torsion  and 

bending-torsion  load  combinations  also  provide  quite  rational  load  combinations  of  shear  and 
bending.  This,  of  course,  is  not  the  case  in  all  areas  of  the  structure. 

» 

A  procedure  is  proposed  in  Reference  13,  the  "equal  probability  technique,”  to  define  design 
conditions  for  load  combinations  that  specifically  account  for  the  effect  of  three  (or  more)  load 
inputs  on  the  design  stress  levels.  For  example,  if  shear,  bending,  and  torsion  all  contribute 
significantly  to  a  design  stress,  the  critical  load  combinations  are  defined  as  a  three  dimensional 
ellipsoid,  rather  than  a  two  dimensional  ellipse.  In  comparing  the  equal  probability  technique 
with  the  matching  condition  technique  it  is  stated  that  the  two  methods  are  equivalent  for 
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combinations  of  two  loads,  but  doubts  were  expressed  concerning  the  validity  of  the  matching 
condition  approach  when  three  or  more  load  inputs  are  significant. 

The  doubts  expressed  appear  to  be  pointed  more  toward  the  definition  of  proper  design  load 
combinations  than  toward  the  fundamental  matching  condition  concept.  Effective  application  of 
the  matching  condition  procedure  is  dependent  on  the  determination  of  design  load  combinations. 
Unfortunately,  procedures  that  have  been  applied  to  define  load  combinations  as  input  to  the 
matching  condition  program  in  the  more  complex  loading  areas  have  not  been  fully  described  in 
past  references. 

In  complex  loading  areas,  such  as  the  interface  between  wing  and  fuselage,  horizontal  or 
vertical  tail  and  fuselage,  or  wing  engine  and  wing,  a  number  of  additional  loads  are  considered 
that  include  selected  internal  stresses  and  the  forces  and  moments  acting  on  concentrated  mass 
items.  In  addition,  some  areas  of  the  structure  require  the  direct  consideration  of  load  phasing 
between  shear  and  bending. 

Information  is  provided  by  the  Stress  Department  to  assist  in  selecting  additional  load 
combinations  that  are  potentially  critical.  The  information  provided  includes: 


•  Stress  results  for  a  number  of  “study"  conditions 

•  Margins  of  safety  for  current  design  conditions 

•  Unit  load  distributions  for  selected  internal  stresses 


•  Estimates  of  relative  significance  of  various  external  loads  in  producing  specific  internal 
stresses 


From  this  information  a  number  of  additional  load  combinations  are  defined.  For  example,  in 
the  interface  area  of  the  wing  and  wing  engine,  load  combinations  are  defined  that  relate  engine 
farces  and  moments  to  each  other  and  to  wing  shears,  bendings  and  torsions.  Load  combinations 
representing  three  external  load  inputs  are  obtained  using  a  concept  similar  to  application  of  a 
design  octagon.  For  example,  conservative  conditions  for  load  combinations  of  shear,  bending  and 
torsion  can  be  defined  from  the  design  octagons  relating  shear-bending,  shear-torsion,  and  bending- 
torsion.  These  conditions  are: 

•  Maximum  shear  with  related  bending  and  related  to  torsion 

•  Maximum  bending  with  related  shear  and  related  torsion 

•  Maximum  torsion  with  related  shear  and  related  bending 

This  produces  24  load  combinations  that  circumscribe  a  design  ellipsoid  in  a  manner  similar  to 
a  design  octagon  circumscribing  a  design  ellipse.  A  sketch  of  the  concept  is  given  in  Figures  13 
and  14.  These  load  combinations  are  slightly  more  conservative  relative  to  the  ellipsoid  than 
those  produced  by  the  design  octagon  relative  to  a  design  ellipse.  After  applying  the  information 
provided  by  the  Stress  Department  very  few  of  these  types  of  load  combinations  are  actually 
required  for  design.  Simultaneous  consideration  of  more  than  three  load  quantities  has  not  been 
necessary.  It  is  quite  possible,  however,  to  define  load  conditions  that  include  directly  phasing 
the  external  loads  with  specific  internal  stresses. 

Elementary  Distributions 

The  mission  analysis  exceedance  curves  are  surveyed  to  determine  the  dominant  profiled)  and 
dominant  flight  eegment(s).  Critical  design  envelope  flight  conditions  are  implicitly  defined. 
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SHEAR 


Figure  13.  Design  Ellipse  3D  Projection  of  Shear-Bending,  Shear-Torsion  and  Bending-Torsion 

Elementary  distributions  are  obtained  for  each  critical  flight  segment  or  design  envelope  condition 
selected.  For  vertical  gust  analysis  these  distributions  include: 

•  Static  aeroelastic  loads  due  to  a  one-g  static  discrete  gust  These  are  effectively  the  loads 

due  to  an  arbitrary  angle  of  attack,  airloads  balanced  by  plunge  and  pitch  inertia 

•  Static  aeroelastic  (or  rigid)  loads  due  to  a  unit  pitch  rate 

•  Loads  due  to  unit  inboard  and  outboard  aileron  angles,  balanced  by  inertia 

•  Loads  due  to  unit  generalized  elastic  mode  acceleration,  q  ,  combined  with  airloads 
occuring  at  the  associated  modal  displacement,  q 

The  predominant  contribution  of  an  elastic  mode  to  airplane  response  occurs  at  its  resonant 
frequency,  so  the  airloads  associated  with  q,  calculated  on  a  zero  frequency  basis,  are  multiplied 
by  the  real  part  of  the  lift  growth  function  at  the  resonant  frequency.  The  loads  per  q  are 
generally  not  in  equilibrium  and  are  balanced  by  use  of  the  rigid  airplane  plunge  and  pitch 


inertia.  The  aerodynamic  term,  proportional  to  modal  displacement,  is  then  directly  added  to  the 
inertia  load,  proportional  to  modal  acceleration,  to  form  the  combined  distribution. 


In  the  lateral  gust  analysis  the  following  elementary  distributions  are  used. 


•  Static  aeroelastic  loads  per  unit  sideslip. 

•  Static  ae  oelastic  loads  per  unit  yaw  rate. 

•  Static  aeroelastic  loads  per  unit  roll  rate. 

a  Static  aeroelastic  loads  per  unit  rudder  angle. 

•  Loads  due  to  unit  generalized  elastic  mode  acceleration.  The  airloads  due  to  modal 
displacement  are  neglected. 


LIST  OF  COORDINATES 


POINT 

SHEAR 

BENDING 

TORSION 

1 

1.0 

ib 

i  a 

2 

1.0 

ib 

ib 

3 

1.0 

ia 

i  a 

4 

1.0 

ia 

ib 

S 

i  b 

1.0 

i  a 

6 

ib 

1.0 

ib 

7 

i  ■ 

1.0 

ia 

8 

i  a 

1.0 

ib  , 

9 

ib 

ib 

1.0 

10 

ib 

i  a 

1.0 

11 

i  a 

ib 

1.0 

12 

ia 

ia 

1.0 

Figure  14.  Coordinates  of  Load  Combinations  That  Circumscribe  •  Design  Ellipsoid 


Bach  of  the  elementary  diatributioni  noted  above  ia  calculated  for  each  of  the  critical  mission 
legmen ta  or  design  envelope  conditiona  selected.  They  are  formed  initially  at  panel  loads,  then 
integrated  to  obtain  the  appropriate  loads  or  response  quantities. 

Design  conditions  should  reflect  rational  levels  of  each  contributing  elementary  distribution. 
For  a  design  envelope  condition  the  maximum  realistic  amount  of  a  modal  distribution  is  defined 
by  the  A  value  for  that  distribution  times  Ufl  .  For  a  predominant  mission  segment  the 
maximum  value  is  set  equal  to  the  value  read  from  the  frequency  of  exceedance  curve  for  that 
distribution.  This  is,  however,  an  approximation  and  some  latitude  is  allowed. 

Generation  of  Matching  Conditions 

Each  design  condition  consists  of  a  portion  a,  of  distribution  E„  a,  of  distribution  E,,  etc. 
A  set  of  coefficients  {a.}  defines  both  a  complete  set  of  panel  loads  and  integrated  loads  (response 
quantities)  throughout  the  airplane.  The  objective  of  the  matching  procedure  is  then  to  define 
sets  of  coefficients  {a^}  that  match  the  design  load  combinations  previously  discussed.  The 
problem  is  represented  as  one  of  linear  optimization  of  the  form: 


Design 

Load 

.Combination 


Eq.  (14) 


The  above  expression  is  developed  in  Reference  12  and  is  comprised  of  both  equalities  and 
inequalities.  The  inequalities  are  obtained  from  the  equations  that  define  the  design  octagon 
boundry  lines.  Equalities  correspond  to  the  equations  that  define  the  phased  design  load 
combinations  for  which  a  design  condition  is  desired.  Constraints  are  applied  to  the  problem 
such  that: 

•  All  loads  will  lie  within  their  design  octagon  boundary  lines 

•  The  allowable  magnitude  of  any  elementary  distribution  is  limited  based  on  psd  results 

•  Each  solution,  set  of  (aj)  coefficients,  contains  the  minimum  number  of  elementary 
distributions  with  the  smallest  contribution  possible  to  obtain  a  design  condition. 

Conditional  Probability  Method 

The  Conditional  Probability  Method  is  a  typical  variation  of  the  Matching  Condition  procedure 
that  is  applied  by  one  of  the  US.  manufacturers. 

Load  combinations  for  a  number  of  response  quantities  are  defined  based  on  the  statistical 
dependency  relationships.  These  load  combinations  are  then  matched  by  manipulating  "generic’’ 
external  force  distributions,  rather  than  elementary  distributions,  to  produce  design  conditions  that 
are  applied  to  stress  models  to  develop  internal  loads  or  to  work  with  transformation  matrices 
that  relate  internal  to  external  loads. 

PSD  gust  analysis  flight  conditions  are  selected  based  on  the  criterion  applied.  Mission  Analysis 
or  Design  Envelope.  The  phased  loads  are  defined  as  the  expected  values  of  load  distribution 
that  will  occur  when  a  selected  envelope  load  is  a  maximum.  These  distributions  are  defined 
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on  the  conditional  probability  density  function.  Under  the  condition  that  y  takes  on  the 


value  ym  .  the  ex  pension,  shown  in  Reference  13,  Appendix  A,  is 
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Eq.  (15) 


The  expected  value  of  load  x  given  load  y  is  them 
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Eq.  (16) 


This  is  equivalent  to  selecting  load  combinations  corresponding  to  points  “T”  on  the  design 
ellipse  shown  in  Figure  10. 

When  the  design  envelope  approach  is  used,  ym  represents  a  specific  load  for  a  specific 
configuration  and  flight  condition.  For  the  mission  analysis,  the  loads  reflect  a  weighted  average 
of  all  the  individual  segments  that  compose  the  complete  mission.  Therefore,  the  phasing  formula 
is  modified  by  the  fraction  of  total  mission  time  in  each  segment  to  account  for  the  phasing 
associated  with  each  individual  mission  segment. 

When  a  complete  set  of  external  loads  in  equilibrium  is  required  for  use  with  a  finite 
element  stress  model,  phased  external  loads  are  used  to  define  discrete  force  distributions  for 
application  at  selected  nodes  of  the  model.  The  number  of  node  points  available  to  apply 
external  forces  typically  exceeds  the  number  of  external  loads  calculated  by  PSD  gust  analysis. 
The  algorithms  that  define  the  node  external  forces  assume  generic  distributions  for  the  undefined 
degrees  of  freedom.  Different  distributions  are  normally  assumed  for  the  aerodynamic  as  opposed 
to  the  inertia  forces  in  the  matching  process.  Three  types  of  phased  PSD  loads  are  calculated: 
net  loads,  aerodynamic  loads,  and  inertia  loads.  The  matching  procedure  algorithms  are  subject  to 
the  constraint  that  the  integration  of  the  applied  nodal  forces  must  be  in  equilibrium  ar.d  must 
reproduce  the  phased  PSD  load  combination. 

Equivalent  Discrete  Gust 


One  of  the  older  concepts  of  defining  external  load  distributions  that  represent  psd  results  is 
that  of  an  equivalent  discrete  gust  To  the  best  of  the  author’s  knowledge,  this  method  has  not 
been  used  as  the  primary  design  procedure  by  any  company  in  meeting  continuous  turbulence 
criteria  for  commercial  transports.  This  approach  should  not  be  confused  with  the  Statistical 
Discrete  Gust  (SDG)  concept  that  has  been  developed  by  J.  G.  Jones  as  a  possible  alternative  to 
PSD  analysis. 

The  procedure  is,  basically,  to  define  discrete  gust  conditions  that  match  the  airplane 
translational  and  rotational  accelerations  about  the  center  of  gravity  as  defined  by  PSD  analysis. 
The  resulting  time  history  solutions  are  assumed  to  produce  properly  phased  external  loading 
distributions.  These  conditions  are  then  factored  to  match  specific  PSD  loads  at  various  locations 
on  the  airplane. 

For  a  relatively  “stiff"  airplane,  with  minimal  dynamic  modal  response,  this  procedure  produces 
a  rational  aet  of  external  load  distributions.  However,  for  more  flexible  aircraft,  a  method  of 
defining  the  probable  external  load  pattern,  similar  to  that  discussed  in  Reference  14  for  SDG 
application,  would  seem  to  be  necessary. 


The  procedure  is  used  primarily  during  the  preliminary  design  phase  to  assist  in  defining  those 
Mess  of  the  structure  that  are  potentially  critical  for  gust  loading.  It  can  also  be  used  to 
provide  balanced  external  load  dis’ributions  for  use  in  various  matching  procedures. 

Internal  Load  Method 

As  the  name  implies,  the  Internal  Load  Method  directly  computes  internal  load  responses  in 
the  PSD  gust  analysis.  Unit  load  coefficients  are  defined  that  relate  the  selected  internal  stresses 
to  unit  external  loads  applied  at  each  individual  structural  grid  point  The  procedure  is  really  no 
different  than  the  computation  of  integrated  external  loads.  Usually  two  internal  loads  or  stresses 

are  sufficient  to  size  a  structural  panel  so  the  problem  of  phasing  three  or  more  external  load 

quantities  is  eliminated.  Correlation  coefficients  are  applied  to  define  a  design  ellipse  whirh  can 
then  be  directly  applied  by  the  Stress  Department 

As  was  stated  earlier,  the  number  of  psd’s  required  to  apply  this  method  throughr.it  the 
airplane  is  very  large.  Application  of  the  procedure  is  normally  limited  to  local  areas  oi 
structure  with  complex  loading  patterns  that  are  also  quite  sensitive  to  gust  loading. 

An  internal  load  method  called  the  Joint  Probability  Technique  is  developed  in  Reference  4. 

This  method  is  not  currently  in  use  and  will  not  be  discussed  here. 

Approach 

The  general  approach  to  sizing  structure  from  psd  results  is  basically  that  of  j  pyramid  or 
hierarchy  of  increasing  complexity. 

At  the  lowest  level,  areas  of  the  structure  that  are  subjected  to  potentially  critical  gust  loading 
are  identified  by  simply  comparing  the  loads  or  stresses  obtained  from  the  psd  analysis  directly 
to  design  envelopes  obtained  for  other  types  of  load  conditions,  such  as  maneuver,  static  and 
diaciete  gust,  or  dynamic  landing  and  taxi  conditions.  Because  of  differences  in  fuselage 
pressurization,  care  must  be  taken  when  comparing  psd  results  with  ground  load  conditions.  The 
maximum  design  values  of  each  load  or  stress  are  applied  without  regard  to  phasing,  only  the 
design  values  are  used.  This  comparison  will  establish  both  areas  of  the  structure  and  the 
quadrant  in  which  the  gust  loading  is  potentially  critical.  The  procedures  discussed  above  arc 
then  applied  to  the  extent  required  to  define  the  design  stress  levels. 

Of  these  procedures,  Matching  Conditions,  or  a  variation  such  as  the  Conditional  Probability,  is 
applied  by  all  of  the  U.S.  manufacturers  that  were  surveyed.  The  Equivalent  Discrete  Gust 
approach  is  used  primarily  as  a  tool  to  provide  additional  information  on  areas  of  the  structure 
that  are  potentially  sensitive  to  gust  loading  and  to  provide  external  foice  distributions  for 
“matching.1'  The  Internal  Load  Method  or  the  use  of  three  dimensional  load  combinations  in  the 
Matching  Condition  approach  are  applied  only  to  areas  of  the  structure  with  complex  internal 
load  patterns  for  which  gust  loads  are  actually  the  critical  design  condition. 

ACTIVE  CONTROLS  CONSIDERATIONS 

The  active  control  system  (ACS)  developed  for  the  L-1011  extended  span  configuration 
involved  symmetric  motion  of  the  outboard  ailerons.  This  was  phased  wim  eg  acceleration  to 
relieve  static  gust  and  maneuver  loads  and  with  wing  tip  velocity  to  increase  elastic  mode 
damping.  Installation  of  the  ACS  was  intended  to  offset  the  load  increase,  due  to  the  increase  in 
wing  span,  in  order  to  minimize  the  extent  of  structural  changes. 
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Load*  due  to  corrective  roll  control  in  turbulence  were  explicitly  included  in  the  dynamic  gust 
analysis  for  the  L-1011  airplane  configuration  with  active  controls.  The  following  background 
information  contributed  to  this  decision. 

Background 

Analysis  of  flight  test  data  from  the  L-1011  development  flight  test  program  in  1971 
indicated  that  the  measured  wing  torsions  were  much  greater  than  predicted  by  theory.  The 
ratio  of  measured  to  predicted  torsions  was  largest  Just  inboard  of  the  outboard  aileron  and  just 
inboard  of  the  inboard  aileron.  At  these  locations  the  measured  torsions  were  l.S  to  3.0  times 
the  theoretical  values-  Bending  moments  and  shears,  however,  were  more  in  line  with  theory. 
The  increased  torsions  were  not  coherent  with  the  measured  gust  velocity,  which  suggested  the 
presence  of  other  inputs.  Although  the  first  wing  antisymmetric  bending  mode  appeared  to  be  a 
contributor,  high  coherencies  between  wing  torsion  and  aileron  angle  indicated  that  corrective  roll 
control  was  also  a  major  source  of  the  increase  in  torsions. 

Fortunately  there  was  sufficient  strength  available  in  the  basic  L-1011  airplane  to 
accommodate  the  measured  torsions.  The  effect  of  the  increased  torsions  was  thereafter  included 
in  all  derivatives  of  the  L-1011  up  to  the  development  of  the  active  controls  system  by 
applying  an  empirical  "flight  test  torsion  increment".  Because  the  active  controls  aircraft 

represented  a  significant  departure  from  prior  derivatives,  a  rational  approach  to  account  explicitly 

for  the  loads  due  to  corrective  roll  control  in  turbulence  was  found  to  be  necessary. 

Flight  tests  were  conducted  during  1977  and  1978,  using  the  L-1011  flight  test  airplane,  as 
part  of  a  NASA-Lockheed  funded  program  to  evaluate  the  use  of  active  controls  for  loads 
reduction.  These  tests  were  made  first  on  the  baseline  span  airplane  without  active  controls  and 
then  with  the  span  increased  9  feet  by  means  of  wing  tip  extensions.  The  results  of  these  tests, 
reported  in  Reference  13,  showed  that  the  active  controls  greatly  reduced  the  loads  due  to  gust, 
Lq,  by  as  much  as  50%  in  the  outer  wing.  But  they  had  no  effect  on  the  loads  due  to  roll 

control.  As  a  result  the  roll  control  loads  became  much  more  conspicuous.  Bending  moments 

and  shears  as  well  as  torsions  were  seen  to  be  involved. 

Loads  due  to  corrective  roll  control  have  always  been  present  when  roll  control  is  by  means 
of  the  outboard  ailerons.  In  the  past,  these  loads  have  been  small  relative  to  loads  produced 
directly  by  the  vertical  gusts.  As  a  result,  they  have  been  adequately  provided  for  by  an 
implicit  conservatism  in  the  design  gust  velocities  (or,  in  a  mission  analysis  determination  of  gust 
loads,  in  the  design  frequency  of  exceedance).  But,  as  noted  above,  an  active-controls  wing  load 
alleviation  system  substantially  reduces  the  vertical  gust  loads,  while  leaving  the  corrective  roll 
control  loads  unchanged.  Consequently,  the  percentage  effect  on  loads  of  the  corrective  roll 
control  increases. 

Further,  this  increase  is  accentuated  by  the  fact  that  the  two  loadings  -  gust  and  roll  control 
-  add  on  a  root-sum-square  basis.  For  example,  if  roll  control  leads,  without  active  controls  are 
30  percent  of  the  gust  loads  then; 
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Although  the  roll  oontrol  load  is  30%  of  the  vertical  gust  load,  the  combined  load  is  Increased 
by  only  12  percent  If  active  controls  are  used  and  wing  vertical  gust  loads  are  reduced  by 
30%  the  result  is: 

With  Active  Controls 


Loads  due  to  vertical  gust 
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05 

Loads  due  to  roll  control 
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05 

Combined  loads 
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The  increase  with  active  controls  is  substantial,  relative  to  Lq  increases  by  a  factor  of 

2.0,  but  the  increase  in  increment  in  net  load  due  to  roll  control  ia  a  much  larger  factor  of 
0414/0118  -  051. 

From  this  example  it  is  clear  that  the  increase  in  gust  loads  due  to  corrective  roll  control 
must  be  explicitly  considered  for  any  airplane  that  accomplishes  roll  control  by  means  of  the 
outboard  ailerons  and  for  which  an  active  control  system  is  used  to  reduce  wing  loads.  In  the 
more  general  context  it  can  be  stated  that  an  active  control  system  can  significantly  alter  the 
contribution  from  vertical  or  lateral  gust  to  combined  loads.  Special  procedures  may  then  be 
required  to  accurately  predict  the  resulting  gust  design  load. 

The  concepts  and  analysis  methods  that  were  used  to  include  corrective  roll  control  effects  in 
the  gust  design  wing  loads  for  the  H01 1  Tristar  with  active  controls  are  summarized  below. 
Based  on  all  available  information,  which  includes  operational  flight  data,  flight  test  data,  flight 
simulator  data,  and  additional  theoretical  studies,  they  provided  a  realistic  set  of  design  loads. 

Other  than  for  the  L-1011,  no  attempt  was  made  in  these  studies  to  explicitly  determine 
when,  or  whether,  or  to  what  extent  roll  control  effects  should  be  included  in  a  dynamic  gust 
loads  analysis.  Therefore,  this  discussion  is  presented  primarily  to  increase  awareness  of  the 
subtleties  involved  in  the  use  of  an  active  controls  system  on  the  determination  of  design  gust 
loads. 

General  Approach 

Design  gust  loads  for  the  basic  L-1011  were  obtained  using  the  mission  analysis  approach. 
The  procedure  to  explicitly  include  loads  due  to  corrective  roll  control  for  the  active  controls 
configuration  was  then  defined  for  mission  analysis  application. 

Results  of  the  basic  L-1011  mision  analysis  showed  that  the  cruise  segments  were  the  greatest 
contributor  to  wing  loading.  This  was  expected,  since  cruise  represents  approximately  80  percent 
of  the  time  in  flight.  In  addition  all  of  the  cruise  segments  exhibit  similar  wing  response 
characteristics.  The  cruise  segment  that  contributed  most  heavily  to  the  wing  shear  and  bending 
loads  at  critical  wing  locations  was  then  selected  for  explicit  roll  control  analysis.  This  is 
referred  to  as  the  predominant  segment  The  effect  of  other  segments,  climb  and  descent,  was 
estimated  using  a  parametric  procedure  and  simplified  exceedance  curves.  Factors  were  developed 
to  adjust  shear,  bending  and  torsion  at  a  number  of  locations  over  the  wing  span.  These  loads 
were  phased  using  a  combined  correlation  coefficient  method,  and  design  load  conditions  were 
developed  for  stress  analysis  using  the  matching  condition  process. 
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3-D  GUST  ANALYSIS 

Extensive  use  was  made  of  a  three-dimensional  gust  analysis  computer  program  that  utilizes 
the  work,  of  Dr.  Frederick  D.  Eichenbaim  of  Lockheed-Georgia  Company,  References  16,  17  and 
18.  An  existing  Lockheed-California  one-dimensional  program  was  modified  to  include  the 
computation  of  gust  input  cross  spectra  and  computation  of  output  spectra  utilizing  these  input 
cross  spectra  and  the  computed  transfer  functions. 

The  term  “three-di  meiu>  onal"  refers  to  the  number  of  position  coordinates  upon  which  the  gust 
velocity  is  assumed  to  depend.  The  most  important  variation  of  the  gust  velocity  is  along  the 
flight  path.  This  is  normally  the  only  variation  considered,  which  results  in  a  "one-dimensional” 
gust  analysis.  The  three-dimensional  analysis  considers  also,  on  a  statistical  basis,  the  s panwise 
variation  of  vertical  gust  velocity.  It  also  considers,  although  these  are  less  important,  the 
vertical  variation  of  the  vertical  gust  velocity,  and  the  vertical  and  lateral  variations  of  the 
lateral  gust  velocity.  In  addition,  it  combines  vertical  and  lateral  inputs  into  a  single  analysis. 

Basically,  the  three-dimensional  gust  analysis  consists  of  three  steps: 


1.  Determination  of  transfer  functions  relating  the  various  airplane  loads  to  gust  velocities 
acting  on  specific  streamwise  gust  strips. 

2.  Determination  of  the  power  spectra  of  gust  velocity  on  the  various  strips  and  cross 
spectra  cf  gust  velocity  for  all  strip  pairs.  The  gust  power  spectrum  on  any  individual 
strip  is  simply  the  usual  one-dimensional  gust  power  spectrum.  Three-dimensional  effects 
arr  brought  in  by  the  cross  spectra. 

3.  The  t?  isfer  functions  and  the  gust  velocity  spectra  and  cross  spectra  are  combined  to 
determine  response  psd’s,  cross  spectra  between  pairs  of  response  quantities,  and  cross 
transfer  functions  that  relate  various  responses  to  the  gust  velocity  at  a  gust  probe,  for 
comparisons  with  flight-measured  transfer  functions  and  coherencies.  These  are  then 
used  as  in  a  one-dimensional  analysis  to  provide  A’s,  No’s  and  correlation  coefficients. 


Steps  2  and  3  were  accomplished  using  Eichenbaum’s  basic  equations  as  presented  in  Reference 
18.  The  only  simplification  was  to  drop  the  fore-aft  component  of  turbulence.  (Although  the 
fore-aft  component  may  become  significant  in  landing  approach,  it  has  only  a  small  effect  at 
the  higher  airplane  speeds  that  produce  the  critical  gust  loads.)  These  equations  not  only  retained 
provision  for  variation  of  vertical  and  lateral  gust  velocities  in  the  y  and  z  directions,  but  also 
provided  for  orientation  of  the  individual  lifting  surface  segments. 

The  program  allowed  20  gust  input  strip;  for  a  half-airplane  (15  were  used)  and  included 

provision  for  arbitrary  dihedral  angles  and  arbitrary  locations  in  the  y-z  plane.  The  gust 

velocity  was  considered  to  vary  linearly  from  a  maximum  at  the  strip  centerline  to  zero  at  the 

adjacent  strip  centerlines.  Computation  of  the  various  gust  velocity  cross  spectra  was  facilitated 
by  applying  the  planar  and  nonplanar  coherencies  tabulated  in  Reference  18.  The  arrangement 
of  the  gust  strips,  defined  for  the  L-1011  studies,  is  shown  in  Figure  15. 

The  three-dimensional  gust  analysis  was  applied  to  determine  the  effect  of  roll  control  by  the 
autopilot  It  was  also  used  to  help  select  the  appropriate  autopilot  mode.  (Some  pilots  prefer  a 
variation  of  the  CWS  (control-wheel-status)  such  as  attitude  hold,  rather  than  the  turbulence 
mode,  when  in  turbulence;  for  the  L-1011  the  roll  control  rms  aileron  angle  values,  Og ’s  were 

about  10%  higher  in  CWS  status  than  in  the  turbulence  modej  * 
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Loads  obtained  from  the  3-D  gust  analysis  were  somewhat  lower  than  the  corresponding  one-D 
loads.  It  was  therefore  necessary  to  design  to  slightly  higher  gust  velocities  in  order  to  retain 
the  concept  of  equivalent  strength  used  in  setting  the  gust  criteria  design  velocities  which  were 
determined  from  a  one-D  gust  analysis.  Based  on  a  comparison  of  several  L-1011  cases,  it  was 

concluded  that  the  3-D  loads  should  be  increased  by  a  factor  of  1.07.  It  is  noted  that  this 

factor  has  much  broader  applicability  than  just  to  the  determination  of  corrective  roll  effects.  If 
3-D  analysis  were  to  become  routine  for  gust  loads  determination,  the  value  of  this  factor  would 
be  of  primary  importance. 

In  the  original  design  study,  3-D  analysis  was  performed  for  a  predominant  mission  analysis 
cruise  segment  with  an  altitude  of  32,000  feet  and  Mach  Number  of  .85.  Later  studies  included 
additional  cruise  segments  and  some  climb  and  descent  segments.  Studies  were  obtained  to 
determine  the  effect  of  ACS-on  and  off  with  and  without  the  autopilot,  the  effect  of  different 
autopilot  modes,  the  effect  of  dihedral,  and  the  effect  of  the  number  of  structural  airframe 
modes  on  the  computation  of  stability  derivatives. 

RMS  Aileron  Angle 

The  parameters  required  to  determine  loads  due  to  both  autopilot  and  pilot  roll  control  are  the 

rms  aileron  angle,  the  aileron  angle  pad,  and,  for  each  response  quantity  to  be  analyzed,  the  ratio 

of  the  roll  control  load  component  to  the  rms  gust  velocity,  A^. 

The  magnitude  of  corrective  roll  control  aileron  angles  in  turbulence  depends  upon  the 
turbulence  intensity  and  the  flight  condition.  It  is  reasonable  to  expect  the  rms  (root-mean- 
square)  aileron  angle,  o^,  to  vary  in  proportion  to  rms  gust  velocity.  A  hey  parameter  in  the 
analysis  is  then  the  rms  aileron  angle  as  a  ratio  to  rms  gust  velocity,  or  A  ,  which  can  be 
expected  to  vary  with  flight  condition. 
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Figure  15.  Gust  Strip  '"Centerlines''  Used  in  3-D  Gust  Analysis 
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In  addition  the  amount  of  corrective  roll  control  is  likely  to  vary  from  pilot  to  pilot, 
between  pilot  and  autopilot,  and  from  one  turbulence  encounter  to  another  for  the  same  pilot. 
As  a  result,  selection  of  a  design  Ag  involves  judgment  in  evaluating  all  of  the  available  data. 

Data  on  the  L-1011  rms  aileron  angle  were  obtained  from  a  number  of  sources. 

•  L-1011  gust  response  flight  tests  and  flight  simulator  tests 

•  L-1011  operational  data 

•  L-1011  3-dimensional  gust  analysis  with  autopilot  roll  control 

Data  from  three  flight  test  programs  were  available.  Seven  bursts  of  from  one  to  four 
minutes  duration  were  available  from  the  original  L-1011  flight  test  program  conducted  in  1971. 
Three  bursts  were  available  from  the  1977  program  and  five  from  the  1978  program.  Roll 
control  was  by  the  pilot,  except  for  a  limited  amount  of  data  from  the  1978  tests. 

Flight  simulator  data  were  available  for  three  different  pilots.  These  data  were  taken  as 
averages  of  six  5-minute  samples  for  pilot  No.  1  and  four  each  for  the  other  two  pilots.  This 
program  was  conducted  in  1978  utilizing  the  1977  flight  test  conditions. 

Operational  data  were  obtained  from  British  European  Airways  (BEA)  AIDS  tapes.  For  the 
original  analysis  only  five  samples  were  available,  all  from  the  same  flight.  The  data  were 
questionable  and  the  analysis  was  quite  crude.  However,  for  later  studies,  50  additional  records 
were  available  from  British  Airways  and  Gulf  Airways,  also  obtained  by  means  of  the  AIDS 
system.  Twenty-two  of  these  records  were  analyzed  for  the  roll  control  study.  As  a  condition 
of  CAA  certification,  British  Airways  continued  to  obtain  such  data,  of  which  17  records  were 
reduced  for  use  in  the  roll  control  analysis. 


Tabulated  values  from  the  AIDS  system  included  CG  acceleration  with  a  sampling  rate  of  8 
per  second  and  aileron  angle  once  per  second.  Rms  values  of  these  quantities  were  calculated 
from  the  time  histories.  Rms  gust  velocity  was  then  determined  from  the  rms  CG  acceleration 
by  means  of  a  procedure  utilizing  an  extensive  set  of  curves  developed  by  Lockheed  in  prior 
studies,  in  which  the  continuous  turbulence  gust  response  factor,  Ka  ,  was  computed  and  plotted 
as  a  function  of  four  dimensionless  rigid  airplane  parameters. 

These  curves  are  based  on  simple  theory.  The  airplane  is  considered  rigid,  but  is  allowed  to 
pitch  as  well  as  plunge.  The  effect  of  gust  penetration  on  pitch  is  neglected,  and  unsteady  lift 
growth  is  accounted  for  only  with  respect  to  the  gust  input,  not  the  airplane  motions.  The  Von 
Karman  shape  of  gust  pad  is  assumed. 


The  curves  were  adjusted  to  agree  with  results  given  by  the  more  sophisticated  analytical 
methods  used  for  the  L-1011  by  back-figuring  K0  for  a  number  of  mission  analysis  flight 
segments.  It  was  found  that  applying  a  factor  of  1.13  to  the  plunge-only  curve  developed  from 
simple  theory  provided  a  good  representation  for  the  L-1011,  Figure  16.  The  relation  of  oAn  to 
ow  is  then  given  as: 
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Eq.  (17) 
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Three-dimensional  theoretical  data  were  obtained  for  evaluation  of  the  autopilot  primarily  from 
the  predominant  mission  segment.  Because  the  autopilot  is  expected  to  do  essentially  the  same  job 
as  the  pilot,  it  should  require  approximately  the  same  amount  of  aileron  as  the  pilot  to  do  it. 
The  pilot  response,  however,  is  not  as  likely  to  be  Gaussian. 

In  all  cases,  the  rms  aileron  angle  was  computed  and  related  to  the  associated  rms  gust 
velocity.  Both  the  rms  aileron  angle  and  gust  velocity  were  adjusted  to  give  “effective"  values 
corresponding  to  a  reference  flight  condition. 

The  use  of  an  effective  value  allows  direct  comparison  of  data  from  various  flight  conditions. 
It  is  based  on  the  concept  that  as  the  flight  condition  and  turbulence  vary,  the  roll  control  5 

A 

applied  by  the  pilot  will  be  such  as  to  give  a  rolling  moment  proportional  to  the  rolling 
moment  caused  by  the  turbulence.  It  is  believed  that  the  gust  induced  rolling  moments  are 
primarily  due  to  lateral  gust  working  through  the  rolling  moment  due  to  sideslip  stability 
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derivative,  Ct^.  This  parameter  is  then  a  key  value  in  obtaining  the  effective  values  for  o^. 
The  reference  flight  condition  was  defined  as: 


h  -  sea  level  altitude 

Ve  -  300  Knots 

M  -  045 

W  -  350,000  lb 


which  yields: 


and 


where 


V^o  V.  \  Ko  °w 
°w-«  *  1  300  knots  f,  K„  1 


Eq.  (18) 


Eq.  (19) 


density  of  the  atmosphere,  pQ  -  sea  level  value 
equivalent  airspeed 

rate  of  change  of  rolling  moment  coefficient  wit'  aileron  angle 
rate  of  change  of  rolling  moment  coefficient  with  sideslip  angle 

ratio  of  rms  sideslip  angle  produced  by  continuous  turbulence  to  rms  “sharp-edge” 
gust  sideslip  angle 


The  operational  data  (except  for  the  very  early,  questionable  five  data  points)  is  plotted  in 
Figure  17.  (The  flight  test  and  simulator  data  points  are  not  shown,  but  they  displayed  a 
similar  trend.)  The  solid  line  represents  ^  for  the  autopilot  as  determined  from  the  3-D 
analysis.  The  dashed  line  is  a  least  square  approximation  for  the  pilot  related  data.  The 
effective  value  of  A,  for  the  autopilot  was  0.240  deg/fps  and  for  pilot  control  0.210  deg/fps. 
In  the  original  design  study  0.240  was  used  for  both  pilot  and  autopilot  Actual  values  of  A^ 
for  a  given  flight  condition  are  then  obtained  by  reverse  application  of  the  above  parametric 
relationship. 


Aileron  Angle  PSD  and  ARC 

Loads  due  to  autopilot  roll  control  could  be  obtained  directly  from  the  3-D  analysis. 
However,  for  pilot  roll  control,  it  was  necessary  to  establish  an  aileron  angle  pad  for  computation 
of  the  roll-control  loads.  This  was  done  by  examining  the  ped  shapes  obtained  from  various  test 
flights,  flight  simulator  tests,  and  theoretical  (autopilot)  cases.  Four  typical  psd  shapes  were 
selected,  shown  in  Figure  18.  Ratios  of  rms  wing  load  to  rms  aileron  angle  for  pertinent  wing 
load  quantities  were  calculated  for  each  psd  shape,  on  the  basis  of  aileron  input  only.  The  two 
most  severe  sets  of  loads  were  then  averaged  to  give  a  set  of  ratios  for  use  in  design. 
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Figure  17.  Operational  Data  on  Relation  of  Effective  RMS  &a  to 
Effective  RMS  Gust  Velocity 


Multiplying  these  ratios  by  A,  yields  the  ratio  of  the  roll  control  load  component  to  rms  gust 
velocity,  ,  that  is: 
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Eq.  (20) 


Loads  Due  to  Roll  Control,  Lg^ 


The  roll  control  component  for  each  load  quantity  can  then  be  calculated  as: 


lRC  _  Uo  *  aRC 


Eq.  (21) 


In  the  original  design  analysis,  roll  control  effects  were  explicitly  obtained  for  one 
predominant  mission  segment  For  this  approach,  was  calculated  seperateiy  for  each  load 
quantity  from  the  expression: 


where 


A  « 


u 


o 


Eq.  (22) 


Limit  design  value  of  the  load  quantity  form  the  mission  analysis  exceedance  curves 
One-g  value  of  the  load  quantity  for  this  mission  segment 
Ratio  of  rms  load  to  rms  gust  velocity  for  this  mission  segment 


4-35 


For  the  predominant  segment,  which  was  selected  based  on  the  contribution  of  wing  shears 
and  bending  moments,  U0  averaged  about  105.  This  means  that  in  order  for  design  gust  loads 
to  be  encountered  during  this  mission  segment  a  design  gust  velocity  of  105  fps  is  required. 
The  Ufl  value  was  then  limited  to  105.  Torsions  generally  required  a  higher  design  gust  velocity 
to  reach  their  design  levels,  indicating  that  this  particular  flight  segment  was  not  critical  for 
torsion. 


Combined  Loads  -  Phasing 

The  combined  loads  due  tc  gust  and  roll  control  were  computed  from  the  following  equation: 

lG  ,  RC  '  NAs  +  LrC  E*  03) 


where  Lq  is  the  statistically  defined  gust  increment  value  Le  -  L^,  factored  by  the  effects  of 
control  system  saturation  and  unavailability  (reliability). 


Tnc  effect  of  the  increment  in  load  due  to  roll  control  was  included  in  the  phasing  by 
adjusting  the  expression  for  computing  correlation  coefficients  to  account  for  a  second  unccrrelated 
input  For  any  two  load  quantities  the  expression  becomes: 
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Eq.  (24) 
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Figure  IS.  Faired  PS  Us  of  Roll  Control  Aileron  Angle 
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where 

p  -  correlation  coefficient 

o  -  A  ,  rms  value 

x,  y,  u  and  v  denote  load  quantity  and  input  source,  for  example,  correlations  for  shear  or 
bending  to  torsion  are  defined  as  shown  in  the  following  table.  However,  the  expression  is 
valid  for  any  two  load  quantities,  not  just  shear  or  bending  to  torsion. 

Load  Quantity 

S  or  B  T 

Input  Gust  x  y 

RC  u  v 

The  effect  of  roll  control  on  the  correlation  of  wing  bending  and  torsion  is  illustrated  in 
Figure  19.  The  effect  on  the  correlation  of  wing  shear  and  torsion  was  similar.  The  correlation 
of  these  leads  tended  to  become  more  negative,  shifting  loads  on  a  design  shear-torsion  or  bending 
torsion  envelope  from  quadrant  I  toward  quadrant  II,  and  broadening  the  load  envelope  in 
general. 

Other  Segments  Factor 

Primary  emphasis  for  determining  the  effect  of  roll  control  was  placed  on  the  predominant 
mission  analysis  segment.  However,  the  climb  segments  are  at  lower  altitude  with  a  higher 
equivalent  airspeed  and  tended  to  yield  a  higher  ratio  of  to  Lq  .  For  these  segments  the 

A^’s  were  estimated  using  the  parametric  approach  as: 


The  effect  on  loads  was  then  determined  by  constructing  simplified  exceedance  curves  for  a 
limited  number  of  segments  within  each  profile.  The  exceedance  curves  for  the  individual 
segments  were  then  adjusted  to  reflect  the  relative  severity  of  the  roll  control  effects.  The  sum 
curves  with  and  without  adjustments  were  then  compared.  The  resulting  factor  was  applied  in 
addition  to  the  Lq  +  ^  /Lq  factor  obtained  based  on  the  predominant  segment 

Other  Considerations 

The  effect  of  control  system  saturation  and  reliability  should  be  accounted  for  in  the  design 
loads.  It  was  noted  earlier  that  this  was  done  by  applying  suitable ’factors.  From  prior  studies 
the  effect  of  system  unavailability  on  loads  was  estimated  as  a  factor  of  1.01. 

Control  system  saturation  is  an  entire  subject  by  itself.  The  methods  used  to  account  for 
saturation  of  the  I.-1011  active  controls  system  are  reported  in  an  AIAA  paper,  Reference  19. 
With  respect  to  roll  control  a  separate  study  indicated  that  saturation  tended  to  reduce  the 
increase  in  loads  due  to  corrective  roll  control.  This  effect  was  not  included  in  the  study. 
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A  BEND) HQ  A  TORSION  P 


OUST  1.00  1.00  0 

ROU  CONTROL  1.00  1.00  -1.00 

COMBINED  1.414  1.414  -  SO 


Figure  19.  Effect  of  Roll  Control  on  Bending-Torsion  Gust  Load  Octagons 


Roll  Control  Accountability  Factor,  RCAF 

Continuous  turbulence  design  gust  criteria  are  based  on  the  concept  of  equivalent  strength  with 
past  airplanes  that  have  a  satisfactory  service  life.  Three  reference  airplanes  were  used  in 
establishing  the  criteria,  the  Lockheed  Model  749  Constellation,  the  Lockheed  Electra  (Model  188), 
and  the  Boeing  Model  720B.  Although  the  Boeing  airplane  limited  the  use  of  outboard  ailerons 
for  roll  control  to  the  flaps  extended  configurations,  the  Lockheed  models  used  outboard  ailerons 
for  roll  control  throughout  their  flight  envelopes.  As  a  result  the  design  gust  velocities,  U0  , 
and  the  2  x  10'^  design  frequency  of  exceedance  include  the  effects  of  roll  control  by  means  of 
the  outboard  aileron  for  aircraft  without  active  controls.  Ther  fore,  the  effect  of  loads  due  to 
corrective  roll  control  is  included  only  to  the  extent  that  the  percent  increase  due  to  roll  control  i 

!  ! 


with  active  controls  exceeds  the  percent  increase  without  active  controls  for  the  seme  airplane, 
that  is 


Design  Load  - 


Load  ACS-on  with  roll  control 
Load  ACS-off  with  roll  control 


Load  ACS-off  without  roll  control 

In  the  above  expression,  the  numerator  represents  various  phased  load  combinations  at  a  given 
location  with  the  full  effect  of  roll  control.  The  denominator  is  referred  to  as  the  "roll  control 
accountability  factor”  or  “RCAF.”  This  is  the  factor  by  which  the  loads  would  increase  due  to 
roll  control  if  the  airplane  were  designed  without  an  active  control  system.  This  effect  was 
included  as  follows. 


First,  gust  loads  were  obtained  with  the  active  control  system  present  and  the  full  effect  of 
roll  control  included.  These  loads  were  then  divided  by  the  "roll  control  accountability  factor.” 
If,  for  example,  roll  control  increased  the  gust  loads  by  66  percent  with  active  controls  and  31 
percent  without,  the  with-active-control  loads  would  be  divided  by  1.31.  In  this  example  the 
RCAF  is  the  factor  IJI. 


The  RCAF  was  evaluated  separately  at  each  wing  station.  It  was  found  to  vary  from  1.01 
at  the  wing  root  to  1.31  at  Buttline  833  (85  percent  semi-span).  This  variation  is  shown  in 
Figure  20.  It  is  seen  that  the  roll  control  effects  are  relatively  small  at  the  root  but  become 
large  near  the  wing  tip. 

Additional  studies  were  later  performed  with  the  intent  of  developing  general  procedures 
separately  for  use  with  mission  analysis  and  design  envelope  analysis  that  eliminated  use  of  the 
RCAF.  The  primary  approach  has  been  to  identify  conservatisms  in  the  methods  and  data  used 
to  obtain  the  loads  due  to  corrective  roll  control.  It  appeared  that,  with  only  fairly  modest 
reduction  in  the  roll  control  loads,  perhaps  in  combination  with  a  very  small  reduction  in  design 
gust  velocities,  the  net  loads  (gust  plus  roll  control)  might  be  sufficiently  reduced  so  that  the 
RCAF  would  not  be  necessary.  Although  these  studies  have  shown  some  promise,  they  have  not 
yet  yielded  a  design  procedure  that  is  ready  for  general  application. 

As  a  result  only  the  original  design  procedure  has  been  presented  here;  more  work  remains  to 
be  done  before  a  more  general  procedure  is  ready  to  use. 


81 

Figure  20.  S panwise  Variation  of  Lq  +  ^  /Lq  for  Wing  Bending 
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SUMMARY  COMMENTS 

Both  the  Mission  Analysis  and  Design  Envelope  criteria  were  developed  from  the  concept  of 
equivalent  strength  with  past  airplanes  that  have  a  satisfactory  service  life.  If  an  airplane  is 
operated  in  a  manner  that  is  considerably  different  than  these  past  airplanes,  the  design  gust 
loads  on  the  new  airplane  will  be  different  depending  on  which  of  the  two  criteria  is  used.  It 
is  important  to  understand  the  differences  in  criteria  to  assure  that  the  criterion  selected  for 
analysis  will  indeed  give  realistic  loads  for  the  new  airplane. 

The  most  common  method  of  determining  internal  loads  and  stresses  for  the  design  and  sizing 
of  structure  from  the  pad  results  is  that  of  "Matching  Conditions”  or  variations  of  this  approach. 
The  selection  of  load  combinations  to  be  matched  and  the  use  of  rational  external  force 
distributions  in  the  development  of  the  design  conditions  are  of  primary  importance.  The  general 
approach  is  to  identify  areas  of  the  structure  that  are  potentially  critical  for  gust  loading  based 
on  a  conservative  comparison  of  the  PSD  gust  loads  with  other  types  of  conditions.  The  more 
complex  analysis  procedures  are  then  applied  to  these  areas  of  the  structure. 

Application  of  the  psd  continuous  gust  design  criteria  to  airplanes  designed  with  active  controls 
and  gust  load  alleviation  systems  presents  a  challenge.  This  is  an  area  where  considerable 
analysis,  development  effort,  and  possibly  new  methods  will  be  required  in  the  coming  years. 
The  use  of  three-dimensional  gust  analysis  procedures  is  likely  to  be  beneficial. 
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Abstract: 

Depending  on  the  country  of  certification  different  gust  models  and  means  of  compliance  of 
the  airworthiness  requirements  have  to  be  covered  in  structural  design  of  civil  transport  aircraft. 

The  influence  on  aircruft  design  from 

-  gust  models 

-  A/C  modelling 

-  control  systems/laws 

is  demonstrated  on  example  of  a  short  to  medium  range  transport  aircraft. 

Recommendations  for  future  harmonised  approaches  in  gust  methodes  and  modelling  will  be  given. 


A  b  b  r  a 

viations 

AA 

Airworthiness  Authority 

A/C 

Aircraft 

CAA 

Civil  Aviation  Authority 

dgac 

Direction  Generale  de  1 'Aviation  civile 

FAA 

Federal  Aviation  Agency 

LBA 

Luftfahrtbundesamt 

RLD 

Rij  ksluchtvaartdienst 

STPA 

Service  Technique  des  Programmes  Aeronautiqui 

FAR  25 

Federal  Airworthiness  Requirements  (US) 

JAR  25 

Joint  Airworthiness  Requirements  (Europe) 

NV 

National  Variant 

AC 

Advisory  Circular 

SC 

Special  Condition 

IM 

Interpretative  Material 

EDP 

Electronic  Data  Processing 

EFCS 

Electrical  Flight  Control  System 

FBW 

Fly  By  Wire 

LAS 

Load  Alleviation  System 

CT 

Continuous  Turbulence 

MA 

Mission  Analysis 

DEA 

Design  Envelope  Analysis 

SDEA 

Supplementary  DEA 

DC 

Discrete  Gust 

DTG 

Discrete  Tuned  Gust 
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MoC  Moans  of  Compliance 

TAS  True  Airspeed 

BAS  Equivalent  Airspeed 

Ude  Derived  Gustvelocity  (fps,  EAS) 

U«  Guatvelocity  (fps,  TAS)  in  CT 

c  Guatgradient 

QF  quasi  flexible 

FF  full  flexible 

NZc  Loadfactor 

Vc  Crusing  speed 

9  rate  of  attitude 

0  Bank  angle 

0c  bank  angle  Indices: 

0c  rate  of  bank  angle  c  commanded 

r  roll  velocity 

6  sideslip  angle 


1.0  Introduction 


Different  from  military  fighters  the  large  transport  aircraft  up  to  relatively  large  grossweights 
are  designed  by  gusts  in  big  parts  of  their  structural  components. 

Being  so,  the  component  weight  and  the  overall  standard  of  safety  during  flying  in  gusts 
will  depend  on 


-  Gust  Models 

-  A/C  Modelling 

-  Systems  Introduction 

used  in  static  design  work.  But  the  manufacturer  is  not  free  in  the  choice  of  models,  neither 
in 


Gusts 


nor  in 


Aircraft 


He  is  rather  guided  by  existing  regulations  or  their  related  and  by  authorities  accepted  means 
of  compliance. 

Although  flying  around  in  the  same  atmosphere  all  over  the  world,  different  countries  require 
through  their  airworthiness  authorities  different  Gust  and  A/C  models. 

AA  normally  defend  their  since  centuries  nearly  unchanged  position  by  stating: 

THE  EXISTING  LEVEL  OF  SAFETY  MUST  NOT  BE  ERODED 

Nobody  will  and  can  oppose  to  that,  the  question  is: 

WHEN  IS  EROSION  STARTING? 

The  author  believes,  based  on  discussions  with  several  AA  in  this  field  that  sometimes  the  improve¬ 
ments  against  the  past  in 


-  more  and  more  sophisticated  investigations  of  gusts  in  relation 
to  masses,  c.G.,  massdistributions,  Mach,  altitude  and  dynamic 
pressure  using  extensively  EDP 

-  physical  and  numerical  methods 

-  correction  of  data  by  ground  and/or  flight  tests 

-  extensive  service  experience  with  former  designs 

-  weather  radar  and  improved  weather  forecast  by  satelites 

-  extensive  use  of  qualified  equipment  and  systems  on  board 
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is  not  adequately  taken  into  account  in  this  judgement. 

Therefore,  while  accepting  the  requirement  for  a  high  safety  standard,  there  is  on  the  other 
hand  the  danger 


TO  PENALIZE  A/C  DESIGNS  BY  UNNECESSARY  CONSERVATISMS 

The  truth  will  be  -  as  always  -  a  very  narrow  path.  Therefore  very  openminded  discussions 
will  be  needed,  to  get  real  progress  in  this  field.  But  let  us  first  review  the  today's  situation. 

3.0  Comparison  of  Requirements  and  Interpretations 


2.1  Gust  Requirements 


The  requirement  situation  for  large  transpoit  /»/C  and  the  different 

interpretations  represented  by  the  following  national  airworthiness  authorities  will  be  reflected 
(Fig  1)  . 


COUNVHY 

REGULATION 

AA 

Germany 

JAR  25 

LBA 

r  nee 

JAR  25 

STPA/DGAC 

Great  Britain 

JAR  25  ■»  NV 

CAA 

NETHEHLAnDS 

JAR  25  +  NV 

RLD 

USA 

FAR  25 

FAA 

Fig.  1.  Airworthiness  Require'1  onto  and  AA  in  differen  t  countries 


N\  stands  for  national  variant  und  means  that  in  this  country 

additional  tc  the  basic  rule  a  special  gust-requirement  exists.  Also  in  the  case  where  the  require¬ 
ment  basis  is  similar,  different  means  u  compliance  and  interpretations  from  one  country  to 
another  might  be  valid,  as  late*-  will  ba  shown.  (Fig,  2) 
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Fig.  2,  Gust  requi  remen  t  s  in  different  regulations 
*  Floor  Level  for  Supp  1  emen  t  a  ry  PEA,  Uo  -  60  fps 


2.2  Aircraft  Modeling 

Besides  the  already  existing  differences  in  gust  requirements  as  shown  in  chapter  2.1,  different 
A/C  models  are  required  by  the  forementioned  AA. 

These  differences  are  mainly 

-  how  flexibility 

-  if  and  how  unsteady  aerodynamic  forces 
are  taken  into  account. 

Whilst  in  the  requirements  exists  a  clear  statement  how  flexibility  has  to  be  considered,  it 
is  normally  left  to  the  manufacturer  how  unsteady  lift  is  introduced. 

So  all  AA  are  in  agreement  that  in  CT-investigations 

-  dynamic  response  for 

-  the  full  flexible  A/C  taking  into  account 
steady  and  unsteady  lift 

have  to  be  considered. 

The  question  how  many  modes  and  frequencies  and  which  unsteady  lift  theory  will  be  taken 
to  have  an  optimal  representation  of  the  A/C  lies  in  the  hands  of  the  manufacturer. 

The  situation  is  different  in  the  required  discrete  gust  models.  The  QF-approach,  A/C  regarded 
as  frozen  under  deformation,  is  accepted  (see  figure  2)  from 

-  FAA  for  Pratt-Formula  and  the  (1-cos)  gust. 

-  RLD  for  their  N.V.  of  the  negative  gust. 

-  Basic  JAR  25  as  "Quasi  Pratt",  but  in 
combination  with  a  FF  discrete  (1-cos)  gust. 

For  DTG  the  CAA  and  for  the  DG  with  90%  Ude  all  other  European  authorities  require  the  FF-A/C. 

The  term  "Quasi  Pratt"  is  used  here  for  a  complete  A/C  response  calculation  of  a  (1-cos)  gust, 
but  suppressing  the  effect  of  the  flex  modes  and  frequencies. 
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Effect  of  Control  Systems  and  Control  Laws 


In  history  the  only  systems  effect  came  in  by 

-  autopilots 

-  autostabilizers 

-  yawdampers. 

Normally  for  design  the  A/C  was  conservatively  regarded  without  these  systems  and  the  effect 
had  only  to  be  demonstrated  by  a  supporting  study. 

This  situation  has  changed  since  modern  transports  have  EFCS  (FBW)  with  Control  Laws  and 
sometimes  Load  Alleviation  Systems. 

It  is  self-evident  that  AA  in  this  situation  will  request  that 

-  systems  have  to  be  introduced  in  the  A/C  modelling 

to  prevent  that  the  loads  situation  does  not  become  worse  than  for  the  conventional  A/C. 

Accepting  this,  the  interest  of  the  manufacturer  is  that 

-  systems  effects  arc  taken  into  account 
also  t.iere,  where  it  means  benefits 
for  tlie  structure 

and  that  this  is  adequately  accepted  by  AA. 

Naturally  the  problem  starts  with  the  different  opinions  represented  by  AA  and  industry  about 
the  meaning  of 

adequately. 

The  existing  requirements  do  not  cover  this  area  at  all. 

The  FAA  gives  in  its  Advisory  Circular 

AC  25.672 

"Active  Flight  Controls" 

some  advice. 

The  JAA  have  together  with  European  industry  in  the  frame  of  an  actual  A/C  design  established 
a  special  condition  to  JAR  25 


SC  -  A  2.1.1. 

together  with  interpretative  material 

1M  -  A  2.1.1. 

"Certification  Criteria  for  an  A/C  designed  with  systems  interacting  with  structural  Performance" 

The  European  paper  goes  further,  as  can  already  be  seen  by  the  heading,  and  covers  -  to 
the  opinion  of  the  author  -  the  problem  more  appropriate. 

Both  papers  are  derived  from  actual  necessities  in  the  different  countries  and  will  need  - 
having  now  more  experience  -  a  critical  review,  update  and  generalisation.  A  lot  of  conservatism 
is  in  the  early  papers,  taking  into  account  the  lack  of  experience  and  the  fear  of  the  AA 
to  erode  safety. 

In  between  the  term  "adequate"  needs  a  revised  interpretation,  not  to  destroy  the  attraction 
of  the«=e  new  systems,  which  also  brought  an  improvement  in  handling  for  the  manufacturer 
by  insisting  on  old  unjustified  conservatisms. 

Fig.  3  to  6  give  for  control  systems  of  modern  transport  A/C  in  schematic  sketches. 


Pitch  control 


FIG.  4  CONTROL  LAW  REALISATION  POR  LOADS  -  PITCH  AXIS 
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EFCS  installation 


EFCS  installation 


\ 


FIG.  6  CONTROL  LAW  REALISATION  FOR  LOADS  -  LOAD  ALLEVIATION 
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3.0  Influence  of  Requirements  on  Static  Design  Load  Levels 


The  influence  of  the  different  national  requirements  and  its  means  of  compliance  on  the  static 
design  load  level  of  a  large  transport  A/C,  which  is  designed  in  major  parts  by  gusts, 
is  demonstrated  in  the  following  figures. 

As  representative  quantities  were  chosen, 
for  Vertical  Gust 

wing  bending 
tail  bending 

for  Lateral  Gust 

-  fuselage  bending 

-  fin  bending 

Further  the  influence  of  introduction  of  control  systems,  control  laws  and  load  alleviation 
on  the  Vertical  and  Lateral  Gust  Load  Level  is  shown. 


3.1  Vertical  Gust 


The  following  gustmodels  and  MOC  are  compared: 

-  Vertical  Gust  Load  Level  due  to  European/American  requirements  (Fig.  7  and  8) 
All  system  effects  inclusive  load  alleviation  are  included. 

CT  (85  fps.  TAS)  -  reference  level 
with 

CT  (75  fps.  TAS) 

DG  +  DTG  /  100%  resp.  90%  Ude/FF 
DG  /100%  Ude  /  QF 


For  the  wing  (Fig.  7)  it  is  found  that  the  DG  and  DTG-Luvel  with  the  European  interpretation 
to  be  calculated  FF  gives  a  load  level  even  higher  than  the  .1AA  load  level  due  to  CT  (85 
fps). 

The  FAA  approach  with  the  chance  of  CT  (75  fps)  in  special  cases  and  the  DG  interpreted 
QF  gives  a  marked  lower  load  level. 

For  the  tailplane  the  highest  gust  load  level  is  found  also  for  the  JAA  gust  interpretation 
with  the  DG  interpreted  QF  next,  both  higher  than  CT. 


-  Vertical  Gust  Load  Level  due  to  Britisch  National  Variant 
(Fig.  S  and  10) 

Systems  as  before. 

DG  /  90%  Ude/FF  -  reference  level 
with 

DTG  /  shaping  law  with  max.  90%  Ude/FF 


Fig.  9  and  10  show  the  effect  of  the  gust  shaping  due  to  the  UK-National  Variant.  It  is  found 
that  at  wing  and  tailplane  depending  cn  the  spanwise  station  the  National  Variant  produces 
the  higher  load  level.  This  means  that  at  components,  where  CT  is  below  DG,  the  DTG  will 
give  the  design  case. 


-  Vertical  Gust  Load  Level,  effect  of  pitch  control  law  and  load  alleviation 
(Fig.  11  and  12) 

DG  /  90%  Ude/FF  pitch  control  law  and  LAS  incl  reference  level 
with 

DG  /  90%  Ude  FF  pitch  control  law  included 
DG  /  90%  Ude  FF  no  systems  included. 

For  the  wing  (Fig.  11)  it  is  found  that  the  introduction  of  the  pitch  control  law  gives  the 
same  level  as  the  conventional  A/C.  This  is  different  for  the  tailplane  where  the  pitch  control 
law  introduction  leads  to  an  increase  of  tailplane  loads  in  relation  to  all  "systems  introduced" 
where  the  conventions?  A/C  leads  to  an  underestimation. 
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3 . 2 Lateral  Gust 


The  following  gustmodels  and  MOC  are  compared: 

-  Lateral  Gust  Load  Level  due  to  European /American  requirements  (Fig.  13  and  14)  The  influence 
of  systems  and  control  laws  was  not  introduced. 

CT  (85  fps,  TAP)  *  reference  level 
with 

CT  (7L  fps,  TAS) 

DG  +  HTG  /  100%  rasp.  90%  Ude/FF 
DG  /  100%  Ude/QF 

The  JAR  25  requirement  basis  inclusive  the  NV  and  its  usual  M'Xl  lead  to  a  higher  load 
level  as  well  in  continuous  turbulence  as  in  discrete  gust  tha  FAR  25  and  its  usual  interpre¬ 
tations. 


-  Lateral  Gust  Load  Level  due  to  British  National  Variant 

(Fig.  15  and  16).  The  influence  of  systems  and  control  laws  was  not  introduced. 

DG  /  100%  Ude  /  FF  -  reference  level 
with 

D’t'G  /  shaping  law  with  max  90%  Ude  /  FF 

This  NV  gives  a  higher  to  equal  load  level  for  the  fuselage  and  a  lower  one  than  discrete 
for  the  fin.  In  this  design  the  CT-level  is  dominating  all  components  in  lateral  gust  which 
must  not  always  be  the  case. 


-  Lateral  Gust  Load  Level,  effect  of  systems  and  control  laws 
(Fig.  17  and  18) 

CT  ( 8G  fps,  TAS)  without  EFCS  -  reference  level 
with 

CT  )85  fps,  TAS)  with  EFCS,  100%  operating 

The  introduction  of  systems,  especially  yawdamper  function  and  lateral  control  laws  leads 
in  general  to  reductions  of  gust  loads  because  of  its  damping  effect  of  the  dutch  roll.  It 
is  obvious  that  if  this  level  is  taken  for  design,  system  failure  cases  have  to  be  investigated. 
IM-A2.1.1  defines  a  safety  factor  as  function  of  failure  probability  of  occurrence  to  calculate 
ultimate  load.  (Fig.  19  and  20)  These  figures  illustrate  the  high  importance  of  a  reliable 
system  on  structural  weight. 
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Fig.  1 9  Static  strength  factor  at  time  of  failure  /  Active  Failure 


Fig.  20  Static  strength  factor  after  failure  /  Passive  Failure 


4.0  CONCLUSIONS 


5- 


The  shown  situation  can  be  regarded  as  reprasantative  for  a  major  gust  load  designad  transport 


It  was  shown  that  more  severe  gust  requirements  or  interpretations  in  Europe  lead  to  higher 
gust  design  loads. 

It  is  doubtfull  whether  this  is  necessary  to  guarantee  an  acceptable  level  of  safety.  But 
there  is  no  doubt  that  this  situation  distortes  the  international  competition. 

The  objections  of  AAs  to  this  statement  will  be  that  also  Import-  A/C  entering  the  European 
market  will  have  to  cover  these  requirements. 

Although  this  is  true  it  cannot  satisfy  a  manufacturer  that  own  A/C  designs,  flying  in 

the  same  atmosphere  as  those  of  a  competitor,  have  to  cover  different  severe  gust  requirements, 

depending  on  the  country  of  certification. 

A  critical  review  ot  tho  total  gust  requirement  situation  is  therefore  requested,  ending  up 
in  a  harmonisation  of  the  conditions  to  be  covered,  independent  from  the  country  of  certifica¬ 
tion  . 

We  must  refrain  from  the  former  artificial  load  conditions,  coming  as  near  as  possible  to 
the  roal  physical  conditions. 

New  statistical  material  lor  better  representation  of  the  atmosphere  is  needed,  only  this 
will  guarantee  a  good  level  of  safety  and  will  at  the  same  time  give  us  the  chance  to  use 
all  tho  benefits  we  can  get  also  in  the  structures  field  from  introduction  of  modern  electronic 
systems,  to  improve  overall  economics. 

Industry,  Laboratries  and  Authorities  are  requested  to  take  over  this  task  in  the  interest 
of  all  of  us. 
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Summary  :  A  program  of  in-situ  measures  using  the  Nord  260  plane  equipped  with 
accelerometers  has  allowed  to  compare  the  predicted  and  the  measured  responses  of  the 
flexible  aircraft  to  the  turbulence.  It  shows  a  good  agreement  between  the  two  sets  of 
results  and  it  emphasizes  the  better  modeling  of  the  turbulence  using  the  isotropic  model 
rather  than  the  cylindrical  one. 


INTRODUCTION 


A  flying  program  using  the  plane  Nord  260  has  allowed  to  compare  the  experimental  and  the 
calculated  responses  of  the  flexible  airplane  to  the  low  altitude  atmospheric 
turbulence.  The  results  ot  these  comparisons  give  a  good  fitting  between  calculated  and 
measured  transfers,  spectral  densities  and  coherence,  once  the  adequate  turbulence  length 
scale  and  standart  deviation  have  been  found. 

Those  comparisons  show  also  that  the  responses  calculated  by  using  an  isotropic  model 
in  order  to  represent  the  atmospheric  turbulence  are  much  more  closer  to  the  experimental 
responses  than  those  obtained  using  a  cylindrical  model,  which  shows  the  more  realistic 
behavior  of  ths  isotropic  model. 

After  recalling  in  a  first  part  how  the  responses  are  calculated,  we  show  by  drawing 
the  curves  representing  the  transfers,  the  spectral  density  and  the  coherence,  the 
comparison  between  the  ei  periment  and  the  calculus 


I  Theorlcal  responses  to  the  turi  'ence 

1.1.  Non  stationary  forces  due  to  the  turbulence 

The  non-stat ionary  pressure  created  by  the  atmospheric  turbulence  at  point  M  and  time  t 
on  a  flexible  aircraft  is  given  by  the  relation  : 

(1)  ASP(M,t>  -  /  £<M,M',t)  «0  <M',t)  dM' 

where  G(M,M',t)  represents  the  Green  function  of  the  problem  calculated  by  the  doublet 
lattice  method  [1],  and  a(M',t)  is  the  local  angle  of  attack  induced  at  point  M'  and  time 
t  by  the  turbulence,  S  beeing  the  wing  surface 

The  generalized  non-stationary  forces  relative  to  the  displacement  field  K(M)  are  then 
given  by  : 


(2)  y  H(t)  -  Ja  ASP(M,t)H(M)dM 
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In  the  nodal  basis  of  the  plane,  the  vector  q  of  the  generalized  coordinnates  is  given 
by  the  relation 

(3)  q(t)  -  a  *  F(t) 


where  a  is  the  impulsional  response  and  F  is  the  vector  of  the  generalized  forces 
components.  Ne  denote  A  -  A  the  admittance  mattrix. 

Knowing  the  vector  h(M)  of  the  displacements  (h  (M> ) 

1  i 

of  the  eigen  inodes  at  point  M,  we  can  deduce  the  vertical  displacement  Z  at  that  poin- 
using  the  relation 


(4)  Z  (M,  t)  -  ^  h  (M)  q^t) 

In  order  to  achieve  the  calculus  it  is  necessary  to  know  the  expression  of  the  local 
angle  of  attack  <X(M,t)  in  function  of  the  turbulence.  This  expression  will  be  different 
according  to  the  model  chosen  to  represent  the  atmospheric  turbulence. 

1.2.  Cylindrical  turbulence 

The  cylindrical  model  is  a  simplified  model  of  the  atmospheric  turbulence  which  suppose 
that  it  is  constant  spanwise  and  which  is  entirely  defined  as  soon  as  its  value  at  a 
single  point  M  ^  is  known. 

More  precisely,  for  an  airplane  flying  along  the  x  direction  at  speed  v,  the  local 
angle  of  attack  <x(M,t)  induced  by  the  turbulent  field  "VffM.t)  is  given  by  : 

(5)  a  (M,  t)  =  Iir  <M,t)  /  v  -Ur  (M  ,t  -  <*.-  *) /v)  /  v 

where  M  (xfl)  is  the  point  where  the  turbulence  is  measured  and  which  is  located  at  the 
nose  of  the  plane. 

Denoting  A^(M ,<d)  the  Fourier  transform  of  the  non-stationary  pressure  AS(M,t)  at  print  M, 
we  get  easily  : 

(6)  A?(M,o>)  = -u r  <M,to)  /  V  Js^.(M,M',<n)exp(-jO)(x0-Xi(  )  /  V  )  d  M' 

Therefore  we  see  that  the  pressure  is  obtained  by  a  linear  filtering  of  the  cylindrical 
turbulence,  the  frequency  response  beeing  given  by  the  previous  integral. 

In  the  same  way  the  non-staticnary  forces  and  therefore  also  the  vertical  displacement 
Z(M,t)  are  obtained  by  a  linear  filtering  of  the  cylindrical  turbulence  and  we  consider 
the  frequency  response  that  we  shall  call  transfer  at  point  M  by  : 

(7)  T  <M,0»  =  Z(M,0»  /W  (M  ,01)  =  1/V  h1  A  F 

Vi  ° 

We  can  then  construct  : 

<I>  (M,w;  -  I  T  (M,0))  1*0  (to) 
zz  »z  Mrur 

and  the  coherence  : 

y*(M,to)  =  |0  <M,a>)  l*/0  (to) O <M, <a) 

V2  WV  ZZ 

1.3.  Isotropic  turbulence 


The  vertial  atmospheric  turbulence  is  modeled  by  a  second  order  stationary  continuous 
and  zero  mean  gaussian  process  wich  is  isotropic  and  homogeneous.  In  this  context,  the 
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local  angle  of  attack  is  given  by  the  relation 

(8)  a(M, t)  -v  (M,  t)  /  V  , 

and  the  non-statlonary  pressures  and  the  responses  of  the  aircraft  are  not  given  any  more 
by  a  linear  filtering  of  the  turbulence  at  the  single  point  M„  . 

In  order  to  get  the  transfer,  we  are  going  to  construct  directly  the  spectral  density 
cd)  of  the  response  Z 

at  point  M  and  the  interspectral  density  between  the  turbulence  measured  at  the  reference 
point  M  and  the  response  Z  at  point  M,  <I>  (M„,M, a>). 

To  calculate  the  generalized  forces,  the  surface  S  is  discretized.  we  denote  , 

G(t)  «  (G^(t))  ■  ($(  Mi,M,t)>i  the  mattrix  of  the  discretized  kernel  ^.(t)  , 

H  •  (  h^(K))  the  mattrix  which  columns  contain  the  displacement  of  the  mode  i  at  point  M 
of  the  lattice,  W(t)  the  vector  which  componants  are  the  value  of  the  turbulence  at  time 
t  and  at  the  points  of  the  lattice  M  . 

The  vertical  displacement  at  point  M,  Z(M, t)  is  therefore  given  by  the  relation 

(9)  Z(M,t)  -  hT(M)  A  »  (  H  G  *  W/V  )  (t) 

The  process  Z  is  thus  obtained  by  the  composition  of  two  linear  filtering  of 
process  W/V  ,  the  first  impulse  response  beeing  HG  <t) ,  the  second  one  beeing  h 
general  theory  of  second  order  processes  ( 4 ]  show  then  that  the  process  z  is  a 
order,  zero  mean  continuous  gaussian  process  which  is  stationary  ,  the  spectal 
which  beeing  given  by  the  formula  : 

(10)  0  (M,cd)  -  hT(M)  A  (01)  H  G  (00)  S  (CO)  G*(CD)  H*  A*(fl>)  h  (M) 

ZZ  H/V 

The  quantity  3  ^(cd)  the  mafrix  spectral  measure 
of  the  process  W/V  which  elements  are  defined  by  : 

(11)  S  (co)  =  [  S  (M  ,  M  ,a»  ] 

w/v'  W/v  1  j  l,j 

where  ^  (M  ,M  ,co)  represents  the  transverse  spectral  measure  of  the  process  w/v  at  the  two 

points  M  and  M  . 

i  ) 

In  the  same  way,  the  cross-spectral  density  between  the  response  Z  and  the  turbulence 

W/v  at  the  measure  point  M  is  given  by  : 

o 

(12)  *  (M  ,M,0>)  -  hT(M)  A  (CD)  H  G  (CD)  S 

9Z  o  W/V 

0 

where 

(13)  S  (CD)  -  [  S  (M  , M  ,CD)  ] 

W  /V  W/V  0  J  j 

0 

represents  the  column  of  the  transverse  spectral  measure  at  point  M  . 

In  order  to  be  able  to  define  completly  expression  (10)  and  (12),  we  still  have  to 
write  down  tne  transverse  spectral  measure  S  (M',M,®>)  between  two  joints  M  and  M' . 

ID/v 

This  quantity  has  been  calculated  in  reference  [3)  : 

(14)  S  (M,M',CD)  -  &xp(-i(x-X')?C(CD,T|)0  (CD)  ■  O  (CD)®  (0») 

w/v  m  w/v  ww 

(15)  C(cd,ti)  =  2/r<p+i/2)  ((tj/2)p*1/jk  (T1)  -  2  B(tj/2)p*j/j  Kp"1/j(tj)) 

P+1/2 

D  =  p/a  V(ma);/vJ+l 
p  =  V  <y-y' )2+  (*-*' )J 
a  -  iWx  r (p)  /  r<p+i/2) 

B-  (l-(l+(a(D)W)  <<z-z')/p)’>  (l+2(p+l)  (a<o)J/vVl 


the 

T  (M)  A(t)  .  The 
second 
moasure  of 


I" 


L  beeing  the  turbulence  length  scale,  K  the  Bessel  function  of  order  p. 

p 

*  (01)  Is  the  spectral  density  of  the  vertical  turbulence  : 

(16)  QJ.U)  -  o*  l/v2k  (1+2 (p+l)  (ato)J/V2)  <l+(aa))2/v2  )  'P"3'3 
The  Kerman's  model  is  obtained  when  p  »  1/3  . 

Ne  can  then  construct  a  pseudo  transfer  T  between  the  response  Z 

ii 

and  the  process  of  the  vertical  turbulence  : 

(17)  ®  (M,W)  -  I  T  (M,W)  |2®  (0>) 

ZZ  Et  WN 

where 

(18)  |TJM,0»)|2  -  h^(M)  A(0»  H  G<0>)  O^<0»  G  (01)  H  A  (0))  h(M) 

In  the  same  way,  using  relation  (12),  we  can  construct  the  pseudo  transfer  T  : 

(19)  <1*  (M,  n>)  -  T  (M,a>)  ®  (oo) 

NZ  HZ  HH 

and  the  pseudo  coherence 

V  (M,  (0)  =  |®  <M,0>)  I2/®  (<i))®  (M,  fli>) 

1  */v  i  w  zz 

1.4.  Turbulence  length  scale 

The  quantities  C  (ovn)  and  ®  (o))  depend  of  the  turbulence 

WW  2 

length  scale  L  and  of  the  standart  deviation  a  .  . 

If  we  want  to  compare  the  theories!  plane's  responses  to  the  measured  ones,  it  is 
necessary  to  estimate  correctly  these  two  parameters  which  characterize  the  turbulence 
state  during  the  in  flight  measures.  In  order  to  do  that,  using  the  measured  spectral 
density,  we  do  the  ratio  of  two  values  of  the  spectum  (16)  for  two  different 
frequencies  and  we  obtain  a  relation  which  depends  only  of  the  length  scale  L.  We  get  L 
using  an  numerical  iterative  scheme. 

Then  we  get  easily  the  standart  deviation.  This  method  is  however  limited  by  the 
scattering  of  the  experimental  spectrum  values  but  it  gives  a  rough  estimate  that  it  is 
possible  to  improve. 


II  In  flight  measures  for  the  Nord  260. 

II. 1.  Modal  Basis 

The  modal  basis  contain  two  rigid  modes  of  pitching  and  plunging  and  twelve  flexible 
modes  obtained  during  ground  vibrations  measures  [2]  and  whicn  characteristics  are  given 
below  : 


frequencies  Hz 

1  4.52  6.32  8.03  9.5 

10.41 

11.95 

13.30 

15.59  16.17 

23.3  27.8 

generalized  mesa 

6.4  201  659  692  6.7 

135 

345 

431 

1.2 

1070 

88.2  256 

structural  damping 

.35  .009  .014  .011  .034 

.027 

.016 

.017 

.025 

.022 

.015  .018 

Ill  Comparison  measure  and  theory 


We  use  a  length  scale  of  75  m  for  the  turbulence  which  give  a  value  of  .78  for  the 
standart  deviation. 

igure  1  shows  the  comparison  of  the  experimental  and  theorical  spectrum.  The  inertial 
area  is  well  represented. 

From  now  on  we  shall  only  look  at  the  acceleration  Z  (M, t)  of  the  vertical 
displacement.  Figures  2  and  3  show  for  two  different  accelerometers  how  the 
theorical  response  compare  when  either  the  isotropic  or  the  cylindrical  model  is  used.  We 
can  see  a  rather  big  difference  for  the  two  results,  especially  in  the  flexible  modes 
area.  Such  a  difference  had  not  been  found  for  ths  responses  of  the  Mirage  III,  this  last 
airplane  having  a  much  smaller  aspect  ratio  than  the  Nord  260  which  can  be  considered  as 
a  large  aspect  ratio  plane. 

The  5  last  figures  show  the  comparison  between  the  theorical  response  using  the 
Isotropic  model  and  the  experimental  response  for  5  different  accelerometers.  We  see  the 
good  agreement  of  the  results  except  ror  the  figura  7  which  involve  accelerometers 
located  at  the  wings's  tip.  This  is  due  maybe  to  bad  measures  at  these  points. 

The  results  for  the  transfers  are  not  much  modified  when  the  value  of  the  turbulence 
length  scale  varies.  This  is  not  the  case  for  the  power  spectral  densities. 

Of  course  the  agreement  is  not  perfect  between  the  experimental  and  theorical  results  : 
the  turbulence  is  already  anisotropic  at  the  altitude  of  loO  m  where  the  measures  were 
made  and  there  is  always  some  noise  wich  add  up  to  the  turbulence  excitation. 
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1.0  nnwotcriPH 


In  tha  pact  various  gust  load  fcoailac  haw  ban  dewlcpad  for  tha  calculation  of  design  gust  loads  on 
aircraft.  Since  It  woe  first  abltehad  In  1954  (ftaf  1),  ths  al Isolation  factor  nroach  of  Pratt  and 
Miter  has  galnsd  alaost  uni  van  al  aocaptanoa  and  for  aany  yaara  has  bsan  a  faalllar  part  of  ths 
airworthiness  raguirassnta  for  loth  civil  and  military  aircraft. 

lbs  original  ocnoapt  of  ths  Trull  fowls*  waa  to  paadlet  ths  paak  accelerations  due  to  dlaczata  gusts 
an  a  giwn  aircraft  for  flight  throu^i  a  diacrats  gust  or  tin  sans  shape  arid  aapl  ttuda.  ihua  tha 
darlvad  guat  valocity  is  not  so  mtch  an  absolutes  physical  quantity,  but  la  rathar  aora  a  gust-load 
tmfar  factor  dafinad  within  ths  tarns  of  tha  focaula.  As  such  ths  satehod  la  tost  accurate  whan 
rastrlctead  to  usa  on  aircraft  with  vary  stellar  charactearlatica.  baas  confidence  mat  ba  sttaohad  to 
prsdlcfasd  gust  loads  an  aircraft  which  ara  unoonwntlonal  whan  ocsparad  to  tha  data  oollacting  aircraft. 

nrtunsbaly,  as  In  sany  other  f  la  Ids  of  aarcnautloal  solanca,  tha  anthoda  of  calculating  dasign  gust 
loads  haw  baan  constantly  l^prowd  onr  tha  past  30  yaara.  national  and  aocurata  analysis  prcoaduraa 
ara  now  ggsn  wongat  all  tha  aajor  aircraft  sanufacturars.  Thsaa  haw  allowad  all  tha  Irportant 
taaturas  of  nsw  aircraft  to  ba  sodstlsd  and  thus  ham  anaurad  that  diffaranow  In  guat  raaponea 
charactesrlstlos  ara  Silly  accounted  for  In  ths  dssl9«  prooass.  This  Is  Silly  In  sooord  with  s 
fwdwontal  principle  of  airworthinaas  rwquirassnts  which  dsaanda  that  load  lntansitlas  and 
distributions  should  cloaaly  raprsssni  actual  oondltlcna  or  alas  thsy  auat  ba  shown  bo  ba  oonsarvatlvs. 
In  ganartl,  howavar,  only  ana  sat  of  deei<*i  gust  wlodtlas  la  prescribed  and  thasa  ara  has  ad  upon 
aqulvalsnt  valuas  darlvad  using  tha  ahqpla  focaula. 

This  would  ba  raaaonabla  If  all  ths  alawants  of  tha  aodarn  guat  analysis  wars  only  applioahla  bo  sodsm 
aircraft  and  if  tha  affacta  an  ths  previous  ganarstlon  of  aircraft  wavs  inslpilflaant.  Bowever,  this 
is  patently  not  tha  oaaa.  Tha  allavlaticn  factor  focaula  approach  laana  heavily  an  aany  simplifying 
aaaiaqiMtea  which  wan  introduced  puxaly  to  facilitate  ths  aolution  of  tha  aquations  of  notion.  That  is 
not  to  say  that  all  tha  alswnta  of  the  aquations  of  action  which  wars  alternated  by  Pratt  and  Miter 
and  other  rase  archers  wars  naglihila  but  rathar  that  thair  Inclusion  would  haw  lad  to  an  overly  arduous 
analysis  task. 

Tbs  degree  to  which  currant  sew  apt  Iran  oanoaming  tha  NodaUing  of  aircraft  response  to  diacrata 
gusts  oould  haw  affactad  earlier  oonoiusions  vpon  gust  static  tics  darlvad  fro  saaaursd  aooateratlon 
data  la  lha  subject  of  this  paper.  Tha  work  described  waa  dona  by  British  Aerospace,  Mybrldga  inter  a 
rasa  arch  contract  sponsored  by  tha  Civil  Aviation  Authority  of  tha  united  Kingdo. 

Tha  opinions  expressed  in  this  paper  ara  aolaly  those  of  tha  Author. 


2.°  nwanannat  or  vm  anas  cr  taxcmc  mptwbb  grew  pwoicnp  kccgaBEBU 

The  dasign  gust  velocities  in  currant  gust  load  rag  Hr  wants  are  baaed,  for  tha  Boat  part,  an 
statistical  data  collected  on  U.S.  trriwport  aircraft  prior  bo  19S0  and  supported  by  further  data 
ool looted  an  European  transport  aircraft  prior  to  19*0.  Guat  valocity  pratebllittes  wars  derived  frtai 
eaoaadaroa  aounte  of  centre  of1  gravity  acceleration  using  a  staple  guat  allavlaticn  factor  approach 
based  upon  tha  aircraft  aws  pazwstar.  Gust  velocities  darlvad  In  this  way  oat  ba  confidently  used 
for  tesign,  providing  that  ths  aircraft  under  cone  Iteration  has  slat  laor  rapcraa  charactaristics  to  tha 
data  ooliactitg  aircraft.  Mhara  this  is  not  ths  csss  a  aora  realistic  picture  of  tha  distribution  of 
true  gust  valocity  should  ba  sought.  Tteally,  dwlpt  esthete  and  daal^t  i  a  j  ill— il  ■  should  ba 
ocaplaasntaryi  together  producing  tha  required  safety  objsetiva  without  unaja  econoalc  penalty. 

Gust  load  analysis  aethote  how  drwlopad  steadily  mar  raoant  yaara  to  tep  in  stop  with  dtuelcpaant a 
in  aircraft  technology.  The  affacta  of  f  legibility.  wing  wasp,  high  apaad  aarcdynwlos,  transport  delay*, 
lntarfsranoa  bstwaan  lifting  surfacss,  ocnflguratlon  chwgw  (flaps,  slats  spoilers  ate)  and  autopilot 
ara  now  routinely  included  in  tha  teaiqi  analysis  Saner  appropriate.  Soa  of  thaw  faaturaa  ara  also 
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wlwt  to  tin  data  oalkouni  aircraft,  tout  un  aaalutad  f»  tha  data  reduction  prewar  to 
faellltata  a  solution  of  tha  aquattana  of  notion.  In  ortar  to  datanlna  tha  af  facta  that  thaaa  features 
aould  hme  had  on  the  gun  statistics  data,  thanfom,  it  waa  naoaaaary  firat  to  untaretand  how  each 
elan ant  of  a  aim  malyeia  aould  have  lnfluanoad  tha  conclusions  of  aarllar  taaaarch  atudlaa,  had  adi 
analytical  toola  bam  mailable  at  tha  tlaa. 

Fa— 1  nation  of  fefatanw  1  shows  nlna  dlffazant  aircraft  typaa  mm  uaad  to  aollact  V-Q  data  during  tha 
parlod  U1WW0.  ot  thaaa  to  aircraft  which  tlfai  pcreal nantly  haw*  barn  aalactad  4.1  worthy  of 
ranawad  invest  iqat  Ion.  -  Type  K  and  Type  j.  Thaaa  aircraft  warn  aalactad  aa  being  representative  of  tha 
claaa  ef  aircraft  tar  which  tha  Pratt  tanatla  la  apaitad  to  bo  amt  affactlwa. 

Thay  ara  alao  aircraft  typaa  tor  vdilch  reasonable  aatlaataa  of  wing  aaaa  md  stiffnaoa  data 
axlat  in  publlahad  liiarature  (hat  2) .  (Ming  thia  data  it  was  poaalbla  to  aat  up  aathaaatleal 
mdBla  which  would  give  a  reasonable  anginaarlng  rapraaantatlai  of  tha  91st  raapon— 
oharactarlatloa  of  aaoh  typo  In  aouordmoa  with  present  thaary . 

larly  >rtfaa>  aowlarcaatar  data  waa  coUamad  on  a  larga  note  of  different  aircraft 
typaa.  Of  thaaa  cm  (Typa  V)  waa  aalactad  m  fairly  typical,  it  waa  alao  advantageous  to  stuty 
m  aircraft  on  which  reliable  dMigi  dam  waa  ami  labia.  All  «1  Inwart  a  datallod  and  aocurata 
aat  heretical  aodal  of  tha  aircraft  to  ho  (emulated  and  chachad  agalnat  tall  aoala  taat  raaulta. 
Publlahad  data  on  thia  alrcnlt  Oaf  3)  haa  m  addad  raflnaamt  of  a  haight  bmd  bcaatatotm, 
requiring  Investigation  of  guat  load  factors  at  graatar  altltudao  thm  war*  naadad  tar  tha  U.S. 
aircraft. 


3.0  MnOWPT  THDtnClL  MPOBd 

fet  hmal  Inal  aodala  of  tha  aircraft  infer  review  ware  dtvalcpsd  using  tha  latoat  analytical 
techniques  froa  previously  puhliafed  aaaa  md  atiffnaaa  data.  Rigid  aarodyn— lc  load  distributions 
ware  calculated  tar  each  lifting  aurfaca  using  a  ooablnatlon  of  vortax  lattloa  md  mlandar  body 
theory,  in  areas  whara  detailed  Intamatlon  was  not  aval  labia  anginaarlng  judgment  was  amreiaod 
to  provide  a  wafting  approulaatlan  of  tha  nlaslng  data,  Owing  to  thaaa  asauaptlona  it  camot  ba 
clalaad  that  tha  final  thaoratloal  aodal  gam  an  sweet  rapraamtatlon  of  each  aircraft.  Itowavar, 
it  should  ha  is  mart  that  tha  aodala  do  poaaaaa  tha  inherent  71st  raaponaa  character  la  ties  rapraa- 
antati—  of  the  clam  of  aircraft  untar  study.  *•  ouch  each  oat  ba  asm  as  lrdicatlva  of  the 
relative  tranta  in  guat  factors  md  gust  loads  to  ha  associated  with  each  change  In  aodalllng 
technique. 


With  tha  aodal  of  tha  Typa  V  aircraft,  of  course,  this  la  not  tha  case.  This  methane ltloal 
rapraamtatlon  waa  baaed  upon  validated  design  data,  with  tha  structural  nodes  affirmed  by  the 
raaulta  of  tall  aoala  stiftneaa  basts,  and  ground  md  flight  vibration  tests,  with  aorafyniadR 
data  supported  by  wind  twcal  taat  aaaauraaanta. 

a  ooapariaon  of  tha  aajor  design  faaturaa  of  each  typa  of  aircraft  studied  is  given  in  Table  1 . 

In  each  oam  tha  aircraft  raaponaa  aquations  ten  established  In  a  way  which  persitted  an 
individual  variation  In  aaaa  aspect  of  ■odslling  technique.  At  each  stags  tit?  reaponsra  to  a 
range  of"  1  -cosine”  stuped  gusts  wan  calculated  to  datsmlna  the  change  in  incnmntal 
aowlaratlon  at  tha  oantra-of-gravity . 

Tha  aodalllng  aaawptlons  inherent  in  tha  Pratt  fonula  approach  to  raaponaa  analysis  an 
sums rlaad  In  Table  2,  together  with  brief  da  tails  of  tha  laprovaaanta  mini  I  ail  In  currant 
stata-of-tha  art  aotalUng  prooatans.  Analysis  of  som  typical  raaulta  obtained  for  each  claaa 
of  aaroplma  an  aiaaarlaad  In  Flguraa  1-3.  Them  figures  show  hew  tie  predicted  incraamtal 
c.g.  acceleration  cm  vary  tfepandirtg  upon  tha  aaeaqrt  Inna  seta  concerning  each  physical 
characteristic  of  tha  aircraft  lndutad  In  tha  aquations  of  notion. 

Bach  block  In  tha  diagrm  represents  tha  —slant  predicted  oantn  of  gravity  aowlaratlon  In 
response  to  a  ’1-ooalna"  shaped  gust,  feeding  free  left  to  right  rapraamta  tha  Introduction 
of  a  aalactad  igrtwaait  In  aotalUng  tachnlqua.  Tie  far  right  of  tha  figure  faithfully 
represents  tha  ccaplata  design  procedure  In  aooordanw  with  current  UR  practice.  Ir  each  case 
tha  dates  value  shown  la  that  wnlch  would  bn  obtained  free  tha  basic  rigid  Pratt  forsula 
solution.  It  can  ba  asm  that  tha  general  trend  la  towards  a  higher  reaper—  as  tha  inathaaatical 
aodal  la  preogrsesively  refined. 

Tha  final  result  confirm  that  tha  slapla  alleviation  factor  approach  will  seriously  under¬ 
set!— te  tha  true  acceleration  which  would  ba  raised  by  ary  dlacrata  guat  of  a  given  nagnltuda. 

Oonverasly,  derived  guat  va  loci  tlaa  baaed  upon  thaaa  slapla  jodnla  will.  In  tum,  ba  ovar- 
aatlmbw  of  tha  true  guat  velocity. 

4.0  aruaenan  k»  aucrmt  deskm  ipmb 

In  addition  to  tha  calculation  of  aircraft  aowUratloia  In  tha  various  oaaae,  eetlaat—  of  tha 

veluaa  of  so—  Uportmt  load  par —alarm  were  alao  derived  so  that  tha  affacta  of  tha  various 

aodalllng  neflnamnta  upon  daslpt  loads  could  ba  chachad.  The  overall  trends  for  Wing  Root 

Banding  Want  for  each  clam  of  aircraft  Investigated  ara  aiwmi  1—1  in  Figures  4-6.  Tha 

remits  oonflre  that  sedam  analysis  tachnlqu—  are  quits  ocn— rvatlva  whan  oaapamd  with  tha 

Pratt  foraula  approach  and  thaw  have  lad  to  a  gradual  but  slptificmt  lncraam  in  tha  j 

severity  of  tha  dasipi  criterion.  Thaw  hlddan  strength  lncreaaw  aat  ba  aaaa— art  In  tha  light  j 

of  currant  safety  objectives  too  ensure  that  tha  intended  airworthiness  targets  are  aat 

without  undua  eoonwlc  penalty.  1 
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5.0  RHAfiggaspff  CP  BBT  HgffiO  DlgglBUTIOM 

Ota  maximum  value  of  oentre-of-grevlty  acceleration  that  results  frcm  an  encounter  with  a  vertical 
gust  is  largely  dependant  igpon  waighc,  speed,  altitude  and  oentrs-of -gravity  position.  The  first 
time  of  those  paraatsrs  ears  fully  accounted  f-  r  by  Pratt  aid  Mslksr  In  their  assasMsnt  of 
operational  gust  statistics  (Mv  1) .  As  only  the  heave  degree  of-trasdcm  eas  ooneldered  at  that 
tine,  the  position  of  the  oantre-of-gravlty  wee  lees  tar  lei-  This  recent  work  has  also  confirmed 
that  the  aaverity  of  the  laspaoa  to  a  gust  la  also  very  dependent  upon  the  gust  gradient  distance. 
This  la  domcnstrated  In  Figures  7-9,  tdiich  present  calculated  responses  of  each  class  of  aircraft 
used  In  the  current  study  to  a  range  of  discrete  gusts  of  varying  wavelength,  lhe  tuning  effect  Is 
clearly  evident  even  for  these  relatively  rigid  aircraft  where  the  dymedc  effects  of  wing  flexibility 
would  be  expected  to  be  such  lower  than  on  a  modern  high  aspect  ratio  wing,  particularly  one  with 
pylon  mounted  engines. 

As  all  the  early  gust  leeponee  data  we  collected  using  standard  V-G  recorders,  no  Information  on 
gradient  distance  la  te  allabls.  Early  raeaarcharm  circumnavigated  this  problem  by  arbitrarily 
assuming  all  gusts  to  la  of  equal  wavelength  for  each  aircraft  type,  although  varying  from  type  to 
type  owing  to  the  flwd  chord  length  ratio.  The  sensitivity  of  tie  data  collecting  aircraft  to  gusts 
of  different  wavelengths  can  be  Introduced  Into  the  me  earn  smut  procedure,  however,  providing  that  a 
probability  of  occumanos  cm  be  established  for  each  gust  length. 

Alternatively,  research  by  Janes  into  the  statistical  nature  of  discrete  gusts  has  shown  (Refs  4-5) 
that  f sullies  of  rnqp  gusts,  whose  amplitudes  vary  as  "wavelength  to  the  pewer  one-third",  have  an 
equal  probability  of  occurrence.  This  physical  characteristic  has  been  demonstrated  in  turbulence 
analysis  studies  carried  out  by  British  Aerospace  at  Meytrldge  (Ref  6)  for  single  rape  and 
ooebinetiora  of  reaps  with  e  Booth  "1-coelna"  profile.  Adaption  of  this  feature  of  equiprobabllity 
to  more  faadllar  eywmetrlc  discrete  gust  (Figure  10)  will  allow  a  family  of  “1 -cosine"  gusts  to  be 
defined,  each  iiwel  m  of  which  will  be  as  likely  to  occur  as  any  other.  Furthermore  If  the  25  chord  is 
choaan  as  the  datixn  for  definition  of  the  amplitude  healing  law  aa  in  Figure  11,  then  a  tuned  gust  of 
equal  probability  to  the  25  chord  gust  can  be  found.  Thus  factors  for  conversion  of  acceleration 
to  derived  gist  velocity  can  be  established  which  incorporate  the  aonc^rt  of  variable  gust  lengths, 
whilst  maintaining  the  sen  level  of  probability  inherent  in  the  earlier  analyses  of  Pratt,  Walker, 
Bullen  and  others. 

The  effects  of  the  application  of  the  discrete  gust  power  law  on  the  peek  acceleration  can  be  seen 
In  Figures  12-14  which  provide  a  direct  ooepariaon  with  the  constant  amplitude  gust  responses  of 
Figures  7-9. 


The  development  of  acceleration  to  gust  velocity  conversion  factors  appropriate  to  the  latest 
modelling  assumptions,  and  a  comparison  with  the  Pratt  datum  result  Is  summarised  In  Figures  (15-17) . 
These  factors  are  baaed  upon  representative  flexible  aircraft  models  which  for  the  aircraft  type  v 
will  give  the  true  correspondence  between  gust  velocity  and  acceleration  at  the  moan  cruising 
oondltlms.  For  the  MS  aircraft  the  factors  represent  a  reasonable  estimate  of  the  likely  correspond¬ 
ence  due  to  uncertainties  in  the  model  data. 

The  prime  use  of  the  revised  gust  factors  described  above  has  been  In  the  evaluation  of  more  realistic 
derived  gust  velocities  appropriate  to  a  large  part  of  the  American  V-G  data  collected  prior  to  1*)50 
and  the  part  of  the  V-G-H  data  measured  on  European  aircraft  between  1953-1958.  Calculated  valu  a  of 
derived  gust  velocity  had  previously  been  fitted  with  theoretical  extreme  value  dlstrllMtlons  (Ref  7) 
in  order  to  smooth  out  Irregularities  in  the  data  and  to  provide  a  cents  tent  basis  for  extrapolation. 
Far  this  Investigation  therefore  it  was  necessary  only  to  apply  the  relevant  gust  conversion  factors 
to  che  given  frequency  distributions.  This  has  a  powerful  effect  upon  the  gust  probabilities  as  stown 
In  Figures  18-19. 

It  Is  emphasised  that  the  new  result  is  a  product  only  of  revised  assumptions  concerning  the  response 
of  the  aircraft  to  the  gust.  These  asmmptions  were  made  In  the  light  of  accepted  modern  techniques 
of  mathmstlcal  modelling  and  of  recent  theories  of  the  nature  of  the  atmosphere.  Other  aaswptlcns 
pertinent  to  the  operation  of  the  data  collecting  aircraft,  much  as  estimates  of  mean  operating 
weights,  speeds  and  altitudes  have  not  formed  part  of  this  investigation  and  have  been  accepted  aa 
soundly  based. 

6.0  AIRWORTHOCSS  OCMBIDBRXnCMS 

British  Civi  l  Airworthiness  Requirements  were  revised  in  Decaiber  1964  in  an  attmpt  to  anticipate 
the  needs  of  turbine  engined  operations.  The  gust  probabilities  mummed  at  that  time  were  extra¬ 
polated  from  pre-1958  operational  data  and  the  target  limit  gust  velocity  wee  set  at  a  value  that 
would  be  net  once  In  30,000  hours.  It  wee  assumed  that  the  gust  occurrence  rate  would  not  vary  with 
tine  so  that  experience  frcm  30,000  hours  of  turbine  powered  operation  could  be  equated  directly  to 
30,000  hours  of  oiston  powered  operations.  The  minimum  requirement  of  SOft/sec  EAS  at  20,000  feet 
reducing  to  25ft/ as'  at  50,000  feet  therefore  reflected  the  maximal  gust  experienced  In  30,000  hours 
at  250  nph  assuming  a  dynamic  stress  factor  of  1.0.  The  constant  gust  velocity  below  20,000  Met 
represented  a  relative  increased  probability  of  meeting  the  design  guts  but  this  wes  balanced  by  the 
decreased  exposure  tinea  associated  with  climb  and  descent. 

From  1964  the  BCAR  discrete  gust  design  velocities  haw  reflected  those  contained  within  U.S.  Federal 
Aviation  Regulations.  That  is  not  to  say  that  the  requirements  themselves  are  Identical  since  the 
FAR  25  discrete  gust  requirements  specifies  a  fixed  gradient  distance  of  12.5  wing  mean  chords.  In 
the  United  Kingdom  there  Is  concern  over  the  need  to  take  account  of  aircraft  dynamic  response  in  the 
calculation  of  design  gust  loads.  To  this  end  the  concept  of  varying  gust  gradimt  distance  to  find  the 
most  critical  responses  on  different  parts  of  the  aircraft  has  been  introduced.  Although  far  many 
aircraft  the  12.5  chord  gradient  le  close  to  being  the  most  critical  far  response  quantities  such  as 
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wing  root  tending  nant,  mast  aircraft  have  laportant  raapcnsas  affecting  atruotural  load*  which 
tin  to  otter  gradient  dlatanoaa. 

Thus,  while,  from  a  point  of  vlar  of  structure  weight  tte  12.5  chord  guat  ^radiant  la  naar  to  telng 
the  aoat  lapcafetet,  thla  la  not  tte  aaaa  tor  all  aircraft  nor  la  It  tte  eaaa  for  all  parta  on  *v 
aircraft,  fetus  dawslapseata.  aueh  aa  tte  application  of  active  control*  or  otter  advenoa 
acnflguMtlona  own  ba  aoca  sensitive  to  discrete  guat  affects  than  past  daslgna  and  It  tteesfora 
baocaaa  Mens  feettivi  to  oonalter  aora  than  ana  gradlant  dlstanoa. 

Aa  ftnlaaly  item  gust  velocities  developed  tram  sUpla  rs^imee  nadala  axa  severs  when  uaad  with 
a  felly  rational  ^nite  analysis.  ite  CM  tea  rarmyilaail  thla  fact  and  in  lta  National  variant  to 
JAR  25  tea  introduced  a  reduoad  guat  velocity  for  uaa  In  oonjvmctlcn  with  a  dyrsmtr-  analysis.  It 
can  be  aaan  In  Figure  30  that  tte  raductlon  rsoopUsae  two  l^octant  considerations.  Firstly,  that 
■odarn  analytical  aathoda  would  aigporl  lower  pmhahdlltlaa  of  "true"  gust  velocity  than  tte 
alleviation  lector  jgproach.  Secondly,  that  tte  reliability  of  nodam  jet  tranepait  tee  led  to 
greater  expectations  of  cpezatlonal  safety  tten  ware  thought  poaalbla  wtan  tte  turhlna  anginad  aircraft 
wan  first  lntroduoad.  An  airworthiness  objective  of  one  encounter  with  tte  limit  guat  In  50,000 
flying  hours  la  new  thcwrfit  to  ha  ^preprints. 

7.0  CPNCUUBICNB 

1)  In  tte  past  simplified  sodsla  of  aircraft  used  to  assess  operational  guat  statistics  have  led  to 
conservative  eat  isalne  of  derived  gust  sjreedaraaa. 

2)  Itodarn  refinements  in  aircraft  modelling  techniques  have  gradually  introduced  ocneervatien  in  tte 
process  of  calculating  gust  loads.  Oust  statistics  navi  awed  In  tte  light  of  these  sodern  analytical 
methods  support  tte  CAA  view  that  guat  velocities  developed  for  use  with  slaple  rigid  aircraft 
models  are  too  aware  for  use  with  a  aodsm  dynamic  analysis. 

3)  Evan  in  the  light  of  i^ooved  safety  targets,  a  10  par  cent  reduction  in  deslipi  gust  velocity 
cm  be  readily  justified,  further  reductions  may  be  justified  on  tte  basis  of  mission  analysis 
considerations,  or  by  investogetion  of  more  reoit  acceleration  statistics  collected  fay  tte 
current  generation  of  transport  aircraft.  In  the  latter  oaaa  it  will  be  essential  to  aooount  far 
all  relevant  feature#  of  the  subject  aircraft  in  the  derivation  of  gust  velocities  so  aa  to  obtain 
a  true  picture  of  tte  guat  statistics. 
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.  14,  Abstract 

A 

The  ^rge-scale  use  of  flight  recorders  by  commercial  airlines,  coupled  with  the  enhanced  quality  of 
results  offered  by  modem  computer-based  reduction  processes  makes  it  possible  to  broaden 
knowledge  of  the  phenomenon  of  atmospheric  turbulence.  At  the  same  lime,  new  methods  for 
predicting  the  response  of  flexible  aircraft  to  turbulence  are  being  projosed,  and  novel  gust 
alleviation  systems  are  being  designed  and  tested.  The  Structures  and  Materials  Panel  held  two 
Workshops  on  this  topic,  one  in  October  1 986  and  the  other  in  October  1987. 

rresentations  at  the  first  Workshop,  which  concentratedotahgeasurements  of  Turbulence  and  Data 
Collection,  are  recorded  in  the  main  report  —  AGARD-R-734?\his  volume  records  the  presenta¬ 
tions  made  at  the  sccondWorkshop^which  was  held,  under  the  chairmanship  of  , 

Mr  R.F.O’Connell.'rb  dlscuss ^dethods  of  Analysis  and ^gulations.  T^p^s 
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